On the Cauchy problem for gravity water waves 
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Abstract. We are interested in the system of gravity water waves equations 
without surface tension. Our purpose is to study the optimal regularity thresholds 
for the initial conditions. In terms of Sobolev embeddings, the initial surfaces 
we consider turn out to be only of C 3//2 class and consequently have unbounded 
curvature, while the initial velocities are only Lipschitz. We also take benefit from 
our low regularity result and an elementary (though seemingly yet unnoticed) 
observation to solve a question raised by Boussinesq on the water waves problem 
in a canal. We reduce the system using a paradifferential approach. 
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1. Introduction 

We are interested in this work in the study of the Cauchy problem for the water 
waves system in arbitrary dimension, without surface tension. 

An important question in the theory is the possible emergence of singularities 
(see [16, 17, 20, 25, 58]) and as emphasized by Craig and Wayne [30], it is impor- 
tant to decide whether some physical or geometric quantities control the equation. In 
terms of the velocity field, a natural criterium (in view of Cauchy-Lipschitz theorem) 
is given by the Lipschitz regularity threshold. Indeed, this is necessary for the "fluid 
particles" motion (i.e. the integral curves of the velocity field) to be well-defined. 
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In terms of the free boundary, there is no such natural criterium. In fact, the sys- 
tematic use of the Lagrangian formulation in most previous works [9, 54, 55], and 
the intensive use of Riemannian geometry tools (parallel transport, vector fields,...) 
by Shatah-Zeng [49, 50, 51], Christodoulou-Lindblad [21] or Lindblad [41] seem 
to at least require bounded curvature assumptions (see also [24] where a logarithmic 
divergence is allowed). In this direction, the beautiful work by Christodoulou- 
Lindblad [21], gives a priori bounds as long as the second fundamental form of the 
free surface is bounded, and the first-order derivatives of the velocity are bounded. 
This could lead to the natural conjecture that the regularity threshold for the water 
waves system is indeed given by Christodoulou-Lindblad's result and that the do- 
main has to be assumed to be essentially C 2 . Our first contribution in this work is 
that this is not the case and that the relevant threshold is actually only the Lipschitz 
regularity of the velocity field. Indeed (see Theorem 1.2), our local existence result 
involves assumptions which, in view of Sobolev embeddings, require only (in terms 
of Holder regularity) the initial free domain to be C 3 / 2 . 

Our second contribution is an illustration of the relevance of the analysis of low 
regularity solutions in a domain with a rough boundary. We give an application of 
our analysis to the local Cauchy theory of three-dimensional gravity water waves 
in a canal. This question was historically at the heart of the analysis of the water 
waves system and goes back to the work by Boussinesq at the beginning of the 20 th 
century (see [15]). Our result seems to be the first one in this setting. 

Our analysis require the introduction of new techniques and new tools. In [1, 2] we 
started a para-differential study of the water waves system in the presence of surface 
tension and were able to prove that the equations can be reduced to a simple form 



where Ty is a para-product and T 7 is a para-differential operator of order 3/2. Here 
the main step in the proof is to perform the same task without surface tension, 
with T-y of order 1/2. It has to be noticed however that performing our reduction is 
considerably more difficult here than in our previous papers ([1, 2]). Indeed, in the 
case with non vanishing surface tension, the natural regularity threshold forces the 
velocity field to be Lipschitz while the domain is actually much smoother (C 5 / 2 ). In 
the present work, the velocity field is also Lipschitz, but the domain is merely C 3//2 . 
To overcome these difficulties, we had to give a micro-local description (and contrac- 
tion estimates) of the Dirichlet-Neumann operator which is non trivial in the whole 
range of C s domains, s > 1 (see the work by Dahlberg-Kenig [31] and Craig-Schanz- 
Sulem [28] for results on the Dirichlet-Neumann operator in Lipschitz domains). We 
think that this analysis is of independent interest. 

Finally, let us mention that, as we proceed by energy estimates, our results are proved 
in L 2 -based Sobolev spaces and our initial data (770, Vp) which describe respectively 
the initial domain as the graph of the function 770 and the trace of the initial velocity 
on the free surface, are assumed to be in H s+ 2 (R d ) x H s (R d ), s > 1 + |. The gravity 
water waves system enjoys a scaling invariance for which the critical threshold is 
s c = 5 + f ■ In other terms our well-posedness result is 1/2 above the scaling critical 
index. We intend to lower the threshold 1 + | by using dispersive estimates in a 
forthcoming paper [3] (see also [2]), following a strategy successfully developed in the 
framework of quasi-linear wave equations by Bahouri-Chemin and Tataru [10, 53] 
(notice that in [3], the Lipschitz threshold is still relevant as even though the initial 
velocity field is only C 1 ' 6 , the solution itself is still L 2 ((— T, T); C 1+e )). 




dtu + Ty ■ Vu + iTryU = f, 
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1.1. Assumptions on the domain. Hereafter, d > 1, t denotes the time 
variable and x G R d and y G R denote the horizontal and vertical spatial variables. 
We work in a time-dependent fluid domain located underneath a free surface £ 
and moving in a fixed container denoted by O. This fluid domain 

Q = { (t, x, y) G [0, T] x K d x R : (<c, y) G O(t) }, 

is such that, for each time t, one has 

Sl(t) = {(x,y) € O : y<T/(t,x)}, 

where 77 is an unknown function and is a given open domain which contains a 
fixed strip around the free surface 

E = {(t,x,y) G [0,T] x R d x R : y = V (t,x)}. 

This implies that there exists h > such that, for all t G [0, T], 

(1.2) n h (t) := {(x,y) G R d x R : rj(t,x) -h<y < rj(t,x)} C fi(i). 

We also assume that the domain (and hence the domain O(i)) is connected. 

Remark 1.1. (i) Two classical examples are given by O = R d x R (infinite depth 
case) or O = YL d x [— l,+oo) (flat bottom). Notice that, in the following, no 
regularity assumption is made on the bottom T := dO. 
(ii) Notice that T does not depend on time. However, our method applies in the 
case where the bottom is time dependent (with the additional assumption in 
this case that the bottom is Lipschitz). 

1.2. The equations. Below we use the following notations 

V = (d Xi )i<i< d , V x , y = (V,d y ), & = J2 1 < i < d d %i> A Xiy = A + d 2 y . 

We consider an incompressible inviscid liquid, having unit density. The equations by 
which the motion is to be determined are well known. Firstly, the eulerian velocity 
field v : Q — > R d+1 solves the incompressible Euler equation 

(1.3) d t v + v ■ V X)V v + V x , y P = -gey, div Xt yV = in 0, 

where — ge y is the acceleration of gravity (g > 0) and where the pressure term P 
can be recovered from the velocity by solving an elliptic equation. The problem is 
then given by three boundary conditions. They are 

' v ■ n = on r, 

(1.4) < dtv = V 1 + |V??| 2 t; • v on E, 

P = on E, 

v ■ 

where n and v are the exterior unit normals to the bottom T and the free surface £(i). 
The first condition in (1.4) expresses the fact that the particles in contact with the 
rigid bottom remain in contact with it. Notice that to fully make sense, this condition 
requires some smoothness on T, but in general, it has a weak variational meaning 
(see Section 3). The second condition in (1.4) states that the free surface moves 
with the fluid and the last condition is a balance of forces across the free surface. 
Notice that the pressure at the upper surface of the fluid may be indeed supposed 
to be zero, provided we afterwards add the atmospheric pressure to the pressure so 
determined. 

The fluid motion is supposed to be irrotational. The velocity field is therefore given 
by v = 'V Xt y4> for some potential 0: Q — > R satisfying 

^x,y<ft = in S7, d n (p = on T. 
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Using the Bernoulli integral of the dynamical equations to express the pressure, the 
condition P = on the free surface implies that (recalling that V = V x ) 

dtrj = d y (j) — Vr] ■ V4> on S, 



(1.5) { dtij> + -\V Xi y<l>\ 2 +gv = onS, 

d n (j) = on r. 



Many results have been obtained on the Cauchy theory for System (1.5), starting 
from the pioneering works of Nalimov [45] , Shinbrot [52] , Yoshihara [59] , Craig [26] . 
In the framework of Sobolev spaces and without smallness assumptions on the data, 
the well-posedness of the Cauchy problem was first proved by Wu for the case without 
surface tension (see [54, 55]) and by Beyer-Giinther in [13] in the case with surface 
tension. Several extensions of their results have been obtained by different methods 
(see [23, 32, 33, 34, 37, 42, 56, 57, 61] for recent results and the surveys [12, 
30, 38] for more references). Here we shall use the Eulerian formulation. Following 
Zakharov [60] and Craig-Sulem [29], we reduce the analysis to a system on the free 
surface S(t) = {y = r](t,x)}. If is defined by 

ip{t,x) = 4>(t,x,rj(t,x)), 

then (f> is the unique variational solution of 

A Xt y(p = in 0, 0| s = ip, d n (p = on T. 

Define the Dirichlet-Neumann operator by 



(G( V )il>){t,x) = V^+WWdn4>\y =v{t , x) 

= {d y (p)(t, x, rj(t, x)) - Vrj(t, x) ■ (V(f))(t, x, r)(t, x)). 

For the case with a rough bottom, we recall the precise construction later on 
(see §3.1). Now (r],ijj) solves (see [29] or [38, chapter 1] for instance) 



(1.6) 



f Otv - G(r?)V = 0, 

fl/o. ^ 1 \xrn 2 1 ( V? ?-W + W 2 n 
W + gr) + -m fT ^ 1 = 0. 



1.3. The Taylor condition. Introduce the so-called Taylor coefficient 

(1.7) a{t,x) = -(dyP){t,x,i ] {t,x)). 

The stability of the waves is dictated by the Taylor sign condition, which is the 
assumption that there exists a positive constant c such that 

(1.8) a{t,x)>c>0. 

This assumption is now classical and we refer to [12, 21, 22, 39, 54, 55] for 
various comments. Here we only recall some basic facts. First of all, as proved by 
Wu ([54, 55]), this assumption is automatically satisfied in the infinite depth case 
(that is when T = 0) or for flat bottoms (when T = {y = —k}). Notice that the 
proof remains valid for any C 1,Q -domain, < a < 1 (by using the fact that the Hopf 
Lemma is true for such domains, see [47] and the references therein). There are two 
other cases where this assumption is known to be satisfied. For instance under a 
smallness assumption. Indeed, if dt<fi = 0(e 2 ) and V Xjy (j) = O(e) then directly from 
the definition of the pressure we have P — gy = 0{e 2 ). Secondly, it was proved by 
Lannes ([39]) that the Taylor's assumption is satisfied under a smallness assumption 
on the curvature of the bottom (provided that the bottom is at least C 2 ). However, 
for general bottom we will assume that (1.8) is satisfied at time t = 0. 



[1.9) B= \ :" U \ V = V^-BV V . 



1.4. Main result. We work below with the vertical and horizontal traces of 
the velocity on the free boundary, namely 

B := (d y 4>)\ y=v , V := (V x 4>)\y =ri . 
These can be defined only in terms of 77 and tp by means of the formulas 

Vr] ■ Vtp + G(n)ip 
I + IVt/P 

Also, recall that the Taylor coefficient a defined in (1.7) can be defined in terms 
of 77, V, B, i/j only (see Section 1.6 below). 

Theorem 1.2. Let d > 1, s > 1 + d/2 and consider (770, V'o) such that 

(1) 770 G iT s+ i(R d ), r/> G # s+ i(R d ), Vb € -ff s (R d ), S Gi? s (R d ), 
i/iere exists h > suc/i i/iai condition (1.2) /10/ds initially for t = 0, 
i/iere exists a positive constant c suc/i i/iat, /or aZZ x G R d , ao( :c ) > c - 

T/ien i/iere exists T > suc/i i/iai t/ie Cauchy problem for (1.6) mi/i initial data 
(770,^0) /ias a Ttnigue solution (rj,ip) G C°([0, T]; iJ s+ 2 (R d ) x F s+ ^ (R d )) , suc/i tftai 

fij we /ia?je (V,B) G C°([0,T\; H s (R d ) x F s (R d )), 

^ the condition (1.2) ZioWs /or < t < T, mf/i /i replaced by h/2, 

(3) for allO<t<T and for all x G R d , a(t, x) > c/2. 

Remark 1.3. The main novelty is that, in view of Sobolev embeddings, the ini- 
tial surfaces we consider turn out to be only of C 3 / 2 class and consequently have 
unbounded curvature. 

Remark 1.4. Assumption 1 in the above theorem is automatically satisfied if 
770 G H s+ ^(R d ), -00 G H^(R d ), V G H s (R d ), B G L 2 (R d ). 

The only point where the estimates depend on ip (and not only on 77, V, B) come 
from the fact that we consider a general domain without assumption on the bottom. 
Otherwise, we shall prove a priori estimates for the fluid velocity and not for the 
fluid potential (notice that the fluid potential is defined up to a constant). 

1.5. Water waves in a canal. We give here an illustration of the relevance of 
the analysis of low regularity solutions in a domain with a rough boundary. We show 
that the above result (or rather a slight adaptation of this result) allows to prove 
the existence of three-dimensional gravity water waves in a canal having vertical 
walls near the free surface. The propagation of waves whose crests are orthogonal 
to the walls is one of the main motivation for the analysis of 2D waves. It was 
historically at the heart of the analysis of water waves. The study of the propagation 
of three-dimensional water waves for the linearized equations goes back to Boussinesq 
(see [15]). However, there are no existence results for the nonlinear equations in the 
general case where the waves can be reflected on the walls of the canals (except 
the analysis of 3D-periodic travelling waves which correspond to the reflexion of a 
2D-wave off a vertical wall, see Reeder-Shinbrot [46], Craig and Nicholls [27] and 
Iooss-Plotnikov [36]). More precisely, we consider a fluid domain which at time t is 
of the form 

Sl(t) = {(x 1 ,x 2 ,y) G (0, 1) x R x R : b(x) < y < T)(t, a?)} , 

for some given function b. We do not make any regularity assumption on b; again, 
our only assumptions are that: (i) there exists a positive constant h such that 
77(4, x) > b(x) + h and (ii) the Taylor sign condition holds initially. 
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Figure 1. A non rectangular canal with vertical walls near the free surface. 



In this context, we are able to prove that the Cauchy problem for the water waves 
system in this domain is well posed (see Section 7 for a precise statement). The 
analysis is based on an elementary (though seemingly yet unnoticed) observation: 
for any solution of the free boundary Euler equation having non vanishing Tay- 
lor coefficient, the free surface necessarily makes a right-angle with the rigid walls 
(see Figure 2 and Proposition 7.1). This suggests, following Boussinesq (see [15, 
page 37]) to perform a symmetrization process, as illustrated on Figure 3. 




Figure 2. Two-dimensional section of the fluid domain. 




Figure 3. Two-dimensional section of the extended fluid domain. 



Notice on Figure 3 that, to handle non rectangular canals, we have to work with 
rough bottoms (even if initially the bottom is smooth, after symmetry/periodization 
this is no more the case). Of course, the symmetry should yield in general a Lipschitz 
singularity for the symmetrized free surface. However, here the possible singularities 
are weaker since the above physical observation about the right angles at the interface 
implies that d xi r](t,0, X2) = and d xi rj(t, 1, X2) = 0. Therefore, we are in position 
to apply our low regularity Cauchy theory. Namely, we apply Theorem 1.2 for 
some s < 3. Notice however that we need a small adaptation of our main result as 
we consider domains in T x R 2 instead of R 3 . We refer to Section 7 for details on 
how to adapt the result to the periodic framework. 
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1.6. The pressure. The purpose of this paragraph is to clarify, for low regu- 
larity solutions of the water waves system in rough domains, the definition of the 
pressure which is required if one wants to come back from solutions to the Zakharov 
system to solutions to the free boundary Euler equation. This definition will also 
provide the basic a priori estimates which will be later the starting point when estab- 
lishing higher order elliptic regularity estimates required when studying the Taylor 
coefficient a = —d y P |s- On a physics point of view, the pressure is the Lagrange 
multiplier which is required by the incompressibility of the fluid (preservation of the 
null divergence condition). As a consequence, taking the divergence in (1.3), it is 
natural to define the pressure as a solution of 

(1.10) ^x,yP = ~ <^x,y(v ■ V x ,yV), P \y= V = 0. 

Notice however that the solution of such problem may not be unique as can be seen 
in the simple case when Q = (—00, 0) x R d . Indeed, if P is a solution, then P + cy 
is another. Notice also that if P satisfies (1.10), then 



& x , v (p + gy + -\vf) =0. 



Definition 1.5. Let (77,^) G (W 1 * 00 n H^ 2 (R d )) x H^ 2 (R d ). Assume that the 
variational solution (as defined in §3.1 ^ of the equation 

(1.11) &x,y<t> = 0, (j) \ y = v = ip, 

satisfies 

\V XtV <f>\ 2 (x,ri(x)) etfi(R d ). 
Let R be the variational solution of 



A.,..,,./? /;-/ il R \ y=v = - [gv+ ^\^x,y(, 



We define the pressure P in the domain £1 by 



P(x, y) := R(x, y) - gy - -\V x , y <t>{x, y)\ 2 . 

Remark 1.6. The main advantage of defining the pressure as the solution of a 
variational problem is that it will satisfy automatically an a priori estimate (the 
estimate given by the variational theory). 



It remains to link the solutions to the Zakharov/Craig-Sulem system (1.6) to solu- 
tions of the free boundary Euler system (1.3) with boundary conditions (1.4). To 
do so, we proved in [4] that if (r], ifj) is a solution of System (1.6), if we consider the 
variational solution to (1.11), then the velocity field v = V x ^ y (j) satisfies (1.3), which 
is of course equivalent to 

(1.12) P = -d t <fi -gy-^V^l 2 . 

Theorem 1.7 (from [4]). Assume that (rj,ip) G C([0,T];# s+ 5(R d ) x # s+ 5(R d )) ; 
with s > 1 + d/2, is a solution of the Zakharov/Craig-Sulem system (1.6). Then the 
assumptions required to define the pressure are satisfied, and (1.12) is satisfied, and 
the distribution dt4> is well defined for fixed t and belongs to the space -ff 1,0 (f2(i)) 
(see Definition 3.3). 
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1.7. Plan of the paper. At first glance, Theorem 1.2 looks very similar to our 
previous result in presence of surface tension [1, Theorem 1.1]. Indeed, the regularity 
threshold exhibited by the velocity field (namely V, B G H s (H d ), s > 1 + d/2) is the 
same in both results and (as explained above) appears to be the natural one. How- 
ever, an important difference between both cases is that the algebraic nature of (1.6) 
(and its counter-part in presence of surface tension) requires that the free domain 
is 3/2 smoother than the velocity field in presence of surface tension and only 1/2 
smoother without surface tension. This algebraic rigidity of the system implies that 
in order to lower the regularity threshold to the natural one (Lipschitz velocities), 
we are forced to work with C 3 / 2 domains (compared to the much smoother C 5 / 2 
regularity in [1]). This in turn poses new challenging questions in the study of the 
Dirichlet-Neumann operator. Indeed, at this level of regularity the regularity of the 
remainder term in the paradifferential description of the Dirichlet-Neumann opera- 
tor G(rj)ip is not given by the regularity of the function tp itself, but rather by the 
regularity of the domain. This is this phenomenon which forces us to work with the 
new unknowns V, B rather than with ip. 

In Section 2, we wrote a review of paradifferential calculus and proved various tech- 
nical results useful in the article. In Section 3 we study the Dirichlet-Neumann 
operator and prove the main estimate (see Theorem 3.17) that we believe is of in- 
dependent interest. In Section 4, we symmetrize the system and prove a priori 
estimates. In Section 5 we prove the contraction estimates required to show unique- 
ness and stability of solutions. In particular we prove a contraction estimate for the 
difference of two Dirichlet-Neumann operators, involving only the Cz norm of the 
difference of the functions defining the domains (see Theorem 5.3). In Section 6 we 
prove the existence of solutions by a regularization process. Finally, in Section 7 we 
study the Cauchy problem in the canal. 



Let us review notations and results about Bony's paradifferential calculus. We refer 
to [14, 35, 43, 44] for the general theory. Here we follow the presentation by 
Metivier in [43]. 

2.1. Paradifferential operators. For k G N, we denote by W k ' 00 (R. d ) the 
usual Sobolev spaces. For p = k + a, k G N,a G (0, 1) denote by Pv^>°°(R d ) the 
space of functions whose derivatives up to order k are bounded and uniformly Holder 
continuous with exponent a. 

Definition 2.1. Given p G [0, 1] and m G R, T™(R d ) denotes the space of locally 
bounded functions a(x,£) on H d x (R d \0), which are C°° with respect to £ for £ ^ 
and such that, for all a G N rf and all £ ^ 0, the function x i— > d^a(x, £) belongs 
to W p,00 (R d ) and there exists a constant C a such that, 



2. Paradifferential calculus 




Then F^R^) denotes the subspace of r™(R ci ) which consists of symbols a(x,£) 
which are homogeneous of degree m with respect to £. 



Given a symbol a, we define the paradifferential operator T a by 



(2.1) 
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where a(9,£) = fe ix a(x, £) dx is the Fourier transform of a with respect to the 
first variable; x anci ip are two fixed C°° functions such that: 

(2.2) VO?) = for \t]\ < 1, ^(77) = 1 for \rj\ > 2, 

and x(#) ^) satisfies, for < £1 < £2 small enough, 



x(M) 

and such that 



1 if |0|<eiM, x (9, v ) = if \9\>e 2 \v\, 



V(M) : |^x(M)| <C a)j9 (l + H)-l a H^I. 

Remark 2.2 (Choice of the cut-off function). It is known that one can change 
the cut-off function % above, up to adding some smoothing remainders (see [43, 
Proposition 5.1.8]). However, this requires to work with smooth enough symbols. 
Below, we shall need to work with rough symbols, so that the choice of the cut-off 
function can, sometimes matter. It will therefore be convenient to fix the cut-off 
function so that all the computations below can be done with the same choice. 
Since we need to work with paraproducts, we chose a cut-off function \ such that 
when a = a(x), T a is given by the usual expression in terms of the Littlewood-Paley 
operators. Namely, we introduce k S C^°(R d ) such that 



k{9) = 1 for \9\ < 1.1, k(9) = for \9\ > 1.9. 



Then we define 



+00 

£ 

fc=0 



K k - 3 (9)<p k (jl) 



where 

K k (9) 



K{2~ k 9) for k £ Z, 



^0 = K 0> an d fk = K 'k~ K k-1 f° r fe > 1. 



2.2. Symbolic calculus. We shall use quantitative results from [43] about 
operator norms estimates in symbolic calculus. To do so, introduce the following 
semi-norms. 

Definition 2.3. For m e R ; p E [0, 1] and a e r™(R d ), we set 
(2.3) M™(a)= sup sup (1 + K|)H-^ a (-,0 

H <f +l+p |€[>l/2 ^-(R") 

Definition 2.4 (Zygmund spaces). Consider a dyadic decomposition of the identity: 
I = A_i + X^o^<?- ^ s an y rea ^ number, we define the Zygmund class C^(R d ) 
as the space of tempered distributions u such that 

\\u\\ C s := sup 2 qs ||A g u|| Loo < +00. 
<? 

Remark 2.5. It is known that C^(R d ) is the usual Holder space W s '°°('R d ) if s > 
is not an integer. 

Definition 2.6. Let m 6 R. An operator T is said to be of order m if, for all fx € R, 
it is bounded from to H ll ~' m and from to C*~ m . 

The main features of symbolic calculus for paradifferential operators are given by 
the following theorem. 
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Theorem 2.7. Let m G R and p G [0, 1]. 

(i) If a G T™(R d ), i/ten T a is o/ order m. Moreover, for all p G R f/iere exists a 
constant K such that 

(2.4) ||T || HM _, HM _ m < ^M m (a), ||T a || c ^ cr m < M m (a). 

(ti) If a £ T m (R d ),b G r™'(R d ) t/ien T a T fe - T ab is of order m + m'-p. Moreover, 
for all p G R £/iere exists a constant K such that 

\\T a T b - T ab \\ H ^ H ^ m , +p < KM™(a)M^'(b) + KM™(a)M™' (b), 
\\T a T b - T ab \\ c ^ crm . ml+P < KM™(a)Mf(b) + KM^(a)Mf(b). 

(Hi) Let a G T™(R a! ). Denote by (T a )* the adjoint operator of T a and by a the 
complex conjugate of a. Then (T a )* — Ta is of order m — p. Moreover, for all p there 
exists a constant K such that 
(2.6) 

II (T a )* - Ta\\ H ^ H ^ m+P < KM™ (a), \\(T a )* - Ta\\ c ^ c ,-m +P < KM™ (a). 

We shall need in this article to consider paradifferential operators with negative 
regularity. As a consequence, we need to extend our previous definition. 

Definition 2.8. For m G R and p G (— oo,0), T™(R d ) denotes the space of distri- 
butions a(x,£) on H d x (R d \ 0), which are C°° with respect to £ and such that, for 
all a G N d and all £ / 0, the function x i-> d^a(x,^) belongs to C*(R. d ) and there 
exists a constant C a such that, 

(2-7) V l^' \\d?a(;Q\\ cS <C a (l + \t\r-W. 

Then T m (R d ) denotes the subspace of r™(R ci ) which consists of symbols a(x,£) 
which are homogeneous of degree m with respect to £. For a G T m , we define 



(2.8) M m (a) = sup sup (1 + \^\) H ' m d?a(-, £) 



\a\<f+ P +l |ei>i/2 



Cf(R d ) 



2.3. Paraproducts and product rules. If a = a(x) is a function of x only, the 
paradifferential operator T a is called a paraproduct. A key feature of paraproducts 
is that one can replace nonlinear expressions by paradifferential expressions, to the 
price of error terms which are smoother than the main terms. Also, one can define 
paraproducts T a for rough functions a which do not belong to L°°(R d ) but merely 
C~ m (R d ) with m > 0. 

Definition 2.9. Given two functions a,b defined on H d we define the remainder 

R(a, u) = au — T a u — T u a. 

We record here various estimates about paraproducts (see chapter 2 in [11] or [19]). 

Theorem 2.10. i) Let a, /3 G R. If a + f3 > then 

( 2 -9) W R ( a > u )W H *+f)-$ {Rd) < K \\ a \\H-(R" ) \\u\\ H ^ Rd) , 

(2.10) || R(a, u) || c «+/3 (R d) < ^ l|a||c ? (Rd) IMI c f(Rd) > 

(2.11) \\R(d, u)|| H c+/3( R d) < K ||a|| C a( R d) |M|tf/3( R d) • 
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ii) Let m > and s G R. T/ien 

(2.12) ||T|jU|| jy s _ m ^ ^ Holler™ ll^llif 8 ' 

(2.13) ||r a n|| c s-m < K \\a\\ c -m \\u\\ Cs . 

(2.14) l|7>|lcs ^ K \\ a \\h°° \\ u \\c* ■ 

Hi) Let so, s\,S2 be such that sq < s 2 and so < s± + s 2 — \, then 
(2-15) \\T a u\\ HH) <K\\ 

By combining the two previous points with the embedding i^ M (R d ) C C* d ^ 2 (R d ) 
(for any /i 6 R) we immediately obtain the following results that we shall need in 
the sequel. 

Proposition 2.11. Let r,fj, GR be such that r + / u>0. d/7 G R satisfies 

d 

7 < r and 7 < r + [i — — , 

then there exists a constant K such that, for all a G H r (R d ) and all u G H fl (R d ), 
we have 

\\au - T a u\\ Hl < K || 

Corollary 2.12. i) If uj G H s i(R d ) (j = 1,2) with si + s 2 > then mu 2 G 
iF°(R d ) and 

(2.16) ll^i^H^o < ^ H^ill/pi ||n 2 ||^ S2 , 
if 

so < Sj; i = 1) 2, and so < Si + s 2 — d/2. 
ii) (Tame estimate in Sobolev spaces) If s > t/ien 

(2.17) ||niu 2 || Hs < Jf(||ui|| Hs ||n 2 || Loo + ||ni|| ioo ||n 2 || Hs ). 
Hi) (Tame estimate in Zygmund spaces) If s > then 

(2.18) ||uiu 2 || C s < if (||ui|| C s ||n 2 || Loo + ||ni|| LO o ||n 2 || C s). 
iv) Let fj,,m G R be such that fi, m > and m^N. Then 

(2.19) \\uiu 2 \\ Hf , < K{\\ "Uill^oo ||n 2 ||^ M + ||n 2 || ( -,-m ||ni||^ M + m ). 

v) Let (3 > a > . Then 

(2.20) ||nin 2 || c ,- ct < K \\ui\\ c p H^H^-a . 

vi) Let s > d/2 and consider F G C 00 (C iV ) such that F(0) = 0. Then there exists 
a non- decreasing function T : R + — > R + such that 

(2-21) im^)llH.<J r (ll^llLoo)||^|lH., 

for any U G fT 8 (R d ) Jv . 

vii) Let s > and consider F G C°°(C N ) suc/i i/iai F(0) = 0. Then there exists a 
non- decreasing function T : R + — > R + such that 

(2-22) Unfile. <HW U Wl~)\\U\\c V 

for any U G C°{R d ) N . 
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Proof. The first three estimates are well-known, see Hormander [35] or Chemin 
[19]. To prove iv) and v) we write 

Ul u 2 = T Ul u 2 + T U2 ui + R(ui,u 2 ). 

Then (2.19) follows from 

\\T Ul U2\\ H n ^ \\ui\\ L oo \\u2\\ HM (see (2.4)), 

\\ t u2 u i\\h>+ ~ ll n 2|| c -m (see (2.12)), 

\\R(ui,u 2 )\\ Hl t < \\u2\\ c -m ||tti|| H „+ m (see (2.11)), 

while similarly (2.20) follows from 

ll^u2 w illc^ Q ~ ll^llc--" Hollies — ll^llc?-" Ihillcrf > 
\\R(ui,u 2 )\\ c - a < \\R(ui,U2)\\ G p- a < ||«2|| c -« \\ui\\ c f) ■ 

(With regards to the last inequality, to apply (2.10) we do need j3 > a > 0.) Finally, 
vi) and vii) are due to Meyer [44, Theoreme 2.5 and remarque], in the line of the 
work by Bony [14]. □ 

Finally, let us finish this section with a generalization of (2.12) 

PROPOSITION 2.13. Let p < 0, m e R and a E T™. Then the operator T a is of 
order m — p: 

\\T a \\ H s_+ H s-( m -p) < CM™(a), 

(2 ' 23) ll^Ucj^cr^ - CA W- 

Proof. Let us prove the first estimate. The proof of the second is similar. 
Notice that if m = and a(x,£) = a(x), then (2.23) is simply (2.12). Furthermore, 
if a(x, £) = b(x)p(£), then 

T a = T b ( X p)(\D x \), 

where x 1S a cutoff function vanishing near and equal to 1 for |£| > 1. As a 
consequence, in this particular case, we get 

\\T a \\ H s^ H s-(m-p) < C\\b\\ C p\\p \ S d-l 1 1 2,oo. 

In the general case, we can expand, for fixed x, a(x,£) in terms of spherical harmon- 
ics. Let (fej/)i/gN* be an orthonormal basis of L 2 (S d ~ l ) consisting of eigenfunctions 
of the (self-adjoint) Laplace-Beltrami operator, = A s <i-i on L 2 (S d_1 ), i.e. 

By the Weyl formula, we know that X u ~ cv~3. . Setting 



we can write 
a 



(x,£)= y a u {x)h v {£) where a u {x) = I a(x,uj)h u (uj) dco. 



1/6N 

Since 

v 2fc „ i j. \ I a k 



Xfa v {t,x) = \ Ala(x,u)h v (u>)du), 

we have, for all v > 1, 



-i 



3d 



(2.24) \K(-)\\c» < CXu 2 t1 < i/"i-3AC(p). 
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Moreover, 

(rf-l) ! , 

(2.25) HML« - CXv 2 <C^2~2S. 

and the result follows because 

||T a || fl »_ fl .-(»- P ) < Cj^^^M^ip). 

This completes the proof. □ 

We shall also need the following technical result. 
Proposition 2.14. Set (D x ) = (I - A) 1 / 2 . 

i) Let s > | + t| and o~ G R 6e suc/i £/ia£ cr < s. TTien t/iere exists K > suc/i i/iai 
/or a// V 6 W 1 '°°(R d ) n H s (R d ) and u G H a -^(R d ) one has 

\\[(D x ) a ,V]u\\ L 2 iRd) < K{\\V\\ w i,o p L d) + ||V|| H .( R «i)}||u|| H<r _^ . 

ii) Let s > 1 + | and aGRfc such that a < s. Then there exists K > such that 
for all V G H s (R d ) and u G £T ff - 1 (R d ) one has 

\\[{D x ) a ,V]u\\ L 2 (R d) < if||V|| H! .( R d)||u|| H <7-i( R d). 

iii; lei s > § + jj and V G H s (R d ). Then 

\\[{D X )3 , V r ]«|| i ao( R d) < i^||^||jf S ( R d)||li||ioo( R d). 

Proof. To prove z) we write 

f || [(^,^11^ <^ + B, 

\A= ||[(^r , T y ]n|| L2 , B = || [(D !8 ) ff , V - T v ]u\\ L 2. 
By (2.5) we have 

A < K\\V\\ W l,oo\\u\\ H a-l. 

On the other hand one can write 

B< \\(D x y {V -T v )u\\ L 2 + \\(V -T v )(D x y u\\ L 2 =B 1 +B 2 . 

We use Proposition 2.11 two times. To estimate B\ we take j = a,r = s,fi = a — ^. 
To estimate B 2 we take 7 = 0, r = s, [i = —\ and we obtain, 

b < #||y||tf S ||u|| a i . 

- II \\n II ll H a- 3 

To prove u), to estimate i?i (resp. .B2) we use again Proposition 2.11 with 7 = a, r = 
s, = <r — 1 (resp. 7 = 0,r = s, /i = —1). Finally, let us prove Hi). Using (2.5) 
with m = 5, m' = 0, p = \ + e, we obtain 

HK-Da;)^ ,Ty]n||c| < C||V||.H*|| w llc° < C|I^II-H" S II W IU 00 - 

On the other hand, 

[(D a )3 , V - T y ]n = (D,)* (V - T y )n - (F - 7V)(Lg§ u. 
Let ^ < r < s — I so that 

ll^lk < cil^llfl*- 

According to (2.14) and (2.10), V - T v is bounded from L°° to C£ by K ||V||cj and 
according to (2.13) and (2.10), from C* 2 to C* 2 by if[|V||cj, which implies 

||[(D x )3,y-Tvr]tt|| r _l <if||V||H.||u||Loo. 
13 



This completes the proof. 



□ 



We shall need elementary estimates on the solutions of transport equations that we 
recall now. 



(2.26) 



Proposition 2.15. Let I = [0,T] and consider the Cauchy problem 

d t u + V •V« = /, tel, 
u\ t =o = u . 

We have the following estimates 

(2.27) IK*)IIl«(r*) < IML«(R<i)+ / ll/(OIL°° (R d)d<7. 

J 

There exists a non decreasing function T : R + — > R + such that 

(2.28) ||u(t)|| i2(Rd) < J c "(||y|| L i (7;H/ i, O c (R d )) )(||n || L 2 (R d ) + I \\f(t',-)\\ L 2 {Rd) dt'). 

Jo 

If s > 1 + 5 and a < s there exists a non decreasing function T : R + — > R + such 
that 

(2.29) \\u(t)\\ HrT(Rd) < -7 r (||F|| L i ( / ; ^(Rd)))(lko||H-(Rd) + / ll/(*'>-)l|jf-(Rd)^ / )- 

JO 

Proof. The estimate (2.27) follows from the exact formula for the solution, in 
terms of the solution to the differential equation 

X = -V{t,X(t)), 

while (2.28) is obtained by computing ^\\u(t, Oll^^d) and (2.29) follows from (2.28) 
applied to (D x ) a u, Proposition 2.14 and Gronwall inequality. □ 

2.4. Commutation with a vector field. We prove in this paragraph a com- 
mutator estimate between a paradifferential operator T p and the convective deriv- 
ative dt + V ■ V. Inspired by Chemin [18] and Alinhac [7], we prove an estimate 
which depends on estimates on dtp + V ■ Vp and not on Vt t xP- 

When a and u are symbols and functions depending on t £ I, we still denote 
by T a u the spatial paradifferential operator (or paraproduct) such that for all t £ I, 
(T a u)(t) = T a ^u{t). Given a symbol a = a(t;x,^) depending on time, we use the 
notation 



Mtf(a):= sup sup sup (1 + |e|) |a| ~ m ^o(*; •, £) 

te[0,T] H <w +1+p |£|>i/ 2 



L°°(R d ) 



Given a scalar symbol p = p(t, x, £) of order m, it follows directly from the symbolic 
calculus rules for paradifferential operators (see (2.4) and (2.5)) that, 

\\[T p ,d t + T v -V]u\\ m < KiM^id^+M^iVp^WVW^WuW^. 

A technical key point in our analysis is that one can replace this estimate by a 
tame estimate which does not involve the first order derivatives of p, but instead 
dtP + V ■ Vp. 
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Lemma 2.16. Let V G C°([0,T];Cl +£ (R d )) for some e > and consider a sym- 
bol p = p(t,x,£) which is homogeneous in £ of order m. Then there exists K > 
( independent of p,V ) such that for any t G [0,T] and any u G C°([0, T]; H m (R d )). 

(2.30) ||[r p ,5 t + ryV]«(t)|| i2(Rd) 

Proof. Set I = [0,T] and denote by 1Z the set of continuous operators R(t) 
from fF™(R d ) to L 2 (R d ) with norm satisfying 

\\m\\c iHm( n WRd)) < K{M^(p) \\V(t)\\ ci+s +MV(dtp + V.Vp)\\V(t)\\ LOB }. 
We begin by noticing that it is sufficient to prove that 

(2.31) (d t + V ■V)T p = T p (d t + T v -V) +R, ReTZ. 
Indeed, by Theorem 5.2.9 in [43], we have (for fixed t) 

|| (V - T v ) ■ VT p u\\ L2 < HVll^i.oo ||r p u|| L2 , 
and hence, by using the operator norm estimate (2.4), we find that 
\\(V - T v ) ■ VT p u\\ L2 < \\V\\ wl>ao M^{p) \\u\\ Hm , 
which implies that (V — Ty) • VT p G 1Z. 

We split the proof of (2.31) into three steps. By decomposing p into a sum of 
spherical harmonic, we shall reduce the analysis to establishing (2.31) for the special 
case when T p is a paraproduct. In the first step we prove (2.31) for m = and p = 
p(t, x). In the second step we prove (2.31) for p = a(t, x)h(£) where h is homogeneous 
in £ of order m. Then we consider the general case. 

Step 1: Paraproduct, m = 0, p = p(t,x). Notice that in this case 

M° (p) = \\ P \\ L ~. 

We have 

j d t T p u = T dtp u + T p d t u, 
\ V ■ VT p u = V • T Vp u + VT p -Vu=: A + B. 



(2.32) 



Decompose V = Sj^ 3 (V) + S j - 3 (V), with 

k<j-2 k>j-3 

to obtain 

' A = Ai + A 2 , 
A, :=^5 i _ 3 (K)5 i _ 3 (Vp)A i «, 



(2.33) 



A 2 := J^'-^Sj-aCVjOAjt*. 



3 

Let us consider the term A 2 . Since 



\&- 3 (y)\\ LOB < £ IIA*v|| L oc< Yl ^ k{1+£) \\v\\ cl+ , <2-^ \\v\\ cl+s 



k>j-3 k>j-3 
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and ||S , i _ 3 (Vp)|| ioo < 2 J \\p\\ Loo , we obtain 



< 



(2.34) 



II 2 



<M° (p)\\V\\ cl+ , \\u\\ L2 . 
We now estimate A\ = An + A12, with 



(2.35) 



A12 :=X;{[^'" 3 (^).^-3]Vp}A jU . 



Write Sj- 3 (V) = V - Si- 3 (V), to obtain 

An = E ^-3(^ • A ^ " E ^-3{^'- 3 (^) • Vp} A. 



Ty-vpu + I + II 



where 



1 = ~ E( V • Sh{S H (V)p})^. ^ = E Sj- 3 { si ~ 3 (y ■ V)p}A jU . 



Then 



Moreover, 



\\i\\ L * <Y, 2j \\ sj - 3 ( v M\ L ~ H a ^IIl 2 

i 

<E2 j 2-^ 1+£ )||y|| ci+e |b|| ioo ||u|| L2 

3 

< \\V\\ ci+ e IbHioc IHI^- 

I^Hl 2 <EH^ 3 ( w )IU IIpIIl- H a ^IIl 2 



< 



ll^|| C l + e 



L 2 • 



Therefore 
(2.36) 



^-11 = Ty.x] p u + R £ TZ. 



In order to estimate A12, note that one can replace Vp by Sj_3(Vp) where Sj- 3 = 
^(2-0- 3 )D) for some function ^ G C£°(R d ) such that V>(£) = 1 for |£| < 2. Next, 
observe that 

A12 = E( [^ 3 (^),Si-3]V5i-3(p)}A 3 « = 

where is spectrally supported in an annulus \c\l 3 < |£| < C22- 7 }, Cj > 0. These 
annuli have only finite overlap, thus by Plancherel we have 

Il^i2||| 2 < E \\{[S'-\V),S 3 - 3 ]VS^ 3 (p)}A 3 u 12 



L 2 



<E 2^11^11^2^11^11^ \\A jU \\l> 



< 



\\n c i+< 



u 



L 2 > 
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where we used the fact that the commutator [S^ 3 (V), Sj-3] is of order —1 (uni- 
formly in j), since V G C°([0, T]; VF 1 ' 00 ). It follows that A i2 = Ru with R £ 7Z. 
Consequently, we deduce from (2.35) and (2.36) that 

Ai = Ty-vpU + Ru, ReU. 

It thus follows from (2.33) and (2.34) that 

A = Ty.xjpU + Ru, R £ 1Z. 

It remains to estimate the term B introduced in (2.32). Again, we split this term 
as follows: 

B = V ■ (T p Vu) = V Sj- 3 (p)VA jU 

j 

= E ^- 3 (^)^-3(p)A i Vn + ^" 3 (^)S i - 3 (Vp)A i Vn =: B 1 + B 2 . 

j j 

We have 

my <EII^" 3 ( y )IL- \\ s j-3( P )\\ L ~ \\^vu\\ L2 

3 

<J2^ j{1+£) \\v\\ c ^ * \\p\\l~ \Hv ■ 
j 

and hence B 2 = Ru with R G 1Z. To deal with the term B\, let us introduce 

(2.37) C := T P T V -Vn = ^ S,-- 3 (p)A,- E S k-s(V) ■ VA fc u. 

j k 

Since the spectrum of S k - 3 (V) • VA fc u is contained in{(3/8)2 fc < |£| < (2 + l/8)2 fc }, 
the term Aj(Sks(V) ■ VA^n) vanishes unless \k — j\ < 3. On the other hand, 

for \k — j| < 3, Sk-z(V) — Sjs(V) = ± X^=-jv A -£+j^' anc ^ hence we can write C 
under the form 

C = C 1 + C 2 = C 1 + J2Sj-3(p)A j {S j - 3 (V)- £ VA W 

j |fc-il<3 

where Ci is given by 

3 i-2 

d = J^Sj-sW^j E E (A£ +J (y)VA i+J (n) - A, +j _ i (y)VA J _ i (n)}, 

i i=l £=-1 

so that 

l|cy La < E IIpIL- 2 ^' (1+£) II^Hc^ 2J 1Mb . 

3 

which implies that C\ = Ru with R G TZ. To estimate C 2 , as before we write 
C 2 = C 2 i + C 22 where 

C 21 :=E^-3(P)[A„5,_ 3 (F)]- E VA ^> 
j [fc-j'|<3 

C^^E^-s^^-aW-A, Y, VA ^' 

3 \k-j\<3 

where (using frequency localization in dyadic annuli and Plancherel formula) 

llCailli. < E IN£« 2" 2J \\v\\ 2 w ^ 2 2 > Y H A HI^ S Hplli« U y Uci+* 1Mb • 

j |fc-i|<3 
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On the other hand, since Aj Yl\k-j\<3 = Aj, we have 

C22 =J2 s j-z( v ) s j-z(.P) VA i u = s i- 
3 

We thus end up with 

(2.38) II = T P T V -Vu + Ru, Re K. 
It follows from (2.32) and (2.38) that 

(2.39) (d t + V ■ V)T p u = T p {d t + T v ■ V)u + T dtP+v .^ p u + Ru, R e K. 

The symbolic calculus shows that TQ tp+ y.\j p 6 1Z, which proves (2.31) and concludes 
the proof of the first step. 

Step 2 : Higher order paraproducts. We now assume that p(t, x, £) = a(t, x)h{Q 
where h(Q = \£\ m h(£/\Z\) with h G (^(S^ 1 ). Then, directly from the defini- 
tion (2.1), we have T p = T a ip{D x )h{D x ) where ip satisfies (2.2). We have 

[T p , d t + T v -V] = [T a ,8 t + T V -V] ^(D x )h(D x ) + T a [ip(D x )h(D x ),d t + T v ■ V] . 

The norm from H m to L 2 of the first term in the right-hand side is estimated by 
means of the previous step by 

if||a||i,°o || VH^i+e + \\d t a + V • Va||z,°o ||F|| Loo , 

while the norm of the second term simplifies to T a [ifi(D x )h(D x ),Ty • V] and is easily 
estimated using (2.4) and (2.5) by 

||o||l°° ll^ll^i+E (II^-IIl 00 + ll^^ ls <i - 1 

Step 3 : Paradifferential operators. Consider an orthonormal basis (h u ) v ^* 
of L 2 (S a!_1 ) consisting of eigenfunctions of the (self-adjoint) Laplace-Beltrami op- 
erator, A w = Agd-i on L 2 ^- 1 ), i.e. AJi v = \ 2 Ji v . By the Weyl formula, we know 
that \ u ~ cud. Setting h v {£) = |£| m h v {u), w = £/ |£|, £ 7^ 0, we can write 



p(t,x,£)= a v (t,x)h v (£) where a u (t,x)= / p(t,x,ui)h„(u) 

— /SJd-1 



dw. 



Since 

\l k a u (t,x)= A k w p{t,x,u)~h v {uj)duj, 

we deduce 

_ 3d_i 

(2.40) sup \\a y {t)\\ Lac < C\ v 2 A4ff(p). 
Moreover, there exists a positive constant K such that, for all v > 1, 

(2.41) ||^|| L oc < CA? 1 . 
Now we can write 

|| [<9 t + V • V, T p ] u|| L2 < 2 || + V ■ V, T^] n|| L2 . 

So using the estimates obtained in the previous steps for every v > 1 and the 
estimates (2.40)-(2.41), we obtain (2.31), since the sum 

*^+l,-(¥ +1) 



is finite. This completes the proof of the lemma. □ 
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We have also a Sobolev analogue of Lemma 2.16 which can be proved similarly. 

Lemma 2.17. Let s > 1 + d/2 and V G C°([0, T]\ H s (R d )). There exists a positive 
constant K such that for any symbol p = p(t,x,£) which is homogeneous in £ of 
order m G R and all u G C°([0, T]; i? s+m (R d )), 

||[r p ,^ + TvV]n(t)||^ (Rd) 

< tf{Mo» + + v ■ Vp) [|F(t)|| LOO } imt)||^ +m(Rd) . 

2.5. Parabolic evolution equation. Consider the parabolic evolution equa- 
tion 

d z w + \D X \ w = 0, 

where z G R and x G R rf . By using the Fourier transform, one easily checks that 

(2.42) sup \\w(z)\\ Hr + ([ \\w(z)\\ 2 i dz)^ <K|K0)|| Hr) 
and, for r G]0, +oo[\N, 

(2.43) lh(*)[lc°([o,l];Cr) <^IK0)|| Cr . 

The purpose of this section is to prove similar results when the constant coefficient 
operator \D X \ is replaced by an elliptic paradifferential operator T a of order 1 with 
regularity C* in x for some p > 0. 

2.5.1. Tangential paradifferential calculus. Given I C R, Zq G I and a func- 
tion = ip(x, z) defined on H d x /, we denote by (p(zo) the function x h-> y?(x, Zo). 
For I C R and a normed space E, {p £ C°(I\E) means that z h-> </?(z) is a contin- 
uous function from / to i?. Similarly, for 1 < p < +oo, cp G L P Z (I\E) means that 
z h-> ||y?(z)|| E belongs to the Lebesgue space L P {I). We endow the spaces C®(I\E) 
and (J; 2?) with the usual norms. 

In this section, when a and u are symbols and functions depending on z, we still 
denote by T a u the function defined by (T a u)(z) = T a ^u(z) where z G L is seen as a 
parameter. We denote by r™(R d x I) the space of symbols a = a(z;x, £) such that 
z i — y a(z; •) is bounded from / into the space r™(R rf ) introduced in Definition 2.3. 
This space is equipped with the semi-norm 



(2.44) .M™(a) = sup sup sup (1 + |£|) |ahm d?a(z; ., £) 

*ei | a |<M +p+ i iei>i/2 



WP>°°(R d ) 



2.5.2. Estimates in Sobolev spaces. Given (j, G R we define the spaces 

(2.45) , 

= Li(/;^(R d )) +L2(/;^-2(R d )). 

Proposition 2.18. Let r £ K, p £ (0, 1), J = [z , zi] C R and let p G T^R ' x J) 
satisfying 

Rep(z;x,£) > c|£| , 

/or some positive constant c. Then for any f G y(J) and wo G H r (~R d ), there 
exists w G X r (J) solution of the parabolic evolution equation 

(2.46) <9 2 u; + T p w = /, w\ z=ZQ = w , 
satisfying 



\ w \\x r (j) 



< K^\\w \\ H r + ||/||yr (J) } , 
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for some positive constant K depending only on r,p,c and A4 p (p). Furthermore, 
this solution is unique in X S {J) for any s G R. 

PROOF. Let r G R. Denote by (•, -)ur the scalar product in H r (R d ) and chose F\ 
and F% such that / = F\ + F2 with 



llFillzKW + ll^ll^^-i^ll/ll^) 



+ 5, 5>0. 



Let us consider for e > the equation 

(2.47) d z w £ + e(-A + ld)w £ + T p w £ = f, w £ \ z=ZQ = w . 

Then standard methods in parabolic equations show that for any Z\ > zq, this 
equation have a unique solution in 

C°([z , 21]; H r (R d )) n L 2 ((z , zO; /T+ 2 (R d )) 

(here we only used that T p is a Sobolev first order operator). To pass to the limit e —> 
0, we need to establish uniform estimates with respect to e. Taking the scalar 
product in H r , directly from (2.47), we obtain 

" U7e ^' ) "^ r + £ ((- A + ld ) w e(z),w £ (z)) H r +Re{T p ( z) w £ (z),w £ (z)) H r 

< WF^Wnr \\w e (z)\\ Hr + 11^)11^1 \\w £ (z)\\ Hr+i . 

It follows from Garding's inequality (see [43, Section 6.3.2]) that there exist two 
constants C±, C% > depending only on A4 p (p) such that for any u G H r , 

Re{T p{z) u(z),u(z)) H r > Ci ||n(z)||^ r+ i - C 2 H^OH^i^. , 

for each fixed z G J. Therefore, we obtain 

2^ H^( z )ll^ + ^<(" A + ld)w £ (z),w £ (z)) H r + Ci ||w £ (z)||^ r+ i 

< ||Fi(z)||^ H^WIIfl, + II^WH^-i K(*)|| ^ + C 2 11^(2)11^.1^ • 
Integrating in z we obtain that, for all z G [zo, zi], 

^0) : = £ - |K(z )Hh''} + £ f z \\We{z')\\ 2 H r+ldz 

is bounded by 

B ■= II^IIlHJ;^) HIL-^) + M^.^-*) \\We\\ L2{J . Hr+i) 

+ C2||u%|| 2 ,i-p . 



Since 



(2.48) B < 4 \\F x f Ll{J . Hr) + - \\w £ \\ 2 LOO{J . tHr) 
and since 

(2.49) - HWellioojj.^) - g ||w £ (z )||^ + Ci Ikell^^.^jj < sup A(z) < B, 
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one can absorb the second and fourth terms in the right-hand side of (2.48) by the 
left-hand side of (2.49) to obtain (notice that the constants C, C\ are uniform with 
respect to e > 0), 



(2.50) \\w e \\ 2 LX{J . Hr) + d Wl^, ^ < \\M\ 2 H r 

+ ^(ll^illW) + 11*1.11^-4, + Kll^^^j). 

Finally, to eliminate the term ||u> F l| 2 , i- P in the right hand side of (2.50), it is 

" U L2(J;H r+ ^ E ) y ' 

2 

enough to notice that the left hand side controls by interpolation c \\w F \\ , i_„ , 

for some p > 2, hence by Holder in the z variable, there exists k > (depending 
only on p) such that if \zq — z±\ < K, we have 

we consequently obtain 

|2 



(2.51) \\w £ \\ L ^ {J . Hr) + Ci \\ w e\\ L2{J . H r + i } 

<2\\w \\ 2 Hr + C(\\F 1 f Ll(J . Hr) + \\F 2 



2 , 77! 1 1 2 ; n 1 1 2 \ 

L 2 (J;H r --Z)' 



We can now iterate the estimate between z$ + k and zq + 2k, ... to get rid of the 
assumption \z\ — Zq\ < K (and of course the constants will depend on z\). By using 
the equation, we obtain now that sequence (w £ ) is bounded in 

L°°( J; H r ) n L 2 (J; H r+ ^) n C^J; F r - 2 ). 

It follows from the Banach-Alaoglu theorem that, up to a subsequence, (w £ ) con- 
verges in the sense of distributions to w G X r (J), which satisfies the equation 
d z w + T p w = f. Then d z w belongs to Y r (J) which implies that w belongs to 
C°([zo, Zi]; H r (R )). Moreover, by the Ascoli theorem, up to a subsequence, (w £ ) 
converges in C°([zq, zi]; Hj^) for some \i > 0. Since w £ \ z= o = wq we obtain that 
w\ z= o = wq, which completes the existence part in Proposition 2.18. The proof of 
uniqueness follows the same steps and we omit it. □ 

2.5.3. Estimates in Holder spaces. The next proposition is the analog in Holder 
spaces of the previous result. Here we follow a similar strategy previously used 
in [48, 40, 5, 1]. 

Proposition 2.19. Let p g (0,1), J = [z ,zi] c R, p G T^IV 1 x J) with the 
assumption that 

Rep(z;x,£) > c|£| , 
for some positive constant c. Assume that w solves 

d z w + T p w = F 1 + F 2 , w\ z=Z0 =w - 
Then for any q G [1, +oo], (ro, r) G R 2 ro < r, i/ 

u; G L°°(J;C: ), Pi G ^(JjCJ), F 2 G L q {J ; cl~ 1+ ~ q+5 ) with 6 > 0, 
and wo G C£(R d ), u;e /iai>e u; G C°(J;C£) and 

\c°(J:Cr) < K { \\ w o\\ci + \\ F i\\l^(.J;CI) + 11*2 II r-l+i+S + IIHI £«>( J;C£°) f ' 

/or some positive constant K depending only on tq, r, p, c, 5, g and _A/fp(p). 
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Proof. For this proof, we denote by K various constants which depend only 
on r$,r,p,c and A4 p (p). Given y G J introduce the symbol e = e(y,z;x,£) defined 
by 

ry 

e(y,z;x,0 = exp(- / p(s; x, £) ds) (ze[z ,y]). 

J z 

This symbol satisfies d z e = ep, so that 

d z (T e w) = (T ep - T e T p )w + T e F, F = F l + F 2 . 
Integrating on [zo,y] the function ■f^T e ^ y ^ z ^ x ^w{z) 1 we find 



2.52) T lW (y) = T e{z=zo w + (T e F)(z)dz + (T ep - T e T p )(z)w(z) dz. 



2(1 



-0 



(Notice that the paraproduct T\ differs from the identity only by a smoothing oper- 
ator.) Introduce G(y) = T e \ z=ZQ WQ + j^ Q (T e F)(z) dz and the operator R defined on 
functions u: J ->• C m (R d ) by 



ry 

(Ru)(y) = (T ep - T e T p )(z)u(z) dz 
J zo 



so that T\w = G + Rw. Now, by a bootstrap argument, to complete the proof 
we see that it is enough to prove that the function G belongs to L°°(J;C^) and 
that R is a smoothing operator of order —a for some a > 0, which means that R 
maps L°°(J; C*) to L°°(J; Cl +a ). Indeed, by writing 

«, = (I + i? + • • . R N )G - R N+1 w, 

and choosing N large enough, we can estimate the second term in the right-hand 
side in L°°(J; C£) by means of any L°°( J; C£°)-norm of w. 

In the analysis, we need to take into account how the semi-norms M~ m (e(z)) (see 
Definition 2.3) depend on z. Then the key estimates are stated in the following 
lemma. 

Lemma 2.20. For any m > there exists a positive constant K depending only on 
supj Mp(p(-; x, £)) such that, for all y £ (0, — z{\ and all z £ [0, y), 

K 



(2.53) 



m: 



(e(z)) < 



(y - z) r 



This follows easily from the assumptions p G T p , Rep(s;x,^) > c|£| , and the ele- 
mentary inequalities (valid for any a > 0) 

(y-zr|e| a exp((z-y)|e|)<l. 

By using the bound (2.53), applied with m = 0, it follows from the operator norm 
estimate (2.5) that, for any z <y and any function / = fix), we have 

(2.54) \\T eM f\\ Cl <MS(e(y,z))\\f\\ Cl 
This implies that 



< K 



ry 

T e \ z=zo w + (T e F 1 )(z)dz 

J za 



< K [poller + K II-^iIIl 1 



L°°(J-CZ) 



On the other hand, by using the bound (2.53), applied with m = 1 — - — 5, we obtain 
that 



(T e F 2 )(z)dz 



-0 



< K 



i°°(J;Cr) 



Z(l 



\y - z \ 



■\\F 2 (z)\\ c r-mdz 
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which implies by Holder inequality that 

W G \\L°°(J;Cr) < K \\ W o\\ci + K \\ F A\l^(J;CZ) + H^ 2 ' 



r-l+i+4 • 

i 9 (J;C» 9 ) 



It remains to show that R is a smoothing operator. To do that, we first use the 
operator norm estimate (2.5) (applied with (m,m',p) replaced with (— m, l,p)) to 
obtain 

\\(T ep -T e T p )(z)\\ ct ^ ci+m - 1+P <M; m {e(z))M l p {p(z)). 
Taking m = 1 — p/2, it follows from the previous bound and Lemma 2.20 that 

\\(T ep - T e T p )v(z)\\ ct+p/2 < j- * \\v(z)\\ cl . 

Since < m < 1 we have Jq (y — z)~ m dz < +oo and hence 

rv 

(2.55) \\Ru{y)\\ c t+ P /2 < / \\{T ep -T e T p )u(z)\\ c t+ p/ 2 < K\\u\\ Lx{ j. cl) , 

J 

which completes the proof. □ 



3. The Dirichlet-Neumann operator 



A notable technical step in the analysis of the Euler equation with free surface con- 
sists in describing the Dirichlet-Neumann operator in domains with free boundaries 
of limited regularity. Here we shall prove some results about elliptic regularity which 
complement previous works (see Propositions 3.18 and 3.19). To do this we shall 
use a paradifferential approach. 

3.1. Definition and continuity. In this paragraph we recall from [1] the def- 
inition of the Dirichlet-Neumann operator under very general assumptions on the 
bottom. One of the novelty with respect to our previous work is that we first clarify 
the regularity assumptions: assuming only that r\ G W 1,oc (R d ) (which is satisfied 
if 77 £ H s+ ^(R d ) with s > l/2+d/2) and / G H?(R d ), we show how to define G{rj)ij) 
and prove that the map 

ip G H^(R d ) ^ G(r])ip G H~-i(R d ) 

is continuous. Our second contribution is to prove that the map r] 1— > G{rj) is 
Lipschitz (in a proper topology). Finally, we prove also that in some weak sense, the 
Dirichlet-Neumann operator thus defined is a local operator (see Theorem 3.9). We 
also refer to chapter 3 in [38] for an introduction to the analysis of the Dirichlet- 
Neumann operator. 

The goal is to study the boundary value problem 

(3.1) A„^ = 0, 0| s = /, cU| r = 0. 

See §1.1 for the definitions of f2,£,r. Since we make no assumption on T, the 
definition of <p requires some care. We recall here the definition of <f> as given in [1]. 

Notation 3.1. Denote by S 1 the space of functions u G C°°(Q) such that V x , y u G 
L 2 (ft). We then define S>q as the subspace of functions u G Q> such that u is equal 
to in a neighborhood of the top boundary S. 

PROPOSITION 3.2 ([1, Proposition 2.2]). There exists a positive weight g G L^ C (VL), 
equal to 1 near the top boundary ofQ, and a constant C > such that for all «£^o, 

(3.2) / g(x,y)\u(x,y)\ 2 dxdy <C |V Xty u(x,y)\ 2 dxdy. 
Jn Jn 
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Definition 3.3. Denote by H 1,0 (Q) the space of functions u on Q such that there 
exists a sequence (u n ) G £Fo such that, 

V x ,yUn ^x,yU in L 2 (O, dxdy) , u n —¥ u in L 2 (£l, g(x,y)dxdy). 
We endow the space H 10 with the norm 

\\ u \\ = \\Vx,v u \\L2(a.) ■ 

Let us recall that the space H 1 ' ^) is a Hilbert space (see [1]). For later purpose, 
we need to ensure that the functions having compact support in the x variable (at 
least near the surface S) are dense in H '°(0,). By regularizing the function rj (see 
Remark 3.7) it is easy to see that there exists 77* G C£°(R d ) such that n — > 77* 
and 

{(x,y) G ~R d x R ; r?*(x) < y < n(x)} C fl 

Lemma 3.4. TTie se£ 

T>o = |n G £> ! supp(n) n {(x, y); -77* + ^ < y < 77} is compact^ 
is dense in H 1,0 (Q, n ). 

Proof. Let ueT> , and ( G C°°(R) equal to for z < and to 1 for z > h/30. 
Then according to Proposition 3.2, we have ((y — rj^)u G Vq and (1 — ((y — r)*))u G 
Vq n Hq(Q) (where Hq(Q) is the usual Sobolev space). Let v n G Co°(Q) which 
converges to (1 — £(y — r)*))u in Hq(Q) (and hence in H We get that C(y — 

77*)?/ + % £ Do and converges to rx in i7 1,0 (f2). □ 

We are able now to define the Dirichlet-Neumann operator. Let / G if2(R rf ). We 
first define an H 1 lifting of / in $7. To do so let \o £ C°°(R) be such that Xo( z ) = 1 
if z > — \ and xo( z ) = if z < —1. We set 

M*,*) = Xo(z)e z{D ^f(x), x G R d ,z < 0. 
By the usual property of the Poisson kernel we have 

l|Va;,*W| i 2([_ lj0 ] xR d) - C 11/11 #2(Rd) • 

Then we set 

ip(x,y) = ipt(x,- — 7^-), (sc,y)€Q. 

This is well defined since Q C {(x,y) : y < rj(x)}. Moreover since the bottom T is 
contained in {(x,y) ■ y < r](x) — h}, we see that tjj vanishes identically near T. 

Now we have obviously = / and since V77 G L°°(H d ), an easy computation 
shows that ip G H 1 ^) and 

(3-3) IMIffi(n) < ^(1 + \\v\\w^)\\f\\ H ^ Rd y 

Then the map 

v i->- - / Vj,^ • V^u (iriiy 
./n 

is a bounded linear form on i2~ 1,0 (f2). It follows from the Riesz theorem that there 
exists a unique u G -£f 1,0 (fi) such that 

(3.4) Vtj G H 1,0 ({1), / V x .yU ■ V X)V v dxdy = — I ■ V X)V v dxdy. 

Jn Jn ' ~ 
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Then u is the variational solution to the problem 

-A x>y u = A X) yij_ inD'(O), u |s= 0, 3 n u | r = 0, 
the latter condition being justified as soon as the bottom V is regular enough. 

Lemma 3.5. The function (j) = u + ip constructed by this procedure is independent 
on the choice of the lifting function ip as long as it remains bounded in and 
vanishes near the bottom. 

Proof. Consider two functions constructed by this procedure, (f>k = u^+tpi,, k = 
1, 2. Then, by standard density arguments, since ip\ — ip2 vanishes at the top bound- 
ary X and in a neighborhood of the bottom T, there exists a sequence of func- 
tions ip n G C^°(f2) supported in a fixed Lipschitz domain Sled tending to ipi — ip2 
in Hq(Q) and hence also in H 1,0 (Q). As a consequence, ipi ~ 1 p2 £ H 1,0 (Q,) and 
the function <f> = <f>\ — fa is the unique (trivial) solution in H 1,0 (Q) of the equation 
A x , y 4> = given by the Riesz Theorem. □ 

Definition 3.6. We shall say that the function (ft = u + tp constructed by the above 
procedure is the variational solution of (3.1). It satisfies 

(3-5) JjV^l 2 dxdy<K\\f\\ 2 Hi{Rdy 

for some constant K depending only on the Lipschitz norm of n. 



Formally the Dirichlet-Neumann operator is defined by 

=n(x) 



(3.6) G{ri)i> = s/l + WW d n <t> \ y=v(x) = [d y 4> ~ Vr? • V<p] \ y __ 



3.1.1. Straightenning the free boundary. We begin by straightening the bound- 
ary. In this paragraph, we fix s > 5 + f ■ 

We shall assume here that one can find a function 77* such that 
(0 rh + ^H°°(R d ), 



(3.7) 



k (Hi) r C {(x,y) e O : y < ??*(x)}. 



Remark 3.7. Assume that we have a function depending smoothly on the time, 
r/ G C°([0,T],H s+ 2 (R d )), such that rj\ t=0 = 770 and satisfying condition (1.2) with 
7j (such as a solution in Theorem 1.2). Then one can construct 77* = 77* (a?) satisfying 
(z), (m) in (3.7) and for some T' <T 

Indeed set 

^(«) = ~^ + e~' /{Dx) Vo(x) 
where v > is chosen such that 1/ II 770 II „_i_i . < ttt- Then chose T" such that 

II lUUffS+l^d) — 10 

ll^-)-r/o|l W(Rd) <^ Vt€[0,n 

and write 

rj(t,x) - r]*(x) = r/(t,x) - 770(0;) + 770(2;) - e" u ^ Dx) rj + ~. 

25 



Then (ii)' follows from the estimate 
\\ m (x)-e-^rio\\ 



and (Hi) is a consequence of (ii)' ■ Indeed for t G [0,T'] we have 



r](t,x) -l]*(x)<~ + 



h 

< -7 + 



h 
< -. 



therefore 



4 \\y\\L°°([o,Ti]xR<i) ^ 1 1- iiyii iO c {[0)T / ])H s -5Rd)) - 2 

r c {(sc,y) : y < - ^} c {(x,y) : y < J/*(ac)}. 



In what follows we shall set 

' Qi = {(x,y) : x G R d ,??*(x) < y < Tj(x)}, 
n 2 = {(x,y) eO:y< rj*(x)}, 
= HiU fi 2 . 



(3.8) 
and 
(3.9) 



fii = {(x,z) : x G R d ,z G /}, I =(-1,0), 

2 = {(x, z) £~R d x (-oo, -1] : (a?, z + 1 + V*( x )) € ^2}, 

= Oi u o 2 . 



(3.10) 



Following Lannes ([39]), consider the map (x,z) 1— > p(x,z) from f2 to R ^ 1 defined 
by 

J p(x,z) = (1 + z)e 5z{ - Dx) ri(x) - ZTj*(x) if (x,z) G f^, 
[ p(x, z) = 2: + 1 + rj*(x) if (x, z) G f2 2 , 
where 5 is chosen such that 

S 1 1 77 1 1 H s+ 1 : = S is small enough. 

Notice that p is Lipschitz on Q. Moreover since s > 5 + f , there exists a constant 
C > such that 



(3.11) 



and 



l|v^|| LCO(/ ^_i (Rd)) <c(h|| HS+ i (Rd) + 



(3.12) 



\\VxP\\l?(I,H'(R<>)) < C (\\v\\ H s +i{T{d) + \\V*\\ H s +i{nd) 

from which, taking 5q small enough, we deduce 



(3.13) 



(i) d z p(x, z) > min(l, — ), V(x,z) G O, 



It follows from (3.13) (i) that the map (x,z) 1— )■ (x, p(x, z)) is a C 1 -diffeomorphism 
from O to f2. We denote by k the inverse map of p: 

(x, z) G O, (x, p(x, z)) = (x, y) (x, z) = (x, k(x, y)), (x, y) G tt. 
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We set 

<j){x,z) = <p(x,p(x,z)). 

Then we have 

(d y </))(x, p(x, z)) = (Ki4>)(x, z), (X7 x (j))(x, p(x, z)) = (A 2 (j))(x, z) 



Ai = ^-d z A 2 = V x -^d z . 

OzP OzP 



(3.14) 

If <f) is a solution of A x ^(j) = in Q then <j> satifies 

(A? + A^)0 = O in ft 

This yields 

(3.15) {ad 2 z + A x + b ■ V x d z - cd z )4> = 0, 
where 

(3.16) a:= 1+ ,^*f , b:=-2^, c := J- (ad 2 z p + A x p + b ■ V x d z p) . 

[d z pY d z p d z p y 

It will be convenient to have a constant coefficient in front of d 2 v. Dividing (3.15) 
by a we obtain 

(3.17) (d 2 z +aA x + !3-V x d z - 7 9^=0, 
where 

< 3 - 18 » a; =TW' ^=- 2 TT^' t-&®p+<*«+ W)- 

In the coordinates (x,z), according to (3.6) we have 

(3.19) G(rj)^ = U\ z=0 , U = A 1 $-V x p-A 2 $. 
The following remark will be useful in the sequel. We have 

(3.20) 8 Z U = -V x ({d z p)k2$). 
Indeed we can write 

d z ll = d z hi4> - V x d z p ■ A 2 (f> - V x p ■ d z A 2 4> 

= {d z p)K\4> - V x 8 zP • A 2 ^ + (d zP )(A 2 - V X )A 2 4> 
= (d z p)(Ai + Al)0-V x ((d z p)A 2 $). 

Since (A? + A$)<j> = we obtain (3.20). 

3.1.2. Continuity of the Dirichlet Neuman operator. 

Theorem 3.8. Let r\ £ Ty 1 ' oc, (R ci ), / G H^(R d ). In the system of coordinates (x, z) 
defined above, the variational solution of (3.1), <p, satisfies 

(3.21) 0GC 2 °([-l,O];F5(R d ))nCi([-l,O];F-^(R d )). 

As a consequence, the map 
(3.22) 

ip G H l 2(R d ) ^ G{rj)ip = ^YTWWdn^ \ y=r)(x) = [d y <f> - Vt? • V0] | y=r){x) 

= ((l + |Vr?| 2 )aJ- Vn-V4>)\ Z=Q 
is well defined. It furthermore satisfies 
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The Dirichlet Neumann operator is also weakly continuous. This fact will be used in 
Section 6 to prove existence of solutions when passing to the weak limits on weakly 
convergent sequences of suitably regularized systems. 

Theorem 3.9. Assume that (?7 n )neN and (V'n)neN are two sequences such that 

i) the sequence (rj n ,' l Pn)ne'N is bounded in W 1,co (R d ) x i/2(R d ) 5 
ii) there exists n G W 1,00 (H d ) such that r\ n converges strongly to rj in W l ^(R d ), 

Hi) there exists i\) G ff2(R d ) such that (ip n )neN converges weakly to ip in //a(R d ) ; 
iv) there exists 77* G W l '°°(Yl d ), h > such that 

h h 
7](x) - - > V*( x ) > v( x ) ~ h,rj n (x) - - > V*( x ) > Vn(x) -h Vx G R d . 

Then the sequence (G(r] n )tp n ) is bounded in .ff _ 2(R d ) and converge weakly to G(n)^. 

Let us also state the second basic strong continuity of the Dirichlet Neumann op- 
erator. Notice that the map rj 1— > G{rf) is strongly non linear and hence strong 
continuity do not imply weak continuity. 

Theorem 3.10. There exists a non decreasing function T : R + — > R + such that, 
for all r)j G W ,1 ' DO (R <i ), j = 1,2 and all f G H*(B*), 

||(G(m) - G(r] 2 ))f\\ H _i < J r (\\(m,V2)\\w i ^xw i ^)\\m - m II w 1 ^ \\f\\ H i- 

Remark 3.11. We shall only prove Theorems 3.9 and 3.10 as the choice rj n = r], tp n = 
tp in Theorem 3.9 implies Theorem 3.8. On the other hand, the fact that we can pass 
to the limit in G(rj n )ip n under convergence assumptions on (r] n ,ip n ) which are only 
local in space, shows that, in a very weak sense, the Dirichlet-Neumann operator is 
a local operator. 

PROOF of Theorem 3.9. Let M > be such that for all n G N, 
(3.23) 

H^ 71 H_H"3 (R d ) ^^ff3(R d ) ll ? '' 1 llw /1 ' 00 (R' d ) ll^ll VK 1 ' 00 (R d ) \\ r l*\\w 1 - ao (R d ) 

<M . 

Our purpose is to prove that G(rj n )ip n is well defined and bounded in H~2(R d ) by 

(hence uniformly bounded) and converges weakly to G(n)ip in ^"^(R^). We shall 
proceed in several steps. 

Step 1: preliminaries. We start by straightening the boundaries of the domains 
f2 n and £1 using the previous section. We recall that 

n n = {(x, y)GO:y< rj n (x)}, = {(x, y) G O : y < rj(x)}. 

For this purpose we use the diffeomorphisms given by p n (constructed with rj n ) and 
p given by (3.10) and we shall use the vector fields A",Aj,j = 1,2 described in 
(3.14) and we set 

A" = (A?,A£) A = (A l5 A 2 ). 

We construct now a H 1 - extension of ip n . Let x £ C' X) ('R),x(z) = 1 if z > 
and x( z ) = if z < — 1 and set 

(3.24) Mx, z) = x(z)e z ^H n {x), ±(x, z) = X (z)e z ^^(x). 

Then Mx,z) G H l (K d x /) and \fy\ H i {RdxI) < C|lV>n||^i (Rd) < CM 
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We make the same construction for ip. Then it is easy to see that the sequence (ifi n ) 
converges in H 1 ^) to tp, Then we set 

(3.25) (fi n = u n + ipn, (fi = U + *fi_. 

According to (3.5) and the assumptions on r] n and rj we see easily that this implies 

(3.26) \\V x>z u n \\ L2 ^ ) <M 2 , VneN. 

Then (u n ) is a bounded sequence in H 1,0 (fl)) and therefore that, up to a subsequence, 
it converges weakly in this space to u. 

Step 2: passing to the limit for the variational solutions. Setting X = {x, z) € f2 
the variational formulation for u n reads 

(3.27) f A n u n (X) ■ A n ((X)J n (X)dX = f A n MX) ■ A n ((X)J n (X)dX 
Jn Jn 

for all C £ Co°(^)) where J n (X) = \d z p n {X)\. We now want to identify the limit. 
We have the following Lemma. 

Lemma 3.12. For all ( G V (Q) we have 

lim / A n u n (X) ■ A n ((X)J n (X)dX = [ Au(X) ■ A((X) J(X)dX, 
n->+oo J n J n 

lim / A n $ n {X) ■ A n ((X)J n (X)dX = [ A^(X) ■ A((X) J(X)dX. 



Corollary 3.13. The function u(x,y) = u(x,k(x,u)) is the variational solution 
in H 1 '°(Q) of the problem —A X:V u = A Xjy ip and u n converges weakly in this space 
to u. 



(3.28) 



Proof of Lemma 3.12. Notice that 

j A n -A = p n d z , supp/3 n C {{x,z) :x€R d ,zG (-1,0)} and 

[ ll/3n.lU°°(i<r) < J r (\\v\\w i - a: '(R d ))\\'nn - 'nWw^iK) 

Then we can write 

A n u n ■ A n C ■ J n - Au ■ AC, ■ J = A x + A 2 + A 3 + M, 

A x = {A n - A)u n ■ A n C ■ J n , A 2 = Au n ■ A n ( ■ (J n - J), 

A 3 = Au n ■ (A n - A)C -J, A 4 = A(u n - u) ■ A( ■ J. 

It follows from (3.28) that we have 



(3.29) 



J_Ax{X)dX < C , (M )||77n-r/|| w ,i > oc (x) ||^5 n || i2(S) ||VxCI| L2( n))- 



The same estimate holds for the term coming from ^3. Moreover since \\J n — 
J\\l°°(K) < F(Mo)\\rj n — vWw-'^iK) we have for A 2 the same estimate as (3.29). 

Eventually since (u n ) converges to u in the weak topology of H ,0 (Cl) we obtain 



lim / A<x(x, y)dxdy = 0. 



□ 
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Step 3: taking traces. Notice that we have 
(3.30) ((A™) 2 + (A£) 2 K = 0, ((Ax) 2 + (A 2 ) 2 )5 = 0, 

and 

j G{r] n )tp n = (A™ - V x p n • A2)5„| 2 =o =: U n \ z=0 , 
{ G{r,)i) = (Ai - V x p ■ A 2 )u\ z=0 =: l/|*=o. 



Since p n converges to /? in W^^°(Q) the sequence (?7 n ) converges weakly to U in 



(3.31) 
iii 

L 2 (tt). Now using (3.20) we obtain 

d g U n = -V x ({d z p n )A%U n ). 

By the same way we have 

d z U = -V x {(d zP )A 2 v). 

Since V X)X p n — >■ V x>z p in L°°(R d x I) and u n — > u weakly in H 1 (R d x I), the 
sequence (d z ll n ) converges to d z U weakly in L 2 (I, H~ 1 (R d )). 

Now we use the following well known interpolation lemma. 

Lemma 3.14. Let I = (—1,0) and consider u G L 2 (I , L 2 (R d )) such that d z u G 
L 2 (I,H- 1 (R d )). Then u G C°([-l, 0],if"2(R d )) and f/iere exisfe an absolute con- 
stant K > suc/i £aa£ 

l|n|l C0([-l,0];/f-3(Rd)) " K (ll n lli 2 (/;i 2 (R d )) + \\ d M\ L l(I;H-HK d )))- 



It follows from this lemma that the sequence (U n \ z= o) is bounded in H 2 (R 1 
by •? r (l|%llw 1 .°°)ll^n|l ir i, Rd , and converges weakly in H~^(R d ) to U\ z=0 , which 
completes the proof of Theorem 3.9. □ 

PROOF of Theorem 3.10. We use the notations introduced in §3.2.1. Namely, 
for j = 1, 2, we introduce Pj(x, z) and Vj(x, z) defined by (3.10), 

Pj(x,z) = (l+z)(e Sz{Dx) r ]j )(x)-zrj^ if x G R rf ,2 G / := (-1,0) 

Pj(x, z) = z + 1 + r/* , if (x, z) g ri 2 

Notice that we have the following estimates 



(3.32) 



00 dzPj > min(l, -), (x,z) G 



(") IIV^-II^^ < C(l + ||7?i|| HS+ i (Rd) ) 

(iii) \\V x ,z(fil - P2)||L°°(/,L° c (R d )) < C\\VI - r l2\\w' L ' ca ('R d ) ■ 
Recall also that we have set 

(3.33) A \ = ^-d z , Al = V x -^d z . 

o z pi d z pi 

It follows from (3.32) that for k = 1, 2 we have with W 1 ' 00 = W 1 ' 0O (R d ) J 

(i) A\- Al = f3jd z , with supp/3 fc C R d x I, where / = (-1, 0), 
l|A:|lz<*>(/xR d ) < •^(IK??l,%)||Tyi>°°xVl/ 1 .°°)ll ? ?l - ^IIh/ 1 '- 



(3.34) 



Then we set cpj(x,z) = (f>j(x, pj(x,z)) (where ^ x , y (pj = in Qj,cj)j\z = f) and we 
recall (see (3.19)) that 

(3.35) G(r]j)f = Uj\ z=0 , U, = A{^ - V x p 3 ■ A^. 
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Lemma 3.15. Set I = (—1,0), v = (pi — 4> 2 , and A J = (Aj, A2). There exists a non 
decreasing function T : R + — > R + such that 

(3.36) IIA^II^.^Rd)) < J r (\\(rji,r) 2 )\\w 1 '°°xW 1 <°°)\\'ni ~ m\\w^ 



1 . 

Let us show how this Lemma implies Theorem 3.10. According to (3.35) we have 

U 1 -U 2 = (1) + (2) + (3) + (4) + (5) where 
(3.37) (1) = AK (2) = (Ai-A?)0 2 , (3) = -V x ( Pl - p 2 )A 2 ^ 1 

(4) = -(V xP2 )Alv, (5) = -{V xP2 ){A\ - A|)&. 

The L 2 (I,L 2 (H d )) norms of (1) and (4) are estimated using (3.36). Also, the 
L 2 (I, L 2 (R, d )) norms of (2) and (5) are estimated by the right hand side of (3.36) 
using (3.34) and (3.5). Eventually the L 2 (I, L 2 (K d )) norm of (3) is also estimated 
by the right hand side of (3.36) using (3.32) (Hi) and (3.5). It follows that 



(3.38) \\Ui - U 2 \\ L 2 {I:L 2) < J r (||(7?i,r? 2 )||iyi,-»xvy 1 ^)ll ? ?i ~ mWw^WJW^- 
Now according to (3.20) we have 

(3.39) d z (U! - U 2 ) = -V x (d z ( Pl - p 2 )A^i + (d zP2 ){K\ - A 2 )fa + (d zP2 )A 2 v). 
Therefore using the same estimates as above we see easily that 

(3.40) \\d z (Ui - U 2 )\\ L 2( I)H -i) < J r (||(?7i,r ? 2)|| w i.°°xiyi.°°)lhi - rj 2 \\wi>°° \\f\\ H ^- 
Then Theorem 3.10 follows from (3.38), (3.40) and Lemma 3.14. □ 

Proof of Lemma 3.15. We use the variational characterization of the solu- 
tions U{. First of all we notice that §\ — 4> 2 = u\ — u 2 =: v. Now setting X = (x, z) 
we have 

(3.41) / A l Ui ■ A l 6 J,dX = - ! A*/ • A l e Xi dX 
Jn Jn 

for all 6 £ H^°(h), where Ji = \d zPi \. 

Taking the difference between the two equations (3.41), using (3.32) and setting 
8 = v = u\ — u 2 one can find a positive constant C such that 



6 

/ \A l v\ 2 dX <CJ2A k 
Jn k=i 



where 



A x = ! I (A 1 - A^tZaHA^I h dX, A 2 = I \{A l - A 2 )i>||A 2 5 2 | J x dX, 
Jn Jn 

A 3 = f |A 2 2 2 ||A 2 t>| I Ji - J 2 | dX, M = t I (A 1 - A 2 )j\\A 1 v\ J x dX, 

Jn Jn 

A 5 = f I (A 1 - A 2 )v\\A 2 J\ Ji dX, A 6 = f ^J^vW^-MdX. 

Jn Jn 

Using (3.34), (3.5), (3.32) we can write 

\A\\ < ||/3||Loo(/xR d )ll^l|lL°°(/xR d )ll^2^2||L2( /><R <i)||A 1 w|| L2 ^ ) 

(3.42) 

< FiWimim^w^xw^Wm -m\\w^\\f\\ H ^\\^ v \\ L *{ny 

Since A 1 -— A 2 = -^j^A\ the term A 2 can be bounded by the right hand side of (3.42). 
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Now we have \\J% — J2\\L°°(ixn, d ) — Cll^i ~~ r l2\\w 1 -°°(R. d ) an d 

l|A 2 v|| i2(S) < J c "(||(?7i,7/2)||h/i>-xw 1 .°°)II a1 ^II L 2 { ^ ) - 

So using (3.5) we see that the term ^3 can be also estimated by the right hand side 
of (3.42). To estimate the terms to Aq we use the same arguments and also (3.3). 
This completes the proof. □ 

Let us finish this definition section by recalling also the following result which is a 
consequence of [1, Lemma 2.9]. 

Lemma 3.16. Assume that —\ < a <b < — | then the strip S a ^ = {(x,y) G R rf+1 : 
ah < y — rj(x) < bh} is included in Q and for any k > 1, there exists C > such 
that 



< 



1 



3.2. Paralinearization of the Dirichlet-Neumann operator. In the case 
of smooth domains, it is known that, modulo a smoothing operator, G{rf) is a pseudo- 
differential operator with principal symbol given by 



(3.43) X(x, := ^(l + lV^x)! 2 )^! 2 -^^)^) 2 - 

Notice that A is well-defined for any C 1 function r\. The main result of this sec- 
tion allow to compare G{r\) to the paradifferential operator T\ when rj has limited 
regularity. Namely we want to estimate the operator 

R(V) = G( V )-T X . 

Such an analysis was at the heart of our previous work [5] [1, Proposition 3.14] for 
"smooth domains" (77 G H s+ z,s > 2 + |). Here we are able to lower the regularity 
thresholds and precise the dependence of the constants (though this dependence will 
not be essential in the present paper; it is important in [3]). 

Theorem 3.17. Let d > 1 and consider s, r G R such that 

3 d 
S> 4 + 2' r>L 

Consider 77 G H s+ ^{K d ) n d (R d ) and f £ H s (R d ) n Cl(R d ), then 

R(7])f G H s -^(R d ). 

Moreover 

(3.44) \\R(r,)f\\ HS _ h < Hh\\ HS+i , \\f\\ HS ) {l + ||r7|y + , 

/or some non- decreasing function T : (R + ) 2 — > R + depending only on s and r. 

The proof of Theorem 3.17 is given in Section 3.3.4. where we also state two corol- 
laries of the method used to prove Theorem 3.17. The following result, which we 
think is of independent interest will be proved in §3.3.2. It complements previ- 
ous estimates about the Dirichlet-Neumann operator by Craig-Schanz-Sulem [28], 
Beyer-Giinther [13], Wu [54, 55], Lannes [39]. 

Theorem 3.18. Let d > 1, s > I + | and \ < a < s + \. Then there exists 
a non- decreasing function T : R + — > R + such that, for all i] G i^ s+ 2(R d ) and 
all f G H a (R d ), we have G(r})f G H a ~ 1 {K d ), together with the estimate 

(3.45) \\GW\\b-W) ^ HM^^) Wf\\H° m ■ 
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We also prove error estimates (see §3.3.3). 

Proposition 3.19. Let d>l and s > | + f . For any \ < a < s and any 

1 1 d 

0<£<-, e < s , 

~ 2 2 2 

there exists a non- decreasing function T : R + — > R + such thatR(rj)f := G(rj)f—T\f 
satisfies 

\\R(v)f\\ H "-i+H-R d ) - ^(II^H^+^Rd)) H/H^(R d ) • 

3.2.1. Flattening of the free boundary. We shall use the diffeomorphism p defined 
by (3.10) which satisfies the estimates (3.13). Also we have 

v\ z=0 = 4>\y=n(x) = f, 

and 



G(n)f = (A lV - Vp ■ A 2 v)\ z=0 = ( 1+ J VP| d z v -Vp-Vv 

OzP 



2 

2=0" 



3.2.2. Elliptic regularity in Sobolev spaces. In this paragraph we state elliptic 
estimates for the solution v of (3.17) with boundary data / on z = 0. For later 
purpose, we will consider the non-homogeneous case. This yields no new difficulty 
and will be useful later to estimate the pressure (see Section 4.3). We thus consider 
the problem 

(3.46) d 2 z v + aAv + P-Vd z v-jd z v = F , v\ z=0 = f, 

where / = f{x) and Fq = Fq(x,z) are given functions. Recall that for p G R, the 
spaces X^(I),Y^(I) are defined by (see (2.45)): 

X»{I) = C°(I; H' M {R d )) n L 2 Z (I; H^* (R d )), 

Y»{I) = L 1 z (F,H^(R d )) + L 2 z (F,H^-^{K d )). 

Lemma 3.20. Assume that a > | . Then the space X U (I) is an algebra. Moreover 
if F : C N — > C is a C 00 '-bounded function such that F(0) = one can find non 
decreasing functions J 7 , T\ from R + to R + such that 

\\F{U)\\x°(i) < H\\U\\L°°(ixB.d))W\\x*(i) <?i(\\U\\ X '{i))- 

Proof. It is known that H a (R d ) is an algebra and so is C°(L; H a (R d )). On 
the other hand, according to (2.4) and (2.9), we have 



\ U V\\ . < \\u\\ T°o(Tid\ \\v\\ .. + \\v\\ Too(-Dd\ \\u\\ _,1 

H + 2(R d ) \ I H + 2(R d ) \**" > H + 2(R d ) 



□ 



With these notations, we want to estimate the X^-norm of V XtZ v in terms of 
the i7 CT+1 -norm of the data and the y°"-norm of the source term. An important 
point is that we need to consider the case of rough coefficients. In this section we 
only assume that r] 6 iJ s+ 2(R d ) for some s > 1/2 + d/2. An interesting point is 
that we shall prove elliptic estimates as well as elliptic regularity (in other words, 
we do not prove only a priori estimates). Our only assumption is that v is given 
by a variational problem, so that one has an estimate for the i? 1 -norm of v, and 
consequently 

(3.47) HV^dl i „ , < +oo. 

V I II x,z Hx-3([_1,0]) 
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Remark 3.21. In the case where v(x,z) = 4>(x , p(x , z)) with <f> the variational solu- 
tion of 

&x, y <f> = 0, <j>\y=n = f, d n (j) = on T, 
then (3.21) shows that v satisfies this assumption (3.47). 

Proposition 3.22. Let d>l and 

1 d 1 1 

s > — | — , — < a < s . 

2 2' 2 ~ ~ 2 

Consider f G i2" CT+1 (R d ), Fq G y <J ([— 1,0]) and v satisfying the assumption (3.47) 
solution to (3.46). Then for any zq G (—1,0), V x ,z v £ X a ([z$, 0]), and 

ll v ^ll^([ ZO ,o]) ^ HU\ H s + ^) + II^o|Ik ct([ _ 1)0 ]) + H v ^ll x -i ([ _ lj0]) } . 

for some non- decreasing function T : R + — > R + depending only on a. 

To prove Proposition 3.22 we shall proceed by induction on the regularity a. 

Definition 3.23. Given a such that — 1/2 < <r < s— 1/2, we say that the property H. a 

is satisfied if for any interval I <e (—1,0], 

(3.48) 

\\Vx, z v\\x°(l) ^ H\\v\\ H s+l) + H*b||y»([_i,o]) + H V ^ u llx-3([-l,0])} ' 

for some non- decreasing function J- ': R + — > R + depending only on I and a. 

With this definition, note that Assumption (3.47) means that property H—1/2 is sat- 
isfied. Consequently, Proposition 3.22 is an immediate consequence of the following 
proposition which will be proved in Sections 3.3 and 3.3.1. 

PROPOSITION 3.24. Let s > | + |. For any e such that 

1 1 d 

(3.49) 0<e ^2' £<S "2"2' 

if Ha is satisfied for some —1/2 < a < s — 1/2 — e, then H a +e * s satisfied. 

3.3. Nonlinear estimates. Let us fix e satisfying (3.49), a such that 

1 1 

— < a < s e 

2 ~ ~ 2 

and assume that 7i a is satisfied. We begin by estimating the coefficients a, (3, 7 
in (3.18) in terms of [|^[| ,1 , 

Lemma 3.25. Let J = [-1, 0] and s > \ + | . We have 



(3.50) 



a 

16 



2 1 2 ■ 

h 2 



_ L .+m\ X s- i(J) + \h\\ X s- hj <H\\ri\i^ 



X 2 (J) 

Proof. According to (3.11), (3.12) we can write 
h 2 

m z p) 2 = —+G with IIGII , K<F(\\n\\ .1 ^ 
Noticing that 1+ |y p p = 1 - 1+jvpp we obtain 



16 v 16 ; l + |Vp| 2 l + |V/o| 2 

and we use Lemma 3.20 with a = s — \ together with (3.11), (3.12). The estimates 
for j3 and 7 are proved along the same lines. □ 
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Lemma 3.26. There exists a constant K such that for all I C [—1,0], 
(3-51) \\Fi\\y-+-{i) < K ll7ll xS -| (J) \\d z v\\ Xa[I) , 

where F\ = jd z v. 

Proof. We shall prove that, on the one hand, if— 1/2 < a < s — 1 — e then 

(3.52) \\ldzv\\ L i {I . H *+ e) < \\l\\ L 2 {I . HB -^ ||3^|| i2(/ ^ CT+ i } , 

and on the other hand, if — e < cr < s — \ — e then 

(3.53) \^ dzV \\ L l{I;H°-h+c) ~ IMI-L^fp- 1 ) PzVW^^.h,) ■ 

Since s > e + ~ + d/2, if -1/2 < a < s - 1 - s then 

11 Id 
s-l + cr+->0, a + e <a + -, cr + e < s - 1, a + e < s - 1 + a + - - -. 

and hence the product rule in Sobolev spaces (2.16) implies that 

\\l{z)d s v(z)\\ Ha+e < h(z)\\ He -i \\d z v(z)\\ H „ + i . 

Integrating in z and using the Cauchy-Schwarz inequality, we obtain (3.52). On the 
other hand, if — e < cx < s — | — e then one easily checks that 

11 Id 

s-1 + ct>0, <r - -+e <a, o--- + e<s-l, a - - + e < s - 1 + a — -, 

and hence the product rule (2.16) implies that 

\\~f(z)d z v(z)\\ Ha _i +e < h(z)\\ Hs -i \\d z v(z)\\ Ha . 
Taking the L 2 -norm in z, we obtain (3.53). □ 

Our next step is to replace the multiplication by a (resp. /3) by the paramultiplication 
by T a (resp. Tp). 

Lemma 3.27. There exists a constant K such that for all I C [—1,0], v satisfies the 
paradifferential equation 

(3.54) d 2 z v + T a Av + Tp ■ Vd z v = F + F 1 +F 2 , 
for some remainder 

(3.55) F 2 = (T a - a)Av + (Tp - 0) ■ Vd z v 

satisfying 
(3.56) 

h 2 



F 2 \\y«+c { i) <K\ 1 + 



a 

16 



Proof. According to Proposition 2.11, we have 

\\au T' a u||j^ 7 ^ ||g|| jj r ll^ll/fM ) 
provided that r, fj,, 7 G R satisfy 

d 

(3.57) r + /x>0, 7 < r and 7 < r + /i — — . 

Since s > e + 1/2 + d/2, if -1/2 < a < s - \ - e then 

1 1 d 

s + a-->0, a + e<s, a + e < s + a - -- -, 
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and hence (3.57) applies with 



7 = 0" + e, r = s, \i = a 



This implies that if— 1/2 < a < s — h — e then 



(3.58) 



\\(T a - a)Av\\ L1{I . H<7+e) < (1 + 



a 



16 



< 



\\(Tp - (3)Vd z v\\ L1(I . HtT+e) - || . Hi , 
which yields 

< \\F2\\ L <{l + 

This concludes the proof. 



L 2 (I;H S ) 



a 



16 



X s -? (I) 



+ \\PW x s-h {I) }\\v*,M\x« { i) 



□ 



Our next task is to perform a decoupling into a forward and a backward parabolic 
evolution equations. Recall that by assumption 77 G H s+1 2 (R d ) with s > e+l/2+d/2. 
In particular, 77 G Cl +£ (R d ). 

Lemma 3.28. There exist two symbols a, A in r^(R d x [—1,0]) and a remainder F 3 
such that, 



(3.59) 
with 
(3.60) 
and 



(d z - T a )(d z - T A )v = F + F 1 + F 2 + F 3 , 
Ml(a) + Ml(A)<T{\\r,\\ HS+ t), 



\\F 3 \\ L2(rHa _ i+E) <F{\\ V \\ HS+i ) IIV^H^^^^, 
for some non- decreasing function T : R + — > R + . 

Proof. We seek a, A satisfying 
a(z; x, i)A{z; x, £) = —a(x, z) |£| 2 , a(z; x, £) + A(z; x, £) = —i/3(x, z) ■ £. 
We thus set 



1 



1 



(3.61) a = -H/H- v 4 « ler - (0 • o 2 ), ^ = 2 H/K+v 4 « Id - (/3 • 6 2 ) 

Directly from the definition of a, (3 (3.18), note that 



3c>0; V4a|£| 2 -03-£) 2 >c|£|. 

According to (3.50) the symbols a, A belong to T^R ' x [—1, 0]) and they satisfy the 
bound (3.60). Therefore, we have 

(3.62) (d z -T a )(d z -T A )v = d 2 z v-TpVd z v + T a Av + F 3 , F 3 = R v + R lV , 
where 

Ro(z) ■= F a (z)TA(z) - T a A, R\(z) := -T dzA( ^ z y 
According to Theorem 2.7, applied with p = e, Rq{z) is of order 2 — e, uniformly in 
z £ [—1,0]. On the other hand, since 

d t peL°°((-l,0);H B -h), d 2 z pZL°°((-l,Q);H s -i), 
according to (2.16) we have 

d z a,d z (3 e L°°((-l,0);F s -§ ) C L°°((-1,0); CJ" 1 ). 
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Therefore d z A G r^_ x (R d x [-1,0]). As a consequence, using Proposition 2.13, we 
get that R\{z) is also of order 2 — e. We end up with 

< -nihil 

*e[-i,o] H 

Now we notice that, given any symbol p and any function it, by definition of parad- 
ifferential operators we have T p u = T p (l — *$>(D x ))u for any Fourier multiplier (/ — 
ty(D x )) such that = for |£| > 1/2. This means that we can replace ||v(^)|| CT+ 3 
by ||Vu(z)|| CT+ i . We thus obtain the desired result from Lemma 3.27. □ 

3.3.1. Proof of Proposition 3.24- We shall apply Proposition 2.18 twice. At first 
we apply it to the forward parabolic evolution equation d z u — T a u = F (by def- 
inition Re(— a) > c|£|). This requires an initial data on z = —1 that might be 
chosen to be by using a cut-off function, up to shrinking the interval /. Next we 
apply it to the backward parabolic evolution equation d z u — Tau = F (by defini- 
tion He A > c|£|). This requires an initial data on z = (which is given by our 
assumption on /) and this requires also an estimate for the remainder term F which 
is given by means of the first step. 

Suppose that 7-L a is satisfied. Let Iq = [Co>0] such that 

||V^|| x<J(/o) < H\\v\\ H s+i ){ II/IU-+1 + \\Fo\\y°(i ) + ll V ^ll x -| ( [_ li0 ]) }■ 

We shall prove that, for any (i > Co, 
(3.63) 

l|Vx,*«||^+« ([Clj o]) < H\\v\\ H s+i){ + W f o\\y°+hio) + H V ^llx-2([-i j0 ]) }■ 

Introduce a cutoff function x such that 

x(Co) = 0, x(z) = l forz>Ci. 
Set w := x( z )(d z — Ta)v. It follows from (3.59) for v that 

d z w - T a w = F', 

where 

F' = x {z){F + F l + F 2 + F 3 ) + x'(z)(d z - T A )v. 

We have already estimated F\,F2,F$ and Fq is given. We now turn to an estimate 
for {d z — Ta)v. According to (2.4) and (3.60), we have 

ll^'Wo;^) - nM H^0W Vv hm ;H^) ^ ^W^) llv^|| XCT(/o) , 

and similarly 

\\ t av\\l^(i ;H-) < F(.M\ H s+$) \\Vv\\ Loo{Io . Ha) < F(\\t]\\ hS+ i) \\V x>z v\\ Xa(Io) , 
Consequently 

\\{d z -T A )v\\ Xa(Io) < F(\\r]\\ HS+ i)\\V XtZ v\\ Xa{Io) . 

This implies that 

(3.64) IMIx-(j ) < H\\v\\ H s + ll v ^llx-(/ ) ' 

(3.65) H^'Hy^+e^) < F(\\v\\ He+ i) \\V*,zv\\x°(Io) + W F o\\y-+^i ) • 

Since w(x,zq) = and since a 6Tj satisfies Re(— a(x,£)) > c|£|, by using Proposi- 
tion 2.18 applied with J = Iq, p = e and r = a + e, we have 

IMIx^/o) ^ H\\v\\ HS+ {||^||y. +e(/o ) + IIHIjr.+.-i (/o) } > 
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and hence, using (3.64) and (3.65) 

(3.66) \\w\\x<*+e(I ) < H\\v\\ H s+$ ){ l|V^|| XCT(/o) + \\F \\ Y a+e iIo) }■ 

Now notice that on I\ := [Ci, 0] we have x = 1 so that 

d z v — Tav = w for z G I\. 

Therefore the function v defined by v(x, z) = v(x, —z) satisfies 

d z v + T~v = -w for z e h = [0, —&], 

with obvious notations for w and A. By using Proposition 2.18 with J = I\, noticing 
that v\ z= q = v\ z= o = f, we obtain that 

\\v\\ Xa+1+£{Ti) < TihW^+O {\\f\\H«+l+e + My.+l+e^ + IH^-J (fi) } ■ 

Using the obvious estimate 

l|w|| yCT+1+e(?i) = iH| yCT+ i +e(/i) < Ikll^^^+j-Hsj ^ Ikllx^ai) ' 

it follows from (3.66) that 

(3.67) \\v\\ xt7+1+£{h) 

< ^(\\v\\ HS +l){\\f\\ H -+^ + ll V ^llx-(/ ) + H F ollw(/o) + IHI^-^))- 

We easily estimate d z v directly from d z v = Tav + w (by using (3.66) and the 
fact that Ta is an operator of order 1). This completes the proof of (3.63). This 
proves that if T-L a is satisfied then 7i a + £ is satisfied and hence concludes the proof 
of Proposition 3.24 (and hence the proof of Proposition 3.22). 

3.3.2. Proof of Theorem 3.18. Let v be the solution of (3.17) with data v \ z =o = 
/. By definition of the Dirichlet-Neumann operator we have 



1 + |Vp| 



2 



z=0' 



(3.68) 0(7]) f = ' n d z v -Vp-Vv 

o z p 

Now, by applying Proposition 3.22 with Fq = and Remark 3.21, we find that if v 
solves (3.17), then for any / d (—1,0], 

(3-69) l|V^|| XCT _ 1(J) < H\H HS+ II/IIh* • 

According to (3.11) and (2.17), we obtain that 
l + |Vp| 2 



-d z v — Vp ■ Vv 



<n\\v\ 

C0([z ,0];H"-^) 



H s+^) \\J \\H° 



d z p 

As a result, taking the trace onz = immediately implies the desired result (3.45). 

3.3.3. Proof of Proposition 3.19. Let 1/2 < gq < s. It follows from (3.66) applied 
with a = o"o — 1 and Fq = that 

\\x{z){d z v -T A v)\\ X a Q -i+e {h) < ^dMI^ i){ \\^x,zv\\ xao -i {h) }, 

for some cut-off function \ such that x(0) = 1. By using Proposition 3.22, we thus 
obtain 

(3.70) \\d z v - 2>|*=o ||^o-i+* < ^(IMI^+i) ll/IUo • 

The previous estimate allows us to express the "normal" derivative d z v in terms of 
the tangential derivatives. Which is the main step to paralinearize the Dirichlet- 
Neumann operator. 
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Now, as mentioned above, by definition of v, 



G(v)f = 1+ J VP|2 ^ -Vp-Vv\ z _ . 

d z p z - u 



Set 

1 + |VH 2 



According to (3.11 
(3.71) 
Let 



Ci := ' C2 := Vp. 



c S a-^*, + ll6 W^*)- JP(Nl ^* ) - 



# = Ci^ - Ca • V« - (T Cl 0,u - T C2 Vt;). 

1 _ d 

2 2' 



Since e < i and e<s— i — |,we verify that Proposition 2.11 applies with 



1 

7 = er -l + e, r = s--, /u = cr - 1, 
which, according to (3.71) and (3.69), implies 

11^ llc°(J;ii" CT 0- 1 + e ) — ^ r (ll 7 7ll i? s+^) ll/ll_f/ CT o • 

Furthermore, according to (3.70) and (3.71), we obtain 

T a d z v - T (2 Vv I 2=0 - (T Cl T A v - T i(2< v | 2=0 ) = R" , 

with 

\\R"\\ mo - i+e <H\\v\\ H s^) \\f\\ H . . 
Finally, thanks to (2.5), (3.71) and (3.60), we have 

\\ T (i(z) T A(z) ~ T Cl ( z )A(z)\\ H<TQ ^ H <r -i < \\Ci{z)\\ w ^ M\{A) < J"(||»7|| HS+ i), 
and hence 

G{rf)f = T ClA v - T <2< v I z=Q + R(rj)f 

where 

\\R(ri)f\\ H « -i + e<T(\\v\\ HS+ i)\\f\\ H « - 

Let 

a = 1+ £f A - ivp ■ e 1 z=0 = Va + iv^*)! 2 )^ 2 -^)^) 2 - 

Then 

G(r l )f = T x f + R(r ) )f, 
which concludes the proof of Proposition 3.19. 

3.3.4. Tame estimates for the Dirichlet- Neumann operator. This section is de- 
voted to the proof of Theorem 3.17. We thus consider the problem 

(3.72) d 2 z v + aAv + /3 • Vd z v - jd z v = 0, v\ z=0 = f, 

where / = f(x) is a given function. In the sequel we fix indexes 5, s, r G R such that 

,„ 1 d ,1 1,- ./I Id. 

(3.73) 0<5<-, s>l + --S, r>l,-<e=~-6<min(-,s-~-~). 

Let zo G (—1,0). It follows from Proposition 3.22 and the Sobolev embedding that 
we have 

ll v ^llco([, ,o];cr 1 - d/2 ) - ^OMI^) H/llfl" ' 

for some non-decreasing function J 7 . However, since we only assume that s > 3/4 + 
d/2, this is not enough to control the L°° norm of V X)Z v. The purpose of the next 
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result is to prove such bounds under the additional assumption that / belongs to C£ 
for some r > 1. 

Proposition 3.29. Let r > 1 and s > 3/4 + d/2. For any z < z\ < 0, we have 

l|V*,*v|lc*nL»([*i,o]xR<«) ^ HWvW^+O { H/llfl. + WfWc:} > 
/or some non- decreasing function T : R + — >■ R + . 
Remark 3.30. Since u|z=o = / € L^R^), we have also 

IMIl,°°(R<*x[«i,0]) - ll/lli°°(Pt d ) + l z ll ll3« w II.L°°(Rdx[*i,0]) 

<^(NI^ + i){ll/ll^ + ll/lk}- 

Proof. We have already proved that there exists zq S [—1, —1/2] such that 
(3-74) < ^(Ikll^) II/IIh. , 

and using again that t> | 2= o= /, 

(3-75) M\x^ H[zom <H\H\ H s^)\\f\\H^ 

As above, introduce a cutoff function x such that 

X(zo) = 0, x(z) = 1 for z > zi := -(z Q - -), 

and set w := x( z )(d z — Ta) v - Since it is convenient to work with forward evolution 
equation, define the function v by v(x, z) = v(x, —z), so that 

d z v + Trv = —w for z G I\ := [0, — z\\. 

We split v &sv = v\ + V2 where v\ is the solution to the system 

8 z v 1 + T x v 1 =0 for *e/i, 5 1 (0)=5;(0) = / 

given by Proposition 2.18, while ^2 = w — vi satisfies 

d z v 2 + TjV2 = —w for z e h, ^(0) = 0. 

According to (3.66) with a = s — 1 and e defined in (3.73), we obtain, using (3.74) 

IWjr-^-Vo) - ^(II^Uhh-j) l|V*,*w||^-i (/o) < ^(IN^ i) II/IIh* • 
which in turn implies, according to Proposition 2.18, 

(3.76) \\v2\\ Xs+ e {Il) < H\\ri\\ H s + i)(\\V x ,zv\\x°-Hi ) + INI^Ji))- 

Using Proposition 2.19 with r = 1, Tq = —1, we obtain (with / = [zo>0]) 

ll«l|lcO(j;CI) - K X (H/IIC5 + H^llco^C- 1 ))- 
But, according to (3.74), (3.75) and (3.76) 

IHIcH/iOT 1 ) ^ C (\Hc°(I;Hs-i) + INI C o ( j.tf S -l)) < -^(Ihll^ 
which since d z v\ = T-rui , implies also 

<^(IMI^ + i)(ll/ll^ + ll/lk). 

On the other hand, according to (3.66) with a = s — 1 and e defined in (3.73), we 
obtain, using (3.74) 

W w h*-^- s (i ) - -^IMI^) \\v*,zv\\x-m ) ^ FihWjjs+O \\f\\ H s ■ 
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Let us go back to the study of v%. According to (3.76), and since d z V2 = —T^V2 — w 
we get 

(3-77) ||V :Ei2 W2|| Xs -i+ e( 7 l) < F{h\\ HS+ i) ||V a .,*u|| x ._i (Jo) • 

Set 

d 

m = s — I + e — — . 
By the Sobolev embedding we have 

|| ^X,Z V 2\\]^!X)rJ 1 .Qm.'\ ^5 || ^ x,zV2\\xs-i+e(J 1 ) • 

Using (3.77) and our previously established estimate for || Vj^uH^s-wj s, this yields 

\\^x,zV2\\ Loo{ i i;C ^) ^ ^(Nlgf*) ll/llfl» ■ 

By assumption s > 1 + d/2 — 5 with 5 < 1/4, so that 

d 1 P d 1 

m = s-l + e-- = s-l + --6-~>--28>0, 

and hence 

Wx,zV2\\ Lx(R d x[0 _ Zl]) < ^(|kH HS+ i ) \\f\\ H s • 

This completes the proof of Proposition 3.29. □ 
3.3.5. Proof of Theorem 3.17. As already mentioned, by definition of v, we have 



l + |Vp| 



2 



G(v)f = ' r ' d z v -Vp-Vv 

OzP 



z=0' 



To perform the paralinearization of the Dirichlet-Neumann operator, we are going 
to revisit Section 3.3.3, using tame estimate at each step. Let us begin by recalling 
our notations. We have 

(d z - T a )(8 z - T A )v = F 1 + F 2 + F 3 , 

where a, A are given by (3.61) and Fx, i*2, F3 by (3.51), (3.55) and (3.62) Recall that 
according to Propositions 3.22 and 3.29, for any zq G (—1,0), we have 

^ \\ Vx ' zV \\c«([z ,o]-,Hs-i)nL2((z ,oy,H°-?) ~ ^vW^+i) \\f\\ H s . 

II V- 



! »' i/:-i.R' ? <i :u.o)) - • ?r (ll 7 ?ll ff .+i){ll/llH» + ll/llcr} 



The end of the proof of Theorem 3.17 is in four steps. 
Step 1. Tame estimates in Zygmund spaces. 

Lemma 3.31. The coefficients a, (3, 7 defined in (3.18) satisfy 

(3.79) ll( a >£)llcO([_l,0];Ci /2 ) + HtII C -0 ([ _ 1>0 ];C7- 1/2 ) - - 7r (ll ? ?ll^+l/2)(ll??ll C .3/2 + 1). 

PROOF. It follows from (3.10) that 

(3-80) \\^ x P\\c°([-l,0];Cl /2 ) ~ H^llc, 372 ' 

Moreover using (3.7) we get 



(3.81) 



dzP-j 
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On the other hand recall that we have, according to (3.10) 
(3.82) l|V^p|| CO([ _ 1)0];HS _i (Rd)) < K h\\ HS+ , 

=> ll v ^H c o (hli0]; ^i (Rd)) < *(i + M H s +h ), 

and consequently by Sobolev embeddings, 

l|V^ 2 p|| LOO(( _ 1>0)xRd) < K(l + IMi^+i). 

Using H s+l / 2 C C* +1 ^ 2 rf ^ 2 C C*^ 4 , we deduce the estimates for a — 1 and /3 from 
the tame product estimate in Zygmund spaces (see (2.18), (2.22)). The estimate 
for 7 follows from (2.20), (2.22). □ 

Step 2. Estimates for the source terms. 
Lemma 3.32. For any zq G (—1,0), and any j = 1;2;3 we have, 

(3-83) \\Fj 11^,0)^-1) < HWlWg+i' H/M I 1 + II^Hcl + H/Hor} • 

Proof. By using the product rule (see (2.19)) 

llwi^H^- 1 - ll^lllc.-i/a H^ll^s-i + \\u2\\l°° 
we find that F\ = jd z v satisfies 

+ ^c°([z ,o];C- 1/2 ) W dzV hmz ,oy,H s ~h) ■ 

The desired estimate for F\ follows from Lemma (3.31), (3.50) and (3.78). Let us 
now study 

F 2 = (T a - a)Av + (Tp - 0) ■ Vd z v = - (T Av a + R(a, Av) + T vdzV ■ /? + R(fi, V8 z v)) . 
According to (2.12), we obtain 

H^^Q!^)!!^-! < ||Au(z)|| c -i \\a(z)\\ Hs , 

\\TvdM*) -P^Whs-i % \\Vd z v(z)\\c^ HWIIfl.. 
On the other hand, since s — 1 > we can apply (2.11) to obtain 
\\R(a,Av)(z)\\ Hs -i < ||At;(z)|| c -i \\a(z)\\ H . , 

WR^Vd.vMW^ < \\Vd z v(z)\\ c -x Wz)\\ H .. 
Consequently we have proved 

(3.84) \\F 2 || £ 2([ 20i0 ]. HS -i) < \\ A v\\ CO ( [z0)0] . c -^ ||a|| L 2 ([20i0 ].^ s) 

+ W^dzvWcodzoflhc, 1 ) WP\\l2([z ,o];H*) ■ 

Notice that 

ll^llc," 1 ~ W Vv Wc° ~ ll Vw lli,°° » \\^9 z v\\ c -i < \\d z v\\ c o < \\d g v\\ Lao 

and consequently, according to (3.78) and Lemma 3.25, we conclude the proof of the 
claim (3.83) for j = 2. It remains to estimate F3. In light of (3.78) it is enough to 
prove that 

(3.85) ||F 3 || 
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for some non-decreasing function. Directly from the definition of a and /?, by using 
the tame estimates in Holder spaces (2.18), we verify that the symbols a, A belong 
to rJy 2 (R d x I) and that they satisfy 

(3-86) M\ /2 (a)+Ml /2 (A)<T{\\rj\\ HS+ i)\\rj\\ c!/2 . 

Notice for alter purpose that we use here only s > \ + | . Moreover, according 
to (3.86) we have 

MU{d z A)<F{H\ HS+h ), 

and by definition (see the proof of Lemma 3.28), 

R (z) := T a{z) T A(z) - T a A, R x {z) := -T 9zA . 

As in the proof of 3.28, we deduce, using Theorem 2.7, (ii) applied with p = 1/2, 
and Proposition 2.13 with p = e — 1, that 

(3.87) ^sup Q] \\R (z)\\ h ^^ + \\Rl(z)\\ H ^ Hfi < HM H s +h )\\v\\ §, 
which implies (3.85). This completes the proof of Lemma 3.32. □ 

Step 3 : elliptic estimates. Introduce a cutoff function k = k(z),z G [—1,0] such 
that k{z) = 1 near z = and such that n{z\) = (recall that I\ = [zi,0]). Set 

(3.88) W := k(z)(8 z -T A )v. 

Now it follows from the paradifferential equation (3.59) for v that 

8 Z W - T a W = F', 

where 

F' = k{z){F x +F 2 + F 3 ) + K'{z){d z - T A )v. 
Our goal is to prove that 

( 3 - 89 ) » w W^l) ^ F ^\\ HS+h , \\f\\ H .) {l + H\ ci + ll/lbr} • 

We have already proved that 

ll^'L^-!) < HWv\\ HS+h , \\f\\ H s) {l + h\\ c S + \\fWcr} • 

We now turn to an estimate for (d z — T A )v. To do that we estimate separately d z v 
and T A v. Clearly, by definition of the space X s-1 , noting that /i C 4 we have 

On the other hand, as in the previous step, since A4q(A) < C(||r/|| we have 

<-^(ll»7ll H H-i)l|Vx,*w||^.-i ( j o) . 
Now recall from Proposition 3.22 that 

||v x .^|| xs _ 1(/o) < TiWvWjjs+i) \\f\\ H ' ■ 

Therefore, 

H^W^-*) - ll^-^ll^^i, <^(Nl^i)ll/llfl.. 
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We are now in position to estimate W = n(z)(d z — T A )v as well as to estimate the 
last term in the definition of F'. We end up with 



IF'I 



aO{ 1 + IMI c i + 



C: 



}■ 



\W\\ i 



Since W(a;, 21) = (by definition of the cutoff function k) and since a S sat- 
isfies Re(— a(x,£)) > c|£|, by using Proposition 2.18 applied with J = I\, p = e 
and r = s — 1/2, we have 

r0 ^^ll^ll^){ll^ll^ (/l) + ll^ll i2(7i; W)}- 
By definitions, we have 

IIWII , . ^ < IIWII . 1/ , and \\F'\\ vs -ut ^<\\F'\\T2n us -n- 

We thus conclude that W satisfies the desired estimate (3.89). 

Step 4 : paralinearization of the Dirichlet-Neumann. We conclude the proof by 
mimicking the analysis in §3.3.3. Again, we shall only use the following obvious 
consequence of (3.89): ||W| z= o|| s _i is estimated by the right-hand side of (3.89). 



Since W\ 



(3.90) 



2=0 



d z v - T A v\ 



2 = 



we thus have proved that 



\d z v - T A v\ 



2 = I 



Now, recall from S3. 3.3 that we have 



<H\\v\ 



H'){l+\\V\ 



3 + 



c;} ■ 



G(v)f = 1 + ' VP|2 ^ -Vp-Vv 



and 



Ci := 



d z p 

i + |vh 2 



z=0' 



d z p 



(2 := Vp. 



As for the coefficients a, /3, we have 
(3.91) 



Ci- 


4 


77 


Ci- 


4 
77 



+ HC 2 || ioo(hl0] c y 2) < ^(l|r?||^ + i)NI| c 3/ 2) 

L°°([-l,0];Cl' ) 

+ l|C2|| i oo(r_ 1 ,01^-1/2) < F(\\n\\ HS+ i)- 



(3.92) 
Write 

with 

# = r^Ci " 2V„ • C 2 + R((i,d z v) - R{( 2 , Vu) 
According to Proposition 3.29 and (3.92), we obtain 



ClO** - C 2 • = T fl 0,t> - T C2 X7v + 



||T 9z ,Ci - T Vv • C 2 || ioo(Ji . ifS -i ) < ^(IMI^+i , 



^{i + NII^I + H/llcr} 



To estimate the last two terms, we use the rule (2.11), (3.91), the Sobolev estimates 

_3_ 

2ft > 



on d XtZ v proved in Proposition 3.22, and the remark that R((\, d z v) = R((,i — jr,d z v) 



(because constants count only for low frequencies). By so doing (and using that s > 
\ + |), we find that 



7/ 3 + 

C? 
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Furthermore, (3.90) implies that 

T (l d z v - T (2 Vv | 2=Q = T Cl T A v - T iC2< v \ z=Q + R" , 
where \\R"\\ „ 1 satisfies the same estimate as \\R'\\ , „ 1, does. 

Thanks to (2.5) we have 

\\T Ci T a -T CiA \\ hs ^ hS _i < WCiW^Ml^ + WdW^MKA) 

< -^(11^11^+1)11^11^3/2, 

and hence 

G(v)f = T^av - T <2< v | z=Q + R(ri)f 

where 

\\R(v)f\\ H s-i < H\\v\\ HS+i )\\v\\ ci /i ll/llfl.- 

Since G(rj)f = T\f + R(rj)f, this concludes the proof of Theorem 3.17. 

4. A priori estimates in Sobolev spaces 

Recall that the system reads 

( d t r) - G(n)i> = 0, 



(4.1) 



^ + OT + -|v^| = 0. 



As already mentioned, we work with the unknowns -B = B(jj,ijji) and V = V(r],ip) 
defined by 

Vr] ■ Vtp + G(77)V> 
I + IVt/P 

It follows from Theorem 3.18 that, for all s > 1/2 + d/2 and all (77, -0) G # s+ 5, 5 
and V are well defined and belong to H s ~z. Moreover, we shall prove that if they 
belong initially to H s then this regularity is propagated by the equation. We shall 
prove estimates in terms of 

M S (T) := sup • [K^flW.VCr),^))!!^^^^ , 

t£[0,TJ 

(4.2) Ms <° := IIW)^(°)^(°)^( ))ll^l xfl . x ^i > 
Z r (T) := [£{\\v(r)f c r+ i + ||(B(r), V(r))||^ xCr } dr 

The main result of this section is the following proposition. 

Proposition 4.1. Let d>l and consider s, r g]1,+oo[ such that 

3d 1 d 

(4.3) S> 4 + 2' s +4"2 >r>L 

Consider a fluid domain such that, there exists h > such that for all t G [0, T], 

(4.4) I (x, y) G R d x R : r/(t,a;) -h<y < r](t,x)\ C O(i). 
Assume that for any t G [0, T], 

a(i,x) > c , 
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for some given positive constant cq. Then, there exists a non- decreasing func- 
tion J 7 : R + —7- R + such that, for all T 6 (0,1] and all smooth solution (n,tp) of 
(4.1) defined on the time interval [0,T], there holds 

(4.5) M S (T) < T(T(M sfi ) + TT(M S (T) + Z r (T))). 

Remark 4.2. (i) If s > 1 + d/2 then we have Z r (T) < M S (T) for r = s - d/2. 
Then, the estimate (4.5)is an a priori estimate in Sobolev spaces, 
(ii) For s < 1 + d/2, to obtain a closed system of inequalities, one need an a priori 

estimate for Z r . This will be done in the upcoming paper [3]. 
(hi) The assumption (4.4) holds provided that it holds initially at time and 

|t=o ll^s+i < e, 

for some small enough positive constant e. 

4.1. Notations. From now on, we consider a hxed time < T < 1, indexes 
s,r satisfying (4.3) and a smooth solution (r),ip) £ C°°([0, T]; H°°(R d )) of (1.6). 

Given I C R, and two functions f : I —> X and g: I — > Y, we use the notation 
\\f\\x < Nly to say that \\f(t)\\ x < \\g(t)\\ Y for all t £ I. 

Hereafter, T always refer to a non-decreasing function J 7 : R + — > R + depending 
only on s, d and h, cq (being of course independent of T and the unknowns). 

4.2. A new formulation. Since we consider low regularity solutions, various 
cancellations have to be used. We found that these cancellations are most easily seen 
by working with the incompressible Euler equation directly, and hence we do not 
use the Zakharov formulation. This means that we begin with a new formulation of 
the water waves system which involves the following unknowns 

(4.6) C = V7/, B = dy4>\y =V , V = V X (f)\y = , n , O = -dyP\y =V , 

where recall that cj) is the velocity potential and the pressure P = P(t, x, y) is given 
by 

(4.7) -P = d t cf ) + ] 1 \V x ^\ 2 + gy. 

Proposition 4.3. We have 

(4.8) (St + V • V)B = a - g, 

(4.9) (d t + V -V)V + a( = 0, 

(4.10) (d t + V ■ V)C = G{rj)V + (G( V )B + 7 , 

where the remainder term 7 = 7(77, tp, V) satisfies the following estimate : if s > 5+5 
then 

(4-11) ||7|| . 1 < -F(||(r?,^,y)|| s+ i ). 

Remark 4.4. In the case T = 0, one can see that at least formally 7 = 0. 

Proof. For any function / = f(t,x,y), by using the chain rule, we check that, 
with V = Vje, 

(d t + V- V)(f\ y=v{t>x) ) = (d t + V ■ V)f(t, x, r,(t, x)) 

= [d t f + v<a • v/ + d y f(d tV + v ■ vr?)] I y=ri{tix) 

= [{dt + V x>y (f>-V Xty )f] \ y=v{tjX y 
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since dtr] + V ■ Vr/ = B (see (3.22)). Applying d y to (4.7), this identity yields (4.8). 
On the other hand, applying d Xk to (4.7), the previous identity gives 

(d t + v-v)v + (vp)\ y=v = o. 

Since P\ y =r) = 0, we have 

= V(P\ y=v ) = (VP)\ y=v + (d y P)\ y=v V V , 

which yields (4.9). 

To derive equation (4.10) on £ := V77 we start from 

d t r) = B -V • Vt? 

Differentiating with respect to a;, (for i = 1, . . . , d) we find that 5,r/ = d Xi rj satisfies 

d 

(4.12) (dt + V- V)d iV = - ^ diV 3 d,rj, 

i=i 

Starting from the definitions of B and V (B = d y 4>\ y=v , V = V(j)\ y=ri ), and using 
the chain rule, we compute that 

d d 
diB - diVjdjr] = [did y (p + d^d^] \ y=v ~Yl d jv[didj(f> + dirjdjd y 4>] \ y=rj 
i=i 3=1 

d d 

= [ d y d i<t> ~ Yl d jVdidj4>] \ y=r] + d^d^ - ^ djijdjdy 
3=1 ' 3=1 

Therefore 



y=v 



(4.13) (d t + V-V)d iV = [d y d i( f>-Vr]-Vd i( t>] \ y=v + d i r)[d y {d y <t>)-Vr ] -Vd y 
Let now 9i be the variational solution of the problem 

A Xty 0i = in Q, 0<| E = F if <9„A = on T. 

Then 



!/=»? 



G(r?)K = V^+WWq^ ^ = (d y e i -Vr ] -V8 l )\x. 
Then we write 

(4.14) (d y — V77 • V)di(j) = G{rf)Vi + B4, where Ri = (d y - Vr? • V)(d i( f> - ^)| s . 

Due to the presence of the bottom we have to localize the problem near E. 

Let xo S C°°(R), rn G H°°(R d ) be such that xoO*) = 1 if z > 0,xo(^) = if z < -\ 
and 

Set 

^(x, y) = Xo( y ~^ (x) ) - f )(x, y). 

We see easily that 

Moreover Ui satisfies the equation 

(4.15) A Xty Ui= [A^xo Cl^ ^id^-Oi) := F { 
and, with a slight change of notation, we have 

h h 

(4.16) suppFi cSn:= {(x,y) : x G R d ,r/(x) - - < y < rj(x) - -}. 

2'5 L 2 5 
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Moreover by ellipticity (see Lemma 3.16) we have for all a G N d+1 
(4.17) 



D x,y F i\\L™{Si i)nL 2 (Si i) < Call (V, B) \\ i ±. 

2 ' 5 1' 5 



Now we change variables. We set x = x,y = p(x,z) = (1 + z)e Sz( - Dx ^rj(x) — zr]*(x) 
and gi(x, z) = gi(x, p(x, z)). Since we have taken 5\\r]\\ . i < \ it is easy to see that 

H 2 ^ 

Yjj. Now, according to section 3.2.1, f7j is 
(<9 z p) 2 - 



on the image of Si i one has —h<z< 

2'5 

a solution of the problem 

(d 2 z +aA + (3-Vd z - jd z )U t 



l + |Vp| 



Due to the exponential smoothing and to (4.17), on the support of Fi the right hand 
side of the above equation belongs in fact to C z ((—h, 0); if°°(R d )). In particular we 
can apply Proposition 3.22 with / = 0. It follows that 

^ x ' z ^C0( [zo ,0);H^ m) ^ ^(II^H^+^dl^H^d-M) + ll V ^H x ^ a _ li0]) ) • 

Notice that according to the constructions of variational solutions and (3.21), the 
norm of U{ in X~z{\— 1,0]) is bounded by 



Since 



we deduce that 



R. 



l + |Vr/p 



l + d(D x ) V 



1 + \\Vi\\ i 



d z -vri-v)u t 



2 = 



-i<H\\v\ 



H i + \\m H 0<Hh\\ H s + ^\M\H°,\\v\\ H °), 



since s > | + |. We use exactly the same argument to show that 
(4.18) (d y -Vr j -V)dy<j> = G(ri)B + Ro, 

where Rq satisfies the same estimate as R{. This completes the proof of Proposi- 
tion 4.3. □ 

Following the same lines, we have the following relation between V and B. 
Proposition 4.5. Let s > \ + |. Then we have G(n)B = — divV + 7 where 
ll7ll W <^(ll(^F,S)|| Wxifix ^). 

PROOF. Recall that, by definition, B = d y cj)\ y=rj and V = V(f)\ y=ri . Let 9 be the 
variational solution to the problem 



A x , y 9 = 0, 
Then G{r,)B = (d y 9 - Vt] ■ V0) ' 



\y=v 



B, d n 6\ r = 0. 



y=v 



Now let 



(d y 9 - Vn ■ V0) 



y=v 



d v 6. We claim that 



divV. 



Indeed, on the one hand we have 



(d y 9 - Vr, ■ V0) = d 2 y <j) - Vr/ • Vd y <p, 

and on the other hand 



Ki<d 



l<i<d 
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y=n 



Then our claim follows from the fact that A Xt y(p = 0. Now we have 

A x>y (9-9) = 0, (9-9)\ y=v = 0, 
so, as in the proof of Proposition 4.3, we deduce from Proposition 3.22 that 

(d y - v v ■ v)(* - e)\\ HS _ h < HKv, v,b)\\^ xH ^), 

which is the desired result. □ 

4.3. Estimates for the Taylor coefficient. In this paragraph, we prove sev- 
eral estimates for the Taylor coefficient. 

Proposition 4.6. Let d>\ and consider s,r such that 

3 d 
S> 4 + 2' • 

For any < e < min(r — l,s—^ — &), there exists a non- decreasing function T such 
that, for all t £ [0,T], 

(4.19) lk(*)-ffll^^<^(ll(^^^^)(i)ll^ x ^ xffSxJ/s 
and 

(4.20) \\a(t)\\ ci+E + \\(d t a + V-Va)(t)\y 



<mM,V,B)(t)\\ He+ , xHa+hxHsxHs ) 1 + v,B m \\ cr+ ^ r 



We shall see below that the estimate of \\a — g\\ s _i follows directly from the elliptic 
estimates (see Proposition 3.22). The estimate of the C £ -norm of dta + V • Va is the 
easiest one. The main new difficulty here is to prove a tame estimate for ||a||^i +e . 
Indeed, there are several further complications which appear in the analysis in Holder 
spaces. To prove the above estimates on o, we use that 

a — dyP\y=rji 

where 

P = P(t, x,y) = - [d t 4> + i | V^| 2 + ^(d y 4>) 2 + gy) . 

The basic idea is that one should be able to easily estimate P since it satisfies an 
elliptic equation. Indeed, since A x>y (j) = 0, we have 

A x ,yP = ~ |V^| 2 . 

Moreover, by assumption we have P = on the free surface {y = n(t,x)}. Yet, this 
requires some preparation because, as we shall see, the regularity of P is not given 
by the right-hand side in the elliptic equation above. Instead the regularity of P is 
limited by the regularity of the domain (i.e. the regularity of the function rf). 

Hereafter, since the time variable is fixed, we shall skip it. We use the change of 
variables (x,z) h-» (x,p(x,z)) introduced in §3.2.1. Introduce ip and p given by 

<p(x, z) = (j)(x, p(x, z)), p(x, z) = P(x, p(x, z)) + gp(x, z), 

and notice that 

a- 9 = ~-^—d z p \ z=0 ■ 

OzP 
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The first elementary step is to compute the equation satisfied by the new unknown v 
in {z < 0} as well as the boundary conditions on {z = 0}. Set (see (3.14)) 



We find that 



A = (Ai,A 2 ), Ai = 
(Af + A$)<p = in 



OzP 



A, 



V-^dz 

OzP 



AV m 



1< z < 0, 

-Kz <0, 



(Af + A|)p 
(Af + Aj|)p = in z < 0, 

together with the boundary conditions 

P = 9V, Aip = #-a on z = 0, 
A2</3 = V, An/9 = B, on z = 0. 

According to (3.5) and Remark 3.21, we have the a priori estimate 

\\^x,z<f\\ 



'X-2(-|,o) 

while according to Proposition 3.22 
(4.21) 



< Hh\\ HS+ i 



i . 

H2 



iv^pll 



< J" \\K\ 



l X-1(-l,0) — V" "XI (-1,0) 111 rl "X3(-l,0) 

where lZ(x,z) = R(x,p(x,z)) and i? is defined in Definition 1.5. 
Expanding A^ + A|, we thus find that p solves 

^ 4 22 ^ d 2 z p + aAp + (3 -Vdzp-^d z p = F (x,z) for z < 0, 

p = on z = 0, 

where a, /3, 7 are as above (see (3.18)) and where 

(4.23) F = -a |AV| 2 • 

Our first task is to estimate the source term Fq. 

Lemma 4.7. Let d>l and consider s,r e]l,+oo[ such that 



<H\\(v^)\\ H s +h )- 



3 d 

S> 4 + 2' 

Then there exists zq < such that 

\\Fo\\ 1 +||Fn|| 

11 0|l £i([* ,o];ir s -2) T 11 0|l L 2 ([20i0] . c ; 



1 d 
fl+5 -->r>l. 



PROOF. Since [Ai, A 2 ] = we have 

(Af + A|)A 2¥ > = 0, (A?+A|)A l¥ > = 0. 

Since A 2 yj| z= o = V and Anp\ z= o = B, it follows from Proposition 3.22 (and Theo- 
rem 3.8 which guarantees that ^7 x ,zf £ X~2(zq,Q)) that 

||V^A i Hlx.-i([^,o])^^(ll(^V' ) ^S)[|^ 4xWx ^ xH J. 
By using the easy estimate (3.11) 

ft 



Vrff 1 + 

1 xHl] C°([z ,0};H s -'2) 



d z p 



C°([zo,0);H 8 -?) 



^ \\v\ 
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and the product rule in Sobolev spaces, we obtain 
(4.24) ||A J A^|| c0(Ni0]; ^ 1) + UM\ L2MHa - h) 

<T(\\(r],iJ),V,B)\\ s+ i s+ x ). 

Since i? s— 2 is an algebra, according to Lemma 3.25, we obtain 

h 2 I 



< 4 - 25 ) 11^11^,0,^^1 + 



a 



16 



, IIAjA^I 



L2([ Z( „0];ff s -2) 

H s)- 



'H s+ ? xiT s+ 2 x_ff s x 



By interpolation, (4.24) also implies that 

l|A 7 -A fc </>]| An , . a < ^(IKn.^F.B)]! s+ i s ,i ). 

11 J Krll L 4 ([ 2o ,0];H s -3) - VH V /s ) JU H s+^ xH s+^ xHsxHs ; 

3 

Since s > 3/4 + d/2, the Sobolev space H s ~i is an algebra and hence, according to 
Lemma 3.25, 



a |Aj-Afcy>|" 



£»([*G,0];fl— I) 



< 1 + 



o 



16 



L°°([z ,0];H B -Z) 



1 J fc ^ ll L4([ Z0 ,0];H s -5) 



The Sobolev embedding i? s 4 c C* 4 2 then yields 
(4.26) 

This completes the proof. 



i?0| i2 (M]^-^ J ( |IW ' F ' 5)i ^x^xH,J' 



□ 



It follows from Lemma 4.7 and Proposition 3.22 applied with a = s — 1/2 that there 
exists zq such that 

(4.27) llV.^ll^l M) <^(ll(^^)ll^iJl^o||^ {[ad ^ 4) ). 
where we used the estimate (4.21) According to (4.25), this implies that 

(4.28) Wxzp\\ s i„ , < r(\\(T},il>,V,B)\\ s+ i s+ i ). 

which in turn implies that \\a — g\\ s _i is bounded by a constant depending only on 
11(77,^)11^+1 and IKV.B)!!^. 

We shall now deduce the estimate of HaH^i/a+e, which is the most delicate one, from 
the following result. 

Proposition 4.8. Let d > 1 and consider (s, r, r') G R 3 such that 

3d Id , _ 

S >- + -, s+--->r>r >1. 

TTien i/iere exists zq < swc/i i/iai 
(4.29) 



c,. 



llVxzpll r / i < J r (\\(v,il>, V,B)\\ B ,i s+ i Ml + ltoll 

/or some non- decreasing function T depending only on s,r,r' . 

Proof. It follows from (4.28) and the Sobolev embedding H s ~i (R d ) C C* 4 (R d ) 

II V*,*p|| LO o M;C i/4 ) < ^(|| (V, 1>, V, B) \\ rS+ i xhS+ i xHSxHS )- 
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that 
(4.30) 



To prove (4.29), let us revisit the proof of Theorem 3.17. With the notations of 
Section 3.3 (see (3.88)), W = n(z)(d z — T A )p satisfies a parabolic evolution equation 
of the form 

8 Z W - T a W = F + F 1 +F 2 + F 3 + F 4> 
where the symbols a and A are as defined in Lemma 3.28, Fq is given by (4.23) and 

Fi = jd z p, 

F 2 = - {T Ap a + R(a, Ap) + T vdzP ■ + R(J3, Vd z p)) , 

F 3 = R p + Rip = (T a T A - T a A)p + T dzA p, 

F 4 = K'{z){d z -T A )p. 

In light of (4.30) and Proposition 2.19 (applied with ro = 1/4, q = +oo), one can 

r'— -+S 

reduce the proof of (4.29) to proving that, for some 5 > 0, the L°°([zo,0]; C* 2 ) 
norm of W is bounded by the right-hand side of (4.29). Again, by using Proposi- 
tion 2.19 (with q = 2), the former estimate will resulted from the following lemma. 

Lemma 4.9. There exists zq < and e > such that, for i £ {0, . . . , 4}, 

To prove Lemma 4.9 we begin by recording two easy refinements of previous bounds. 
First of all, since r— 1 < s— 1/2— d/2 we have H s ~z C Cl~ l and hence it immediately 
follows from (4.28) that 
(4.31) 

\\ v xM\c°({ Z1 ,o];c:-') £ ll v ^llx^( K o]) ^(ll^' F '^ll/r+W+Wx/p)- 

The key point is that we now have an L°°-estimate for V IjZ p which does not depend 
on the higher order norms (compare with Proposition 3.29). 

Our second observation concerns the estimates in Holder spaces for the coeffi- 
cients a, /3, 7. Of course, since we now allow estimates to depend on \\r)\\ r+ i instead 

of \\rj\\ 3 , we have the following variant of Lemma 3.31: 
c ^ 

IMI r 1 +WPW r 1 +II7II r 3 <T{M\„ s + h)(h\\ r+1+1)- 

However, for later purpose, we need estimates for a(z),(3(z) in C£ instead of C* 2 

-i r — — r' — — 

and 7(2;) in C£ instead of C* 2 . Similarly, it is not enough to control 7(2;) in C* 2 
since r — I might be negative. We have such estimates to the prize of a lack of 
uniformity in z. Namely, we have the following estimate: 
(4.32) 

II«IIl 2 ([, ,o];c;) + II/^IIl^co];^) + bll^^cr 1 ) - ^(IMI^+iXIMI^+i + 1)- 

Again, such estimates follow from the definition of p (by means of the Poisson kernel), 
the Sobolev embedding H s+ 2 c W 1 ' 00 and tame estimates in Holder spaces (2.22). 
We are now ready to conclude the proof of Lemma 4.9 

Estimate of Fq. Similarly, since 

c s-3/4-d/2 c C r-l^ ^ 2 g^ implieg that 

ll^oll^ohcr 1 ) ^ T ^^v,B)\\ H3+hxHS+hxHsxH5 ). 

Estimate of F\ . The term F\ = jd z p is estimated by means of 

ll^lllLK^o.O];^- 1 ) - K W d ^h^([z o ,0];Cr 1 ) 1 1 ^ 1 1 i| ([^o ,0] x ) • 

The desired estimate then follows from (4.31) and (4.32). 
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Estimate of F<i- According to (2.13) with m = l,s = r, we obtain 

|| r Ap(*)a(*)||oj-i % \\ A P( z )\\c^ Haulier » 
||2v8. pW • /^Wllcj-i < llV^p^H^x ||/3(z)|| cr . 

On the other hand, since r > 1 we can apply (2.10) to obtain 

||#( a) Ap)(2)|| crl < \\Ap(z)\\ c -i\\a(z)\\ C r, 
\\R(J3, Vd z p){z)\\ c r-, < ||V0*p(*)|| c -i \\P{z)\\ Cl . 
By using (4.31) and (4.32), we conclude the proof of the claim in Lemma 4.9 for i = 2. 
Estimate of F%. Using (3.87) with fj, = s— 1/2 we find 

IIWIIcj-i < II^WII^-i <^(hll HS+ i)ll^ll c 3 ||v^ P (z)|| Hs . 

and hence 

ll i? 3ll L2([20 , 0]; cr 1 ) ^ ■PflMflf i)IM c § 11^.^11^-4(^,0]) ' 

by definition of X s_1//2 ([zo, 0]). The desired estimate follows from (4.28). 
Estimate of F4. This follows from (4.31). 

This completes the proof of Lemma 4.9 and hence the proof of Proposition 4.8. □ 

4.4. Paralinearization of the system. Introduce 
(4.33) U = V + T ( B. 

To clarify notations, let us mention that the ith component (i = l,...,d) of this 
vector valued unknown satisfies XJ% = Vi + Tg iT1 B. The new unknown U is related to 
what is called the good-unknown of Alinhac in [5, 1, 6, 8]. 

To estimate (U, Q in Sobolev spaces, we want to estimate ((D x ) s U, (D x ) s ~ 2 Q 
in L°°([0,T];L 2 x L 2 ) where (D x ) := (/ — A) 1 / 2 . However, for technical reasons, 
instead of working with ({D x ) s U, {D x ) s ~2 £) ; it is more convenient to work with 

U s := (D X ) S V + T ( (D X ) S B, 

(4.34) 

( S :=(D X ) S (. 

Proposition 4.10. Under the assumptions of Proposition 4-U there exists a non 
decreasing function T such that 

(4.35) (d t + T v -V)U s + T a Cs = fi, 

(4.36) (d t + T v -V)( s = T x U s + f 2 , 

where recall that A is the symbol 



\(t;x,0 := ^(l + lVrKt,*)! 2 )^ 2 -^^)^) 2 , 
and where, for each time t 6 [0, T], 

(4.37) \\(h(t),Ht))\\ L2xH ^ 

< HUt)\\ HS+h , \\(v,B)(t)\\ HS ) |i + Ut)\\ c:+ i + ll(v,B)(0[lcr} • 
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Proof. The proof is based on our main result about the paralinearization of 
the Dirichlet-Neumann operator (see Theorem 3.17), the Bony's paralinearization 
formula for a product, some simple computations and the commutator estimate 
proved in Section 2.4. 

Step 1: Paralinearization of the equation 

(d t + V ■ V)V + a( = 0. 

We begin by proving 
Lemma 4.11. We have 

(4.38) (d t + T v ■ V)V + T a ( + T ( (8 t + T v ■ V)B = h l 

for some remainder hi satisfying 

\\hi\\ B . < HWv\\ H s + ^ \\(V,B)\\ H .) jl + \\v\\ cr+h + 11(^,5)11^1 . 

Proof. Using (2.11) and (2.4) we have V ■ VV = T v ■ VV + A 1 where A 1 = 
Ej T d 3 vVj + R(djV, Vj) satisfies 

\\M\h* £ \\w\\ Loo \\v\\ Hs . 

Similarly, (a - g)( = T a _ g ( + T^(a - g) + R((, a - g) where 
(4-39) \\R{M\\ H .< ||C||^-i/2||a|| cl/2 . 

and where ||a|| c i/2 is estimated by means of (4.20). 
Since T^g = 0, by replacing a by g + (dtB + V ■ VB) we obtain 

T c a = T c (d t B + V-VB) 

= T c (d t B + T v ■ VB) + T ( (V - T v ) ■ VB. 

As in the analysis of A\ above, we have 

\\(V -T v )-VB\\ HS <\\VB\\ Loa \\V\\ Hs . 

Now we use \\T^\\ HS ^ HS < ||CIL°o ^ IMIirs+va (since s+ 1/2 > 1 + d/2) to obtain 

\\T ( (y - T V ) ■ VB\\ HS < \\rj\\ HS+i \\VB\\ Loc \\V\\ H . . 

This proves (4.38). □ 

Step 2. We now commute (4.38) with (D x ) s = (I - A) s / 2 . The paradifferential 
rule (2.5) implies that 

|| [To, (D x ) S ]\\ H s-l/2^ L 2 < ||a|| W l/2,oo , 
H s-1/2_ ¥L 2 < |[CI[ WVa.ao l 

l|[ryV,(^) s ]|| M2 <||y|| wl ,oo. 
Consequently, it easily follows from (4.20) and (4.38) that 

(d t + T v • V)(D X ) S V + T a (D x ) s C + T ( (d t + T v ■ V)(D X ) S B = h 2 
for some remainder h 2 satisfying 

\\h 2 \\ L 2 ^^(llr/ll^iJK^^IlHOli + ll^ll^ + IK^^IIcr}- 
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On the other hand, Lemma 2.16 implies that 
\\[T c ,d t + T v .V](D x yB(t)\\ L2 

< HUt)\\ H s +h ,\\(V,B)(t)\\H.) I 1 + Mt)\\ c r + l + \\(V,B)(t)\\ Cl 

Here we have used the fact that the L°° norm of dt( + V ■ V£ is, since s > 3 + 5, 
estimated by means of the identity (4.12): 

WdtQ + V • VC|| L oo < \\VB\\ LOO + lici^oo \\W\\ Laa 

< \\VB\\ Loe + \\ V \\ HS+ i \\VV\\ Loa . 

By combining the previous results we obtain 

(d t + T v ■ V)((D X ) S V + T C (D X ) S B) + T a (D x ) s ( = h 

where f\ satisfies the desired estimate (4.37). 

Step 3. Paralinearization of the equation 

(8 t + V- V)C = G{rf)V + (G(ri)B + 7. 

Writing (V - T v ) ■ V( = T V £ • V + Yfj=i R(djC, Vj) and using (2.11) and (2.12), we 
obtain 

(4.40) || (y - t v ) • vcil HS _i < l|vcil c _i livilfl. < \\v\\ c i \\v\\ HS . 

The key step is to paralinearize G{rj)V + (G(r])B. This is where we use the analysis 
performed in Section 3.3.4. By definition of R(rj) = G{rj) — T x we have 

G( V )V + CG(r])B = T X U + F 2 ( V , V, B), 

where 

(4.41) F 2 = [T c , T X ]B + R{r,)V + QR{rj)B + (C - T ( )T X B. 

The commutator [T^,T\]B is estimated by means of (2.5) which implies that 

\\[T c ,T x ]B\\ hS _i <{m (C)M 1 1 /2 (A)+M 1 /2 (C)A/ 1 (A)}||5||^. 

Since 

M° /2 (C) + Ml /2 (X) < IC{\\v\\ HS+i )\\v\\ c 4, and M °(C) + M^X) < IC(\\v\\ H s +i ) 
we conclude that 

IIITc.^ll^i^^Nlfl^i.ll^WNU- 

Moving to the estimate of the second and third terms in the right-hand side of (4.41), 
we use Theorem 3.17 to obtain that the H s ~2-novm of R{rj)V and R(rj)B satisfy 

\\R(v(t))V(t)\\ HS _ h +\\R(r,(t))B(t)\\ H s- h 



< T 



He+h ,\\(y,B)(t)\\ H s) {l+h(t)\\ c § + \\(V,B)(t)\\ C r} 



Since H s 2 is an algebra, the term (,R(r])B satisfies the same estimate as R{r])B 
does. It remains only to estimate (£ — T^)T X B. To do so we write 

(C - T C )T X B = T TxB ( + R(C, T X B). 

Thus (2.11) (applied with a = and f3 = s — 1/2) implies that 

\\(C-T c )T x B\\ HS _ h < \\T x B\\ c o \\C\\ HS _ h . 
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Using (2.4) this yields 

|| (C - T ( )T x B\\ hS _i < Mq 1 (A) \\B\\ cl HCll^x • 
We thus end up with 

l|*2|| HS -i < ^{\\v\\ H s +i A\(V,B)\\ HB ) [l + ||77|| c 3 + \\(V,B)\\ C: } . 

By combining the previous results, we obtain 

(4.42) ( dt + Tv .V)( = T x U + h 3 , 

where 

\\h3\\ Ha ^ < HM h s + ^\\(v,b)\\ hs ) [1 + iny + \\(v,B)\\ C r} . 

As in the second step, by commuting the equation (4.42) with {D x ) s we obtain the 
desired result (4.36), which concludes the proof. □ 

4.5. Symmetrization of the equations. Recall that 

U B := (D x ) s V + T C {D X ) S B, ( s := (D x ) s £ 

satisfy the system 

' (d t + T v -V)U s + T a ( s = f u 
^ {d t + Ty • V)Cs = T X U S + f 2 , 

where, for each time t 6 [0, T], 

W(h(t)j2(t))\\ L2xH -i 

< H\\v(t)\\ H s +i ,\\(V,B)(t)\\ H .) |l + Ht)\\ c r +h + \\(V,B)(t)\\ C r 

To prove an 1? estimate for System (4.43), we begin by performing a symmetrization 
of the non-diagonal part. Here we use in an essential way the fact that the Taylor 
coefficient a is a positive function. Again, let us mention that this assumption is 
automatically satisfied for infinitely deep fluid domain: this result was first proved 
by Wu (see [54, 55]) and one can check that the proof remains valid for any C 1,a - 
domain, with < a < 1. 

PROPOSITION 4.12. Introduce the symbols 

I a 



(4.43) 



7 = VaX, q l( ^ 



and set 9 S = T q Q s . Then 



(4.44) d t U s + T v ■ VU S + T 7 # s = F l , 

(4.45) d t e s + TvV9 s -T 7 U s = F 2 , 
for some source terms F\,F 2 satisfying 

\m(t),F 2 (t))\\ 



< T 

H 



, + 1 , \\(V,B)(t)\\ HS ) ll + Ht)\\ r + ^ + \WB)(t)\\cr} ■ 
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PROOF. Directly from (4.43), we obtain (4.44)-(4.45) with 
Fi :=f 1 + {T J T q -T a )Cs, 

F 2 := T q f 2 + (T q T x - T 7 ) U B - [T q , d t +Tv V] ( s . 
The commutator between T q and dt + Ty ■ V is estimated by means of Lemma 2.16: 

(4.46) \\[T q ,d t + T v -V]( s \\ L2{Rd) 

< K |Mo 3 (9) || F || cJ+e + + y . Vg ) [|F|| LOO (t)} x ||C S || H _ i (Rd) . 

T q f2 is estimated by means of (2.4). The key point is to estimate (T 7 T g — T a )£ s and 
(T q T\ — Try)U s . Since jq = a, the operator T 7 T g — T a if of order —1/2 since 7 is 
a symbol of order 1/2, q is of order —1/2, and since these symbols are m x. 
Similarly, since qX = 7, the operator T q T\ — T 7 is of order 0. More precisely, by 
using the tame estimate for symbolic calculus (see (2.5)), we obtain 



\T,T q - T a \\ H _^ L2 < Ml'l^M- 1 ' 2 ^) + M 1 / 2 (7)M- 1 / 2 ( g ), 



\\T q T x -T4 L ^ L2 < M 1/ 1 2 / \q)M^X)+M 1/2 {q)M{ /2 (A). 

The above semi-norms are easily estimated by means of the C 1 / 2 norms of Q = Vr; 
and a (given by Proposition 4.6). □ 

We are now in position to prove an L 2 estimate for (U S ,8 S ). 

Lemma 4.13. There exists a non- decreasing function J- such that 

(4.47) \\U B \\ L ao {[Q>nvt) + l|0 s || L oo ([OiT];L2) < T(M B>0 ) + VTT(M s (T))Z r (T). 

Remark 4.14. The fact that this implies corresponding estimates for the Sobolev 
norms of 77, ijj, V, B is explained below in §4.6. 



PROOF. Multiply (4.44) by U s and (4.45) by 9 S and integrate in space to obtain 

\u s (t)f L2 + \\e s (t)\\ 2 L*} + + = (ni), 



d_ 

dt 



where 



(I) 
(II) 
(HI) 



(T m • vu s (t),Us(t)) + (T v(t) • ve s (t),e s (t)), 

(T l{t) s (t),U s (t)) - (T l(t) U s (t),9 s (t)), 
(F 1 ,U 8 ) + (F 2 ,e s ). 



Then the key points are that (see point (Hi) in Theorem 2.7) 

\\(T V (t) ■ V)* + T v{t) ■ V\\ L2 ^ L2 < \\V(t)\\ W ^ , 

and 



We then easily obtain (4.47). □ 
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4.6. Back to estimates for the original unknowns. Up to now, we only 
estimated (U S ,6 S ) in L°°([0,T];L 2 x I?). In this section, we shall show how we can 
recover estimates for the original unknowns (rj, tp, V, B) in L°°([0, T]; H s+ ? x H s+ ? x 
H s x H s ). Recall that the functions U s and S are obtained from (rj,V,B) through: 

U s := (D X ) S V + TAD X ) S B, 



The analysis is in four steps: 



(i) We first prove some estimates for [B, V, rj) and the Taylor coefficient a in some 
low order norms. 

(ii) Then, by using the previous estimate of 9 S , we show how to recover an estimate 
of the L°°([0,T],# s+ ^)-norm of 77. 

(iii) Once rj is estimated in L°°([0, T], H s+ ? ), by using the estimate for U s , we 
estimate (B,V) in L°°([0, T]; H s ). Here we make an essential use of our 
first result on the paralinearization of the Dirichlet-Neumann operator (see 
Proposition 3.19). Namely, we use the fact that one can paralinearize the 
Dirichlet-Neumann operator for any domain whose boundary is in for 
some fj, > 1 + d/2. 

(iv) The desired estimate for i/j follows directly from the previous estimates for 
rj, V, B, the identity = V + BVrj and the fact that one easily obtain 
an L°°([0,r];L 2 )-estimate for tp. 

We begin with the following lemma. 

Lemma 4.15. There exists a non- decreasing function T such that for any r > 0, 

(4.48) \\r]\\ L ^ {[0) T];H^c:) + \\{B,V)\\ 1 , 

< T(M sfl ) + VTT(M s (T))(Z r (T) + 1), 

and, for any 1 < r' < r, 

(4.49) IMIiooflo^-i) < HMsfl) + VTF(M s (T))(Z r (T) + 1). 

Proof. The proof is based on the fact that it is easy to estimate the solution w 
of a transport equation of the form 

d t w + V - Vw = F. 

Indeed, by using the estimates (4.19)-(4.20) for o, tame product rules in Sobolev or 
Holder spaces and the identity dtrj + V ■ Vn = B, we readily obtain that there exists a 
non-decreasing function C (depending only on parameters that are considered fixed) 
such that 

\\a-g\\ 1 , = \\d t B + V-VB\\ 1 , <C(t), 
Kll A . 1 = \\d t V + V-VV\\ 1 , <C(t), 

\\d tri + V-Vri\\ C r nHs <C(t), 
H^a + V-Vall^-i <C(t), 



where 



C(t) = C(\\r,(t)\\ HS ^ , \\(V,B)(t)\\ Hs ) |1 + ||r?(t)|| c:+ x + \\(V, B)(t)\\ Cl 
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Let us come back to the proof of Lemma 4.15. To fix matters, we prove the estimate 
for V only (the proofs of the estimates for B, r/ and a are similar) and we begin by 
proving the Sobolev estimate. Using the obvious estimate 

H^lliHto.T]) ^ INL 2 ([0,T]) ' 
note that Fy := d t V + V ■ VU satisfies 

\\F V \\ i ! < VTT(MJT))Z r (T). 

Using Proposition 2.15 with a = s — ^, we estimate V in L°°([0, T]; H s ~z ). 

Let us now estimate the L°°([0,T]; Cy 2 )-norm of V. Let Ui = (L>I)U. We have 

$Vi + V ■ VVi = (Di)F{r + [V, (p\)]V x V. 
where according to Proposition 2. 14 [in], we have 

\\[v,{Di)]v x v\\ L oo < t(\\v\\ h °)\\v x v\\ l ~ < HWvysWWcr. 

As a consequence, we get 

(d t + V -V){D x ) l ' 2 V <VTT(M s (T))Z r (T). 

We next use the usual L°° estimate for transport equation (see (2.27)). This gives 
an estimate for {D x ) l / 2 V in L°°([0, T]; L°°) and hence the desired estimate for V 
in L°°([0,r];Cy 2 ). □ 

Lemma 4.16. There exists a non- decreasing function J- such that 

(4-50) IMUuo^+i) < H^Msfi) + TF(M S (T) + Z r (T))). 

Proof. Chose e and an integer N such that 

0<e<r-l, (N+l)e>-. 
Set R = I — Tx/qTq to obtain 

Cs = Ty g Tq(s + R( s , 
where recall that £ s = {D x ) s C- Consequently, 

( s = (I + R + ... + R N )T 1/q T q ( s + R N+1 ( S . 



By definition of q = Wa/X, Theorem 2.7 implies that, for all (i£R, there exists a 
non-decreasing function T depending only on e and inf( t i3; ) e [o,r]xR d a (^> x ) > such 
that, 

\\R(t)\\ m ^ m +e < F{\\a(t)\\ C e , \\r](t)\\ c i+,) , 

and 



\\ T l/q(t)\\ H „+l/2^ m < FUvWHw 1 

Therefore 

I|Vt7||^_i =IIC s || h -i <H\\ a \\ 

Now it follows from Lemma 4.15 that 

ll a llL°°([0,T];C*f) + IMlL°°([0,T];C, 1+e ) + IKs||i<x>([o,T];if-l) 

< T(M sfi ) + VTT(M s (T))Z r (T). 
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On the other hand, it follows from Lemma 4.13 that 

\\ T iCs\\ L ^ i[0 , W < -W.o) + VTF(M s {T))Z r {T). 
This implies the desired result. □ 

It remains only to estimate (V, B). 

Lemma 4.17. There exists a non- decreasing function J- such that 
(4.51) \\(V,B)\\ Loo{[0tT] . HS) < F{F{M sfl ) + TF{M S {T) + Z T {T))). 

Proof. The proof is based on the relation between V and B given by Proposi- 
tion 4.5. 

Step 1. Recall that U = V + T^B. We begin by proving that there exists a non- 
decreasing function T such that 

(4-52) H^Uuo^-i) < HHM sfl ) + TF(M S (T) + Z r (T))). 

To see this, write 

(D x ) s -l U = {D x y\ {U s + [{D X ) B ,T C ]B} 
and use Theorem 2.7 to obtain 

\\m\T ( ]B\\ H _ k <u\\^\\B\\ HS _ h . 

Since, by assumption, s > 3/4 + d/2 we have 



and hence 



\ c \ < HCll^-i < ||r/|| HS+ i 



^ s _i < \\u 8 \\ H _ k + h\\ HS+ i \\B\\ $ ■ 



The three terms in the right-hand side of the above inequalities have been already 
estimated (see Lemma 4.13 for U s , Lemma 4.15 for B and Lemma 4.16 for rj). This 
proves (4.52). 

Step 2. Taking the divergence in U = V+T^B, we get according to Proposition 4.5, 
Lemma 4.15 and Lemma 4.16: 

div U = div V + div T f B = div V + T div C B + T ( • VB 

= -G{r,)B + T K<+dlvC B + 1 

= -T X B + R{r,)B + T iC<+diYC B + 7 

= T q B + R{7 1 )B + T d ^ c B + 1 

where, by notation, 
and 

According to Proposition 3.19 (with \i = s — |) and Lemma 4.15, we deduce 
(4.53) T q B = divU- T div( B - R( V )B - 7. 

Now write 

B = TiT q B + ( I - TiT g )B 

q \ q J 

to obtain from (4.53) 

B = Ti dwU-Tij + R £ B 

q q 
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where 

(4.54) R_ e := Ti ( -T divC - R(rj) ) + [I-T l T q ). 



3 

Notice now that according to Lemma 4.16, we control div£ = Ar) in H s ~2, and 
since s > 5 + 5, there exists e > (here one can fix e = 1/4) such that T^ v ( is an 
operator of order (both Sobolev and Holder) 1 — e. Finally, q = — A + i£ • £ £ with 
M^(q) < C(\\r)\\ + i) if s > \ + | + e. Moreover, is of order —1 and we have 

M^cT 1 ) <^(IMI^+i)- 

Consequently, according to (2.4) and (2.5), the operator given by (4.54) is of 
order -e. applying T { _ XirK , i y 1 to (4.53), we get 

B = W + R- e B 

where W := Ti div U — Ti 7 satisfies 

IIWUh. < F[F{M 8fi ) +TF[M S {T) + Z r {T))), 

and is an operator of order — e. Since we have already estimated the H s ~ 2-norm 
of B (see Lemma 4.15), we conclude by a classical bootstrap argument that 

\\B\\ H > < ?[F(M afi ) +TF[M S {T) + Z r (T))), 

and coming back to the relation U = V + T^B we get that V satisfies the same 
estimate. □ 

Lemma 4.18. There exists a non- decreasing function T such that 

I^U([o,T];^) ^ H^Msfi) + TF{M S {T) + Z r (T))). 

Proof. Since Vtp = V+BVn and since the L°° ( [0, T] ; H s ~ I )-norm of (V77, V, B) 
has been previously estimated, it remains only to estimate HV'IIl^Qo t]-l 2 )- 

Since 

Vr? • VV» + G(r/)V> 
: " 1 + lVr?! 2 ' 

the equation for ip can be written under the form 

+ 9V + \ |V^| 2 - i(l + iVr?! 2 )^ 2 = 0. 
Therefore, since V = VV> — BVn, 

d t ^J + V-Vip = d t %b + (V^l 2 - BVn • Vip 

= ~9V + \ |VV| 2 + ^(1 + iVr?! 2 )^ 2 - 5V?? • VV 

(4.55) 1 9 1 o 9 1 90 

= -grj+~ |W - BV V \ 2 - -B 2 \Vn\ 2 + -(1 + \V V \ 2 )B 2 
= -gr 1 + l -V 2 + l -B 2 . 
The desired 1? estimate then follows from classical results (see Proposition 2.15). □ 
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5. Contraction 



In this section we prove the uniqueness of solutions. 

Theorem 5.1. Let j = 1,2, be two solutions of (1.6) such that 

{ i/j . i. 'j ■ \ j ■ I > j i G C°([0,T ];tf s+ 5 x H s+ ^ xH s x H s ), 

for some fixed Tq > 0, d > 1 and s > 1 + d/2. We also assume that the condition 
(1.2) holds for < t < To and i/iai i/iere exists a positive constant c such that for 
all < t < Tq and for all x G R d , we /lane a 3 -(i, x) > c /or j = 1, 2, t G [0, T]. Se£ 

■q := rjx - r) 2 , ip = ipi - ip 2 , V := V\ - V 2 , B = B 1 -B 2 . 
Then we have 

(5.1) || (t?, ip, V, B)\\ Lx mTyHS _ hxHS _. x ^ s _ 1 xHs _ 1} 

<K(Mi,M 3 )\\(ti,i/>,V,B) \ t=0 



' H s ~i xH s ~i xH s ~ 1 xH s - 1 ' 

Let us recall that 

' {d t + Vj ■ V)Bj = a, - g, 

(5.2) J (d t + Vj ■ V)Vj + ajCj = 0, 

. (dt + ^ • V)0 = Cin^Vj + QG^Bj + 7i , = VVj, 
where jj is the remainder term given by (4.10). Let 

N(T) -^II^^V.^WII^^j^^. 

Our goal is to prove an estimate of the form 

(5.3) N(T) < IC(Mi,M 2 )N(0)+TIC(Mi,M 2 )N(T), 

for some non-decreasing function K, depending only on s and d. Then, by choosing T 
small enough, this implies N(T) < 21C(Mx,M 2 )N(0) for T\ smaller than the mini- 
mum of To and l/2/C(Mi, M 2 ), and iterating the estimate between [Ti, 2Ti],. . . , [T — 
Ti,Ti] implies Theorem 5.1. 

Remark 5.2. Notice that we prove a Lipschitz property in weak norms. This is a 
general fact related to the fact that the flow map of a quasi-linear equation is not 
expected to be Lipschitz in the highest norms (this means that one does not expect 
to control the difference (r],ip,V,B) in L°°([0, T ]; H s+ ^ x H s+ t x H s x H s )). 

5.1. Contraction for the Dirichlet-Neumann. The first step in the proof 
of Theorem 5.1 is to prove a Lipschitz property for the Dirichlet-Neumann operator. 
This was already achieved in a very weak norm in Theorem 3.10, and here we used 
elliptic theory to improve the result. 

Theorem 5.3. Assume that s > 1 + 5. There exists a non- decreasing function T 
such that, for all 771, 772 G H s+ z and all f G H s , we have 

(5.4) \\[G( m ) - G( m )\ f\\ HS _, < ^(11(^, 772)11^1) Ht?! - ||/||^ . 

Remark 5.4. In the right-hand side, || 771 — 772 1| s _i controls only the C2- norm of 
771 — 772. This is the source of some serious difficulties in the proof below. 
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(5.5) 



Proof. The proof follows closely that of Theorem 3.10 and we keep the nota- 
tions pj,(pj,v = (p\ — 4>2 , A 3 introduced there. 

Notice that, using the smoothing property of the Poisson kernel, we have 

( (i) Al - A| = /3 k d z , with supp/? fc C K d x J, where J = [-1,0], 

\ (»*) \\Pk\\l?(J,H-i(l!L*)) < H\\{Vl,V2)\ HB+ i xH s + i)\\Vl-m\\ H s-i (Rd y 

Recall that 

(5.6) G( m )f = Ujl^o, Uj = A{fa - V xPj ■ A$j. 

Let us set U = U\ — U 2 - According to (3.14), Theorem 5.3 will follow from the 

following estimate 

(5.7) 

\\U\\l2{j,h°-i) + \\d z U\\ L 2^ HS - 2) < 7/2)11^+1 xH s+i) ll/llfl. \\m - V2\\ HS -i ■ 

According to (3.37) and (3.39) the estimate (5.7) will be a consequence of the fol- 
lowing one 
(5.8) 

5 

^2\\Bkh 2 (j,Hs-i) < F(\\(vi,m)\\ H ,+$ xH s+$) WfWw II 7 ?! -m\\ H s-$ where 
/,• i 

B 1 = A\v, B 2 = {V X:Z p 2 )A 2 2 v, B 3 = (A\ - A{)4> 2 , B A = V x , 2 (pi - p 2 )A$i, 
B 5 = (V x , zP2 )(Al - Al)fa. 

Since <f>j is a variational solution, Proposition 3.22 with a = s — 1 show that 

l|V a! ,*^-|Uoo(j i H.-i) + ||A^ J -|| £ o (j )H .-i) < ^(WvjW^+^WfW^- 

Since s > 1 + |, it follows the from (3.34) that 

5 

(5.9) Yl \\ B i\\&V,H-i) < F(\\(Vi,rt2)\\ HS+ i xHS+ i)\\vi ~ mWp-l \\f\\H>. 

1=3 

Since 

II-SiIIl^j,^- 1 ) + \\ b 2\\l*(j,h°-->-) < J r (\\(vii'n2)\ HS+ i xHS+ i)\\^x,zv\\ L 2( tItH s-i- ) 

using the estimate (5.9), we see that (5.8) will be a consequence of the following 
Lemma. Therefore Theorem 5.3 will be proved if we prove the following result. 

Lemma 5.5. We have 

\\Vx,zv\\l2(j,h*-i) < ^(\\(vi,V2)\\ HS+ i xHS+ i)\\vi -m\\ H s-^\\f\\H^ 

Proof. Notice that v = <f)\ — (f) 2 is a solution of the problem 

(5.10) d 2 z v + a 1 Av + (3 1 -Vd z v--/ 1 d z v = F, v\ z=0 = 
where 

F = (a 2 — «i)A^ 2 + (#2 - /3i) • Vd z 4> 2 - (72 - 71)^2 

and otj are given by (3.18). We would like to apply Proposition 3.22 with a = s — |. 
To this end, according to (2.45), we shall estimate the L 2 ( J, H s ~ 2 (R d )) norm of F 
and the X _ a(J) norm of V xz v. 
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Estimate on F: Since s > 1 + i (thus 2s — 3 > 0) we may apply (2.16) with s± 
s — 2, S2 = s — 1, So = s — 2. We get 



(ai — «2)A02 
! >i -/? 2 )-V9 2 02 

(71 - ll)dz4>2 



L 2 (J,H S ~ 2 ) 



L 2 (J,H S ~ 2 ) 



L 2 (J,H*- 2 ) 



<K\\ai-a 2 \\ L 2 {J)HS -i ) 
<K\\jx -72|li2 ( j ji?s -2) 



A0 2 



5; 



z92 



L°°(J,H S - 2 ) 

L°°(J,H S - 2 ) 
L°°(J,H S - 1 ) ' 



Then, using the product rule in Sobolev space (2.16), and (3.11), (3.12) we obtain 

(5.11) ||qi - a2\\ L2{ j tH s-i) + HA - A h 2 (J,H°-i) + hi -72|Il2 WHs - 2) 

< ^(llfri^H^+i^+i) Hr/! -r/ 2 || ffS _i . 
Moreover from from Proposition 3.22 with a = s — 1 we have 



It follows that 

(5.12) \\F\\ L 2 {J , m - 2 ) < HMvi, m)\\ H s +i ) \\m -m\\ H s- h 

Estimate of \\V X)Z v\\ x _x , J = (-1,0). 

We claim that 

(5.13) \\V x , z v\ 



H" 



X-2(J) 



<HKm,m)\\ 



H"+?xH s+ ? 



Since 



Uj + / we have v = u\ — u%. We begin by proving the following estimate. 



There exists a non decreasing function T : R + — > R + such that 

(5.14) ||v^|| i2(J)I/2) <T{\\(m,V2)\\ HS+ i xHS+ i)\\vi-V2\\ HS ^ 

For this purpose we use the variational characterization of the solutions Uj. Setting 
X = (x, z) we have 



H s 



(5.15) 



tfui • A l 6 J* dX 



A7 • A* 6 Ji dX 



for all 6 £ H l ' Q (Q), where J< = |«9 z/ 9;|. 

Making the difference between the two equations (5.15), using (3.32) and taking 
= v = u\ — U2 one can find a positive constant C such that 



[ (A^l 2 dX < C(A 1 + A 2 + A 3 + A 4 ) 
./n 



where 



Ai= / ((A 1 - A 2 )n 2 ||A 1 t;| J\ dX, A 2 = {(A 1 - A 2 )u| |A 2 S 2 | J\dX, 



A, 



|A 2 S 2 ||A 2 w| |Ji - J 2 |dX, 



A 4 = / I (A 1 - A 2 )f\\A 1 u\ Ji dX, 



A 5 = I (A 1 - A 2 )u||A 2 /| Ji dX, 



'_|AV||A^||Ji- J 2 |dX 



It follows from the elliptic regularity theorem that 

Ai < IIA^ii^^ii^ii^^iia^iiLoc^cx,^)) 



.i^lbll^+i^,!!' ii,/-, 
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».<*)) Il*7i -m\\ 



Noticing that A 1 — A 2 = (3(d z pi)A\ where (3 satisfies the estimate in (3.34) we obtain 

A 2 < ll^pill^^ll^ll^^llA^sll^^HA^II^^. 
Using (3.32), (3.34) and the elliptic regularity we obtain 

A 2 < \\A 1 v\\ L2(fi) T(\\(r]i,V2)\\ HS+ ^ xHe+ i) Ihi-^H^.! 
Now we estimate A3 as follows. We have 

A 3 < ||A 2 u2|| Xoo(S) ||A 2 u|| i2(S) ||Ji - J2|| £ a ( n)- 
Then we observe that 

\\Jl- J2\\ L ^ ) <C\\r 1l - m \\ Hh{nd) 

||A 2 t;|| i2(S) < CIIA^H^^ 

and we use the elliptic regularity. To estimate A4 and we recall that / = e z ^ Dx ' f. 
Then we have 

H^/llL 2 (JxR d ) < H/ 3 llL 2 (/xR d )ll^/llL°°(JxR d )- 

Since ||d z /|| L oo (JxR[i) < ||9 z /||l-»(j,^-i(r^) < ||/||h*(r*) s using (3.34) we obtain 

a 4 + a 5 < IIA^H^^^dl (771, 7/2)11^+1) Wm-mWjr-t II/IIh.- 

The term A§ is estimated like A3. Since ^ < s — | this proves (5.14) . 
To complete the proof of (5.13) we have to estimate ||Va;,zU|| ioo l H -\ ■ The esti- 
mate of llVrull 1 follows from (5.14)and from Lemma 3.14. To estimate 

11 n L°°(J,_H~^) v ; 

^ zV ^l°°(j H~i) We ^ ave *° use (^•^) an< ^ * ne ec L ua ti° n satisfied by v. If we prove 
that 

(5.16) \\8%v\\zi(j jH -i) < ^(|| (771, ^ll^+^^+i) \\vi - mWjje-i \\f\\ H ° ■ 

the result will follow again from Lemma 3.14. Recall that v satifies the equation 
(5.10). 

It follows that we have 

^ 5 17 « 11^^11^2(^-1) < HaiAvHia^-i) + • Vd z v\\ L 2 {J:H -i) 

+ ll7i^«IU 2 (J,H-i) + \\F\\l2(j,H-iy 

Since -1 < s - 2 (5.12) yields 

\\ f \\li(j,h-i) < \\ f \\l\j,h*-i) <^{\\{ , num)\\ H s+i xHS +\) \\vi-m\\ H s^ \\f\\ H - ■ 

On the other hand, since s — \ — 1 > and — 1 < s — \ — 1 — | (2.16) show that we 
have 

\\<xi*»\\i*(.J,h-v) < Ilaill^^-^IIV^H^^) 
lift •V5 z .«| £J (j iH -i) < H Loo(J ^.-jJI^«llii»(j,r») 

1171^1^2(^-1) < 117111^(^-3^^1^2(^2). 

Using Lemma 3.25 and (5.14) we obtain eventually (5.16). 

Now Lemma 5.5 follows from (5. 12), (5. 13) and Proposition 3.22 with a = s— |. □ 

Finally Lemma 5.5 together with (5.8) prove (5.7) wich in turn proves Proposition5.3. 

□ 
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5.2. Paralinearization of the equations. We begin by noticing that, as in 
the proof of Lemma 4.18, it is enough to estimate rj,B,V. Indeed, the estimate 
of the L°°([0,T];H s x ' 1/2 ) -norm of ip is in two elementary steps. Firstly, since Vj = 
Vipj — BjVrjj, one can estimate the L°°([0, T]; i? s_3//2 )-norm of from the identity 

= V + BVr^ + B 2 Vr]. 

On the other hand, the estimate of the L°°([0, T]; L^)-norm of ifi follows from the 
equation (4.55). 

An elementary calculation shows that the functions 

C = Cl-C2, V = Vi-V 2 , B = B 1 -B 2 
satisfy the system of equations 

( d t B + Vi • VB + V ■ VB 2 = a, 

8 t V + Vx ■ W + V ■ VV 2 + a 2 C + aCi = 0, 
. dtC + V 2 • VC + V • VCi = G( m )V + CiG(m)B + (G( m )B 2 + R + 7, 
where 

(5.19) R = [G( m ) - G( m )} V 2 + Ci [G( m ) - G( m )\ B 2 , 

and 7 = 71 — 72, 7j are given by (4.10) 

Lemma 5.6. The differences £, B, V satisfy a system of the form 

f (dt + Vi ■ V)(V + &B) + «2C = /i, 
( ' j \ (ft + V 2 ■ V)C - G(t/i)F - (iG( m )B = f 2 , 

for some remainders such that 



(5.18) 



Proof. We begin by rewriting System (5.18) under the form 
( d t B + Vi ■VB = a + R 1 , 

8 t V + V 1 ■ W + a 2 ( + aCi = i? 2 , 
^ dtC + V 2 -X7( = G( m )V + (iG( m )B + R + 1 + R 3 , 

where R is given by (5.19), 7 = 71 — 72 and 

R 1 = -V-VB 2 , R 2 = -V-VV 2 , R 3 = V ■VCi + CG{t] 2 )B 2 . 
From Theorem 5.3 one has 

\\R\\ , 3 < /C(Mi,M 2 )iV(r). 

Similarly, proceeding as in the end of the proof of Proposition 4.3, we have 

On the other hand, since s — 1 > d/2, fP" 1 is an algebra and 

\\V ■ V-B 2 ||fl»-i < K \\V\\ H s-i ||VB 2 ||fl.-i < K \\V\\ H s-i \\B 2 \\ HS 
and similarly 

\\V ■ VV 2 \\ H .-i < K \\V\\ HS ^ \\V 2 \\ HS . 
On the other hand, according to Theorem 3.18 we have 

\\G( m )B 2 \\ H .-i < C(\\ m \\ HS+1/2 ) \\B 2 \\ HS , 
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and hence 

\\(G(r, 2 )B 2 \\ HS _s < C(\\ m \\ HS+1/2 ) \\B 2 \\ H . ||C|| H8 _| • 
To estimate V • V£i we use the product rule (2.16) to deduce 

r-vci|| ffs _ f <ir||V|| H .-i HVCill^i <K\\V\\ H ^ HW^+i- 
Therefore we have, 

\\Ri\\ H s-i + WMh-i + ll^ll^-a < c{|M| hS _i + \\b\\ hs ^ + , 

for some constant C depending only on ||^-||^ s+ i , ||-Bj||jj s , ||V}[|jy s . The next step 
consists in transforming again the equation. We want to replace a£i in the second 
equation by 

(d t B + V 1 -VB-R 1 )C 1 . 

The idea is that this allows to factor out the convective derivative df + V% ■ V . 
Writing 

(8 t B + Vi ■ VB)Ci = (d t + Vi • V)(BCi) " B(d t + Vi • V)Ci 
we thus end up with 

(5.21) (dt + Vi ■ V)(V + CiB) + a 2 ( = i?iCi + B{dt + Vi • V)Ci + R 2 . 
Since 

(d t + Vi • V)Ci = G(rn)Vi + C,iG(r]i)Bi + 71, 

we have 

IK^ + vi-vKill^-x < T(\\( m ,B 1 ,v 1 )\\ HS+ixHsxHs ). 

By using this estimate and our previous bounds for R\,R 2 , we find 

\\RiCi + B(8 t + Vi • V)Ci + R 2 \\ H s-i < C {lhll ff s-i + + llVH^-i} , 

for some constant C depending only on [I^H , , l|V}||# s - Notice that here, 

as we used the equation satisfied by £i, it was important to have (dt + Vi • V) in the 
l.h.s. of (5.21) and not (dt + V 2 ■ V), and this algebraic reduction required some care 
in the previous step. □ 

5.3. Estimates for the good unknown. We now symmetrize System (5.20). 
We set 1= [0,T]. 

Lemma 5.7. Set 

£:=y/X^, ^:=T^(V + CiB), •& := T^(. 

Then 

(5.22) (d t +T Vl -V)<p + T e # = gi , 

(5.23) (dt+T V2 -V)#-Ti<p = g 2 , 
where 

ll^'^H^c/^-Sx^-l) ^ KWi,M 3 )N(T). 

Proof. We start from Lemma 5.6. By using Proposition 2.11, one can re- 
place Vi • V by Ty 1 ■ V and a 2 C by T a2 (, modulo admissible remainders. It is found 
that 

(5.24) (dt + T Vl ■ V)(V + CiB) + T a2 ( = f[, 
for some remainder /{ such that 

||/ilL=c (J . ff .-i) <JC(M U M 2 )N(T). 
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Similarly, one can replace V 2 -V by Ty 2 -V. According to Proposition 3.19, with e = |, 
we have 

\\G( m )V - TmVW^^ + \\G( m )B - T Xl B\\ Lao{I . H s-^ < W^T), 

and according to Proposition 2.11, with 7 = r = s— |,/i = s— |, 
\\CiG( Vl )B - r Cl T Al B|| J . oo(j5fl ._j ) < /C(M!)iV(T). 

We deduce 

(5.25) + Ty 2 • V)C - T Al F - T (l T AlJ B = & 

where 

and 



Ai := Va + \Vvi\ 2 M 2 -(Vvi-0 2 , 



\\f2\\ Laam s-^<^(Mi,M 2 )N(T). 
Now, according to Lemma 2.17, (3.43) and (4.10) and we find that 

(5.26) ||[r Ar ,(o t + r Vl -v)]|| ifs _ 1 ^ s _3 

< £(M x )(.M*(\A0 + M$((dt + Vx • V)v/a7)) < K'iMt) 
and similarly, according to Lemma 2.17 and (4.20), 

(5.27) ir^^ + T^-V)]^^^ 

< K(M 2 )(M° (^) + M° ((d t + V 2 ■ V)V^)) < K!{M 2 ). 

which implies 

(5.28) (d t + T Vl ■ V)T^(y + CiB) + T^T a2 C = fl 

(5.29) (d t + T V2 ■ V)T^( ~ T^{T Xl V - T (l T Xl B) = f», 
where 

W(fiJ2)\\ LX{I . HS -^< JC(M 1 ,M 2 )N(T). 
According to (2.5), (3.43) and (4.19), since s > 1 + | , 

T VM T *2 ~ T VM^ T V^ is of order 0, 



which implies (5.22). On the other hand, according to (2.5) and (3.43) the op- 
erators T^ 1 T\ 1 — Ta 1 ^ 1 and TajT^ — are of order 1/2 (with norm controlled 
by /C(Mi), which allows to commute T/^— and T a2 in (5.29)). Now, according 
to Proposition 2.11 (with 7 = r = s— ^,fi = s — 1) 

||r ClJ B-Ci5|| H8 _i ^(MOHAH^-i. 

Which implies (5.23) (using again (2.5)). □ 
Recall that we have set 

(5-30) ^^^PlI^^V.BJWII^i^j^^. 
Lemma 5.8. Se£ 

AT'(T) := sup {||^)|| HS _a +11^)11^3}. 
(5.31) N'(T) <IC(M 1 ,M 2 )(N(0) + TN(T)). 

68 



Proof. We first prove that 

(5.32) N'(T) <K{Mi,M 2 )(N , (0)+TN{T)+TN'(T)). 
The desired estimate (5.31) then follows from the fact that 

(5.33) N'(0) < )C(M 1 ,M 2 )N(0), N'(T) < K(M 1 ,M 2 )N(T), 

which follows from the continuity of paradifferential operators in the Sobolev spaces 
(see Theorem (2.7)) and the fact that H s ~ 1 (R, d ) is an algebra since s > 1 + |. 

The proof of (5.32) is based on a classical argument : we commute {D x ) s ~ 3 / 2 
to (5.22)-(5.23) and perform an L 2 estimate. Then the key points are that (see 
point (iii) in Theorem 2.7) 

{5U) IK^-vr+^-vii^^ii^iu., 

\\ T e ~ (Ti)*\\ L 2^ L 2 < FiKvi, m)\\w*K°°) 

and that the commutators [Ty^ • V, (D x ) s ~ 3 / 2 ] are, according to (2.5), of order s— |. 

Notice that since (</?,$) £ C 1 ([0, To]; H s ~2 ), we do not need to regularize the equa- 
tions. □ 

Finally, let us notice that an elementary argument allows to control lower norms 
of (V, B) (and hence also of V + Ci-B): 

(5.35) 11(^)11^.^-1) </C(M 1 ,M 2 )(iV(0) + TiV(T)). 

Indeed, (the proof of) Theorem 5.3 implies that (with a = a\ — a 2 ) 
(5-36) IHU^-f) < AC(M l9 Af 2 )JV(r). 

Since 9 4 S + Vi • V.B = a — V • V5 2 , we have 
(5.37) 

\\B\\lo°{I;H°-*) < \\B(Q) \\ H -*+ J (llVi-VBllj.-a + Hallfl.-a + llV-V^llH.-a)^ 

< \\B(0)\\ HS -2 +TJC{M 1 ,M 2 )N(T). 

Similarly, we have 

(5.38) \\V\\ l «> ( i ; h-») < \\V(0)\\ H s-2 +T1C(M 1 ,M 2 )N(T). 
Now we have 

V + CiB = T^p—icp + (Id - T^T^XY + (iB), 

where according to (2.5), the operator Id — T ^—iT^j^ ls °f order— 1/2. Hence, we 
deduce from (5.33), (5.31), (5.35), (5.38) and a bootstrap argument 

(5.39) \\V + &B|U<» ( j ; jf.-i ) < /C(Mi,M 2 ){iV(0) + TN(T)}. 

5.4. Back to the original unknowns. Recall that I = [0, T] (resp. J = 
(—1,0)) is an interval in the t variable (resp. in the z variable). 

Lemma 5.9. 

(5-40) I^Ud^-i) ^ W,M 2 ){iV(0) +TN(T)}. 
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Proof. From the equation dtijj = G(r]j)if)j we have, 

77(t) = 77(0) + f G( m )^(t')dt' + [ (G( m )-G(rj 2 ))Mt')dt', 
Jo Jo 

from which we deduce according to Theorem 5.3, 

(5.41) \\v\\l^(I;H^) < \\rim\H*-* +TK{M 1 ,M 2 )\\r 1 \\ Loo{pHS _ h) . 

Let R = Id — T 1 T/^, which, according to (2.5) and (4.19) is an operator of 
order — ^ (with norm estimated by )C(M 2 )). We have 

V77 = RVn + T 1 d. 

Therefore we deduce from (5.41), (5.31) and a bootstrap argument, 

l|V??ll L-(/;^-§) " ^(M 1 ,M 2 )(N(0)+TN(T)). 
Combining with (5.41) gives Lemma 5.9. □ 

We are now ready to estimate (V, B). 
Proposition 5.10. 

(5.42) || (V, 5)11x00(^.-1) <)C(Mi,M 2 ){N(0)+TN(T)}. 

The proof will require several preliminary Lemmas. We begin by noticing that it is 
enough to estimate B. Indeed, if 

< /C(M l5 M 2 ){iV(0) + TN(T)}. 

then, by using the triangle inequality, the estimate (5.39) for V + CiB implies that 
V satisfies the desired estimate. 

Let v = — <j>2 1 where (fij is the harmonic extension in f2 of the function ipj and set 

, 9z4>2 „ 
b 2 := q , w = v-T b2 p. 

OzP2 

We notice that 

(5.43) w\ z=0 = ip - Tb 2 t]. 
We first state the following result. 

Lemma 5.11. We have 

(5.44) U-T B2 v\\l°°(i-h°) <K(Mi,M^){N(0)+TN(T)} 

Proof. Indeed, the low frequencies are estimated by (5.35), while for the high 
frequencies, we write 

V(V> - T B2 rf) = Vip - T B2 Vr) - T^ B2 r] 

= Vipi - ViJj 2 - T B2 Vt] - TvB 2 rj 

= V X + £iV?7i -V 2 - B 2 V m - Tb 2 Vt] - T V b 2 V 

= V + (B 1 - B 2 )V m + B 2 (V Vl - Vi l2 ) - T B2 Vi] - T VB2 V 

= V + (iB + (B 2 - T B2 )V V - T V b 2 V, 

where we used that, by definition, Vipj = Vj + BjVrjj and Ci = V771. 

The main term V + Q\B is estimated using (5.39), while the two other terms are 
estimated using (5.40), the a priori estimate on B 2 and the product rules (2.9) 
and (2.12). □ 
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We next relate w, p and B. 
Lemma 5.12. We have 



B = 

Proof. Write 

Bi — B 2 



1 



d z pi 



d z w - (b 2 - T b2 )d z p + T dzb2 p 



2=0 



0: 



z9l 







z92 



z=0 



d z pi d z p2 

a — (d z 4>i - d z (f>2) + { 
d z pi 

1 



1 



1 



d z pi d z p 2 







z92 



z=0 



1 d 



z92 



d z p 



d z pi d z p\ d z p 2 

and replace v by w + T b2 p in the last expression. 



z=0 



□ 



Lemma 5.13. Recall that 6 2 := tFOt- For k = 0,1, 2, we have 



i < C 1 1 V>2 1 1 4-1 

C°([-l,0],L°°(/,H s -2- fc )) " r U H S+ ? 



for some constant C depending only on H772II s+ i • 

Proof. We estimate V,,^ in C°([-l, 0], L°°(I, H s ~2 )) by using the elliptic 
regularity (see Proposition 3.22 and Remark 3.21). Now, using the equation satisfied 
by </>2 and the product rule in Sobolev spaces, we successively estimate d z 4>2 and 
This proves the lemma since the derivatives of P2 are estimated directly from 
the definition of p2- □ 

Notice that 77 and hence p are estimated in L°°(I;H s ~z) (see (5.40)). Now, use 
Lemma 5.13 and Proposition 2.11 (applied with s > 1 + d/2, 7 = s — 1, r = s — 1/2, 
p = s — 3/2) to obtain 

||(6 2 -T b2 )d z p\\ Hs - 1 < \\b2\\ HS -i \\v\\h°-i • 

Now, (2.15) implies that 

\\ T dzb 2 P\\H s - 1 ~ IIMIifs-1/2 IMI/p-i , 

and hence, to complete the proof of the Proposition 5.10, it remains only to estimate 
d z w\ z= Q in L°°(I, .fP - " 1 ). This is the purpose of the following result. 

Lemma 5.14. Fort G [0,T] we have 



(5.45) 



\V x ,zW\\c° ([-1,0] 



, Hs -i) < JC(Mi,M 2 ){N(0) + TN{T)}. 



Proof. To prove this estimate, we are going to show that w satisfies an elliptic 
equation in the variables (x, z) to which we may apply the results of Proposition 3.22. 
We have 

d 2 z v + aiAw + fa ■ Vd z v - jid z v = (71 - 7 2 )9 2 02 + F x , 
where (see (5.10)) 

Fi = (a 2 - «i)A0 2 + {fi 2 - h) ■ Vd 2 2 - 
We claim that for t G [0, T] 



(5.46) 



\ F ^-)\\ L2{JiH s-^ < JC(M 1 ,M 2 ){N(0)+TN(T)}. 
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The two terms in F% are estimated by the same way. We will only consider the first 



one. Using the product rule (2.16) with sq 
write for fixed t 



§>si 



1,S 2 



s — | we can 



a 2 - cki ) A0 2 1 1 ^ 2 ( j. )jH - s _ § ^ < C\\a 2 - ax\\ L ^(j,H s - l )\\ A h\ 



L°°(j,H s ~^y 



Then we use (5.11), Proposition 3.22 with a = s — \ and Lemma 5.9 to conclude 
that the term above is estimated by the right hand side of (5.46). 

Now we introduce the operators 



7j 



Pj :=d 2 z + aj A + p r vd z , 

With these notations we have jj = Q^-PjPj and 

L\v = (71 - 7 2 )<9 2 </> 2 + F\. 



0, 



(i = i,2). 



(5.47) 
Moreover 

7i - 72 



1 



-PxPi 



1 



-P2P2 



1 



d z P2 



d z pi d z p 2 

' PXP +J—(Pl- P2 )P2 + 



P1P1 + 



1 



1 



d z P2 d z p 2 

1 / 1 

OzP2 



1 



1 



where 
(5.48) 



F, 



d z px d z p 2 
1 



d z px 
PxPi + F 2 , 



d z p\ d z p 2 

O z p2' 



PlPl 



d z P2 



P2P2 



d z p2 



(ax - a 2 )Ap 2 + (/3i - /3 2 ) ■ Vd z p 2 



Now we observe that 



1 



1 



which implies 



d z px d z p 2 
1 



PxPi 



z p PxPi 



d z p 



d z p2 d z px d z p 2 



7i, 



7i - 72 



-PiP 



d z p 
d z p2 



71 + F 2 



d z p 2 



d z p2 

Plugging this into (5.47) yields 

(5.49) Lxv - b 2 (L lP ) = Fx + (d z fc)F 2 

We claim that for fixed t we have 



L lP + F 2 . 



(5.50) 



L2(JHS _ §) < JC(Mx,M 2 ){N(0) +TN(T)}. 



{d z <h)F 2 {t 

Indeed we first use the product rule (2.16) to write 

By the elliptic regularity the first term in the right hand side is bounded by K(M2). 
It is therefore sufficient to bound the second one. We have, for fixed t 



1 



d z P2 



(ax - a 2 )Ap 2 



'L°°(J,H S ~?) 



< K,(M 2 )\\ax - a 2 || I; 2 (i7)Hs -i)||Ap2| 

l 2 {j,h s -?) 

Using (5.11) and (3.11) we see that the right hand side is bounded by the right hand 
side of (5.50). The second term in F 2 is estimated by the same way. 

To estimate v — T^p we paralinearize in writing 



(5.51) 



h(Lxp) = T b2 Lxp + T Llp b 2 + F 3 . 
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We claim that for t 6 [0, T] 

(5.52) ||F 3 (i, Oll^^-I) < /C(M l5 M 2 ){iV(0) + TiV(T)}. 

To prove it we shall use (2.9) with a = s — |,/3 = s — 2. Then o + /3-f>s-|.It 
follows that, for fixed z and i we have 

\\F 3 (t,;z)\\ H ^s <C\\b 2 \\ HS _i\\L lP \\ Hs -i. 

Therefore 

Now as we have seen before we have II &2 (*•)!! « i < IC(M 2 ) and due to the 

ll *\ '"L° a (J,H a ~Z) - y ' 

smoothing of the Poisson kernel ||-^ip||l 2 (j,J3" s - 2 ) — ^(^i)IMI The estimate 

(5.52) thus follows from (5.40). 

Setting F 4 = Ti ip b 2 we claim that for fixed t we have 

(5.53) \\F,(t, •)ll i2WffS -3 ) < IC(M 1 ,M 2 ){N(0) + TN(T)}. 

To see this we use (2.15) with so = s— |,si = s — 2, s 2 = s — \. We get 

Wm^W^j^ <ll^lKti-)llL»(J,H.->)ll&2(^-)ll LOO(J)fl .-i ) 

and (5.53) follows from estimates used above. 
Now according to (5.49), (5.51) we have 

£i« - T b2 L lP = Ft + (dJ 2 )F 2 + F 3 + F 4 . 
We claim that we have 

L\T b2 p = T b2 Lip - F 5 



with 



ll^i, Oll^^-Jj < /C(M l5 M 2 ){iV( ) + TN(T)}. 



To see this we use (5.13) and (2.15). It follows then that we have 

L 1W = L x (y- T b2 p) =F X + (d z fo)F 2 + F 3 + F 4 + F 5 := F 
where \\F(t, -)|| s _3 is bounded by the right hand side of (5.53). 

Using (5.43) and Lemma 5.11 we may then apply to w Proposition 3.22 with a = s— 1 
to conclude the proof of Lemma 5.14 and thus that of Proposition 5.10. □ 

6. Well-posedness of the Cauchy problem 

Here we conclude the proof of Theorem 1.2 about the Cauchy theory for s > l + d/2 
for the system 

( d t rj + G(jj)i> = 0, 



(6.1) 



w+g^-m -- — rT ^ = o. 



We previously proved the uniqueness of solutions (see Theorem 5.1). To complete 
the proof of Theorem 1.2, it remains to prove the existence. We obtain solutions to 
the water waves system as limits of smooth solutions to approximate systems. This 
approach has been detailed in [1], where we considered the problem with surface 
tension. The analysis is actually easier without surface tension. One reason is that 
with surface tension, we needed in [1] to use some mollifiers with various properties 
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(since we need good estimates for commutators with the principal part of the oper- 
ator). Here it is possible to use a simpler regularization of the equations since the 
reduced paradifferential system involves only operator of order less than or equal 
to 1. 



6.1. The infinite depth case. To explain the scheme of the proof, we first 
consider the case without bottom (r = 0). This case is easier since we know, from 
previous results (see Wu [54, 55], Lannes [39], Lindblad [41]), that the Cauchy 
problem is well-posed for smooth initial data. Then, one can obtain the existence of 
smooth approximate solutions in a straightforward way : by smoothing the initial 
data. 

Denote by J £ the usual Friedrichs mollifiers, defined by J £ = j(eD x ) where j G 
Cg°(R d ), < j < 1, is such that 

j(0 = 1 for |£| < 1, j(0 = for |e| > 2. 

Set V'o = Je4>o and r/g = J £ Vo- Then (V'oj^o) e -ff°°(R d ) 2 and the Cauchy problem 
for (6.1) has a unique smooth solution {ip £ , rj e ) defined on some time interval [0, T*). 
Since s > l + d/2, by applying Proposition 4.1 with r = s — d/2, it follows that there 
exists a function J- such that, for all e G (0, 1] and all T <T e , we have 

(6.2) MKT) < F{F{M sfi ) + T^(M|(T))), 

with obvious notations. Then by standard arguments, we infer that the lifespan 
of (r} £ ,ip E ) is bounded from below by a positive time To independent of e and that 
we have uniform estimates on [0,To]. The fact that one can pass to the limit in the 
equations follows from the previous contraction estimates (see (5.1)), which allows us 
to prove that (n £ , tp £ , B £ , V £ ) is a Cauchy sequence (this argument has been explained 
in [1]). Notice that these estimates were proved under the assumption a{t) > ao/2. 
This actually follows from the a priori bound (6.2), (4.19), (4.20) and a bootstrap 
method. Then, it remains to prove that the limit solution has the desired regularity 
properties. Again, this follows from the analysis in [1]. 



6.2. The general case. In the case with a general bottom, to apply the strat- 
egy explained above, the only remaining point is to prove that, for smooth enough 
initial data, the Cauchy problem has a smooth solution. Namely, for our purposes 
it is enough to prove the following weak well-posedness result (where we allow an 
arbitrarily large loss of N derivatives). 

Proposition 6.1. For all sq > there exists N > such that the following result 
holds. Consider an initial data (r]o,ipo) G H So+N (R d ) 2 such that, initially, the Taylor 
coefficient a is bounded from below by ao > 0. Then, there exists T > and a 
solution (r),ip) G C 1 ([0,T]; H s ° (R d )) 2 to (6.1) satisfying (?](0),^(0)) = (r?o,V>o) and 
a(t) > ao/2. 

Proof. We use a parabolic regularization of the equations and seek solutions of 
the Cauchy problem (6.1) as limits of solutions of the following approximate systems: 



(6.3) 



dtrj — G(n)ip = eAn, 

w+m+^m - 5 — — = eA v>, 

. (j],ip)\ t =o = (J £ r)a,Jsipo), 

74 



where e g]0, 1] is a small parameter and 

(6 - 4) B = i+WW ■ 

Write (6.3) under the form 

u = e £tA u + f e £ ('- r ) A A[u{t)} dr, u = (rj, iji) . 
Jo 

To solve this Cauchy problem we use two ingredients. Firstly we have 

ll^ >^IL 2 ([0,T];.H' s + 1 ) — Wf\\n B ' 

together with a dual estimate, and secondly (for s large enough) 

ll-^MIL^O.T];^- 1 ) < v^'^ r (ll' u lll / oo([o,r];H s )) II u \\l°°([0,T\;H s ) ' 
which follows from Theorem 3.18 and the product rule in Sobolev spaces. Then, 
given e > and (770,^0) £ (H s (R d )) 2 with s large enough, by applying the Banach 
fixed point theorem in L°°([0, T];H S ) n L 2 ([0, T];H S+1 ), one gets that, for T small 
enough possibly depending on e, the Cauchy problem for (6.3) has a smooth solution. 
Moreover, the maximal time of existence T e satisfies: either T e = +00 or 

limsup \\(r],ip)(t, -)\\ Hs = +00. 

We have to prove that these solutions exist for a time interval independent of e G 
]0, 1]. To do so, we begin by computing the equations satisfied by B and V = 
Vip - BVrj. 

We first observe that, since G(r])ip = B — V ■ Vr/, we have 

(6.5) dtv + V ■ Vr/ = B + eArj. 
On the other hand, as in (4.55), we have 

(6.6) dtif) + V ■ W = -m + \\V\ 2 - \b 2 + eAif). 
We next write that 

dtV + V • W = (dt + V • V)(Vip - BVrj) 

= V(d t ip + V ■ Vtp) — (d t B + V • VB)Vr] - BX7(d t T] + V • Vr/) + R 

where R = (i?i, . . . , Rd) with 

R k = ~y~] dkVjdjtp + B^2 dkVjdjT}. 
i 3 
Then, it follows from (6.5) and (6.6) that 

dtV+V-VV = v{-gi 1 +^\V\ 2 -]^B 2 +eAil?j-{d t B+V-VB)Vr]+BV{B+eAri)+R. 

Now, observing that R + ^V\V\ 2 = and simplifying, 

d t V + V • W + (g + d t B + V ■ V5)Vr/ = eVA^ - eBVAi] 
which in turn implies that 

d t V + V • W + aC = eAV + 2e(V-B • V)C, 

where Q = Vr/ and 

(6.7) a = g + d t B + V VB- eAB. 
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Now, for s > 5/2 + d/2, from the proof of (4.38) (using Lemma 2.17 instead of 
Lemma 2.16 as in the proof of Lemma 5.6), we obtain that U = V + T^B satisfies 

(6.8) 3 t U + T V -VU + T a C - sAU - eT VB ■V( = F 1 

where 

ll-^i llif (#») - FdKip,!], V,B,a)\\ L c a{HS+1/ 2 xHS+ l/2 xH s xHS><C l/2)). 



Also, as in the proof of Proposition 4.12, we find that 9 = T q Q (where q = \Ja/X) 
satisfies 

d t 9 + T V -V9- T y U - eA9 - e[T q , A]T q -i9 = F 2 , 
where F 2 is estimated by means of Theorem 2.7 and (2.17): 

ll-^2|li«>(H») < ^(\\(^,V, V,B,a,d t a + V ■ Vd)\\ L ^^ H s+i/2 xH s+i/2 xHsxH s xC i/2 xLX )). 
Using the unknown 9, one can rewrite (6.8) as 
(6.9) d t U + T V -VU + T & C - eAU - eT VB • VT q -i9 = F 1 

where F\ = F\ + eTys ■ V(T 9 T„-i — I)Q and it follows from Theorem 2.7 that 
||r V B • V{T q T r x - I)(\\ HS < \\VB\\ Lao Ml'l^M-^iq- 1 ) H\ H s^ . 



Since M 1 3 / / 2 2 (q)M 3/ \ /2 (q- 1 ) < C(\\r,\\ c5/2 , ||a[| c3/2 ) 

ll-^llif {H°) - -^(H V ' B ^ 5 ')\\L^(H s + 1 / 2 xH s + 1 / 2 xH s xH s xC' i / 2 )" 

Now, since the operator 

'A T VB ■ VT g -i 



we obtain 



A + [T q , A]T q - 

is elliptic (being a perturbation of A of order 3/2), we estimate (U, 9) in L^°{H S xH s ) 
by commuting the equation with (I — A) S//2 and performing an L 2 estimate. To do 
so it remains only to prove that 

||(5, d t d + V ■ Vo)|| Z|0(c .3/2 xX< x. ) < F^WO), rj, V, B)\\ i » (£ p+i/2 xfl *+vaxfl*xfl*)) • 

Here this is much easier than in the above analysis. Indeed, since we do not need to 
work with critical regularity, it is enough to prove that there exist s 2 > s\ > 3/2 + so 
with sq > d/2 and s 2 < s (possibly sq <S si <C S2) such that 



(6-10) ||5 - 3 [| £,00 (/pi) < F\\\\'h ^)\\Lf(H s 2xH s '2) 

(6-11) ll^a|| L oo ( ^o) < r\\\fa$)\\i%>{H*xlP2) 

To prove (6.10) we express a — g in terms of 77,^ (as noticed in [38]). To do so, we 
use two observations. Firstly, by definition of B (see (6.4)), we have 

d t B = TT ^- 2 (vd tV ■ W + V • V<W + dtGm) - 7+|v# ^' 
Secondly, we have 

(6.12) fltGfaty = G(r])(dtip - Bd tV ) - div(Vd t r]). 

The formula (6.12) is proved by Lannes for smooth bottoms (see [39]) or infinite 
depth ([38]). In the case with a general bottom considered here, this follows from [4]. 
Then, replacing dtip and dt"r\ by their expressions computed using System (6.3), we 
are in position to express dtB in terms of rj, only. Setting the result thus obtained 
in (6.7) we get an expression of a — g in terms of rj, ip only. The estimate (6.10) then 
follows from the nonlinear estimates in Sobolev spaces (see (2.17) and (2.21)) and 
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repeated uses of the estimate for the Dirichlet-Neumann proved in Theorem 3.18. 
Repeating this reasoning, we get an expression of dta in terms of r/, tp only and 
hence (6.11). 

Up to now, this gives uniform estimates in e for (U,9) in L°°([0, T]; H s x H s ), 
as long as a(t) is bounded from below by ao/2. Moreover, as in the proof of 
Lemma 4.15 and Lemma 4.16, it follows from (6.11) that one can estimate (77, B, V) 
in L°°([0, T]; H s x H s ° x H so ). Now we can recover estimates for the unknowns 
(r],ip,B,V) in L°°([0, T]; H s+ ^ x H s+ ^ x H s x H s ) by a bootstrap argument as in 
the proof of Lemma 4.17. 

Using these uniform estimates, we get by a bootstrap argument that the solutions 
to (6.3) exist on a time interval independent of e (and satisfy a(t) > ao/2 according 
to (6.10) (6.11)) on this time interval uniform estimates. To conclude it remains 
to pass to the limit in (6.3). For this we use again the uniform estimates which 
hold as long as the Taylor coefficient a(t) is bounded from below by ao/2 (and the 
equation to bound the time derivatives) to extract subsequences converging weakly 
and the fact that the Dirichlet-Neumann operator, though nonlocal, passes to the 
limit thanks to Theorem 3.9. Eventually, we use again a bootstrap argument to 
control the Taylor coefficient using (4.19) and (4.20). □ 

7. Three-dimensional waves in a canal 

Let us recall that we consider a fluid domain which at time t is of the form 

Q(t) = {(xi,X2,y) G (0, 1) x R x R : b(x) < y < rj(t, x), x = (x±, x 2 )} , 

for some given function b. Denote by S the free surface and by T the fixed boundary 
of the canal: 

E(t) = {(xi,x 2 ,y)€(0,l) xRxR : y = V {t,x)}, 
and we set T = dQ(t) \ T,(t) (which does not depend on time). We have 

r = r!ur 2 , 

^ Ti = {( Xl ,x 2 ,y) G (0,1) x T x R;b( Xl ,x 2 ) = y} 

T 2 = { (x 1 ,x 2 ,y) G {0, 1} x T x R ; b(x 1 ,x 2 ) < y < rj{xi,x') } 

Denote by n the normal to the boundary V (remark that near Ti,n = (±1,0,0)) 
and denote by v the boundary to the free surface E. 

We consider an initial surface t]q and an initial velocity field in SI satisfying 

(7.2) div^j, v = 0, curl^ v = 0, v | r -n = 0. 

We write vo = an d define as before 

V = V x <j)(x,r] (x)), B = d y (j)(x,r] (x)). 

We will also need to consider v Xl , v X2 and v y the components of the velocity field 
directed respectively in the xi,x 2 and y direction, Let us recall that we assume 

i) there exists a positive constant h such that tjq(x) > b(x) + h 

ii) The Taylor sign condition (1.8) holds initially. 

We begin with the following (elementary but seemingly new) observation: in the case 
of vertical walls, as long as the Taylor sign condition (1.8) is satisfied, for the Cauchy 
problem to be well posed, it is necessary that at the points where the free surface 
and the boundary of the canal meet (£(£) n T), the scalar product between the two 
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normals (to the free surface and to the boundary of the canal) vanishes : v ■ n = 
on SnT, which means that the free surface E necessarily makes a right-angle with 
the rigid walls (see Figure 4). 

Proposition 7.1. Let (rj,v) G C°([0,T];Ci ([0, 1] x R) x C 1 ^)) 6e a solution of 
System (1.5) suc/i £/taf f/ie Taylor coefficient a is continuous and non-vanishing. 
Then the angle between the free surface, S(i) and the boundary of the canal T is a 
right angle: 

Vt G [0,T],Vx G £(t) fir, n-i/(t,x) = 0, 
which is equivalent to 

(7.3) 5 Xl r/(t, zi, a:') Ui=o,i= 0. 

Remark 7.2. Taking into account that 

d Xl ip = d Xl cp(x,rj(x)) + dy(j)(x,n(x))d Xl r](x), 

and using that on the vertical part of the boundary V, the normal velocity vanishes 
which implies 

d Xl 4>(x 1 ,x / , y) = 0, x 1 = 0, 1, b(x) < y < rj(x), 
and consequently v Xl 1^=0,1= 0. 

PROOF of Proposition 7.1. Since near the free surface, n = (±1,0,0) is con- 
stant, and the vector field v-'V is tangent to T, the boundary condition d n (f> = v-n = 
implies that [dtv + (y ■ V X y)v] ■ n = on Ti, and consequently also by continuity on 
the intersection og T± with E. It follows from the Euler equation that 

V x , y P ■ n = on Ti n E = {(aq = 0, l),x 2 ,y = r/(t, x)}. 

On the other hand, by assumption, the pressure is constant on the free surface 
and hence V XtV P is proportional to the normal to E, v. Notice now that the non 
vanishing of the Taylor coefficient o reads d y P and consequently V XtV P 0. 
This implies that ^-n = 0onEnr. □ 



r 




Figure 4. Two-dimensional section of the fluid domain, exhibiting 
the right-angles at the interface E n T 

Here is our result about the Cauchy problem for water waves in a canal. 
Theorem 7.3. Let s G (2,3), s / \, 

?T((0,l)xR) = i7 s+ 5((0,l)xR)x^ s+ 5((0,l)xR)xi7 s ((0,l)xR)x^ s ((0,l)xR), 
and consider initial data rjo,VQ, satisfying 

(loio,4Bo)eH s ((0,l)xR) 

and 

78 



i) d xi rj Ui=o,i= °> 
n) d xi ip 1*1=0,1= 0, 

Hi) the Taylor sign condition, ao(x) > c > Ois satisfied at time t = and J7o(x) > 
+ /i for some positive constant h. 

Then there exists a time T > and a unique solution (j],v) of the system (1.3) (1.4) 
such that 

i) (n, V>, V, B) € C([0, T);fi s ((0, 1) x R)), 

m) i/ie Taylor sign condition is satisfied at time t and rj(t) >b + h/2. 

Remark 7.4. Assumptions (i) and (ii) are according to Proposition 7.1 and Re- 
mark 7.2 necessary to solve the water waves system. In the case of a flat bottom 
(say b(x) = —1) we do not need assumption (iii) which is in this case always satisfied. 
Also, this condition is satisfied under a smallness assumption. Finally, our result 
excludes the case s = | for technical reasons. It would be possible (but unnecessarily 
complicated) to include this case. 

7.1. Reduction to a result for periodic data. In this section, we shall 
denote by T = R/2Z the torus and consider a domain 

n = { (t,x 1 ,x 2 , y) € [0, T] x R x T x R x R : ( Xl ,x 2 ,y) G Sl(t) }, 

such that, for each time t, one has 

Q(t) = {(x, y) eO : b(xi, x 2 ) < y < rj(t, xi, x 2 )} , 

where b is a fixed continuous function, rj is an unknown function and O is a given 
open domain which contains a fixed strip around the free surface 

S = {(t,x,y) S [0, T] x T x R x R : y = r](t,x)}. 

This implies that there exists h > such that, for all t G [0,T], 

(7.4) b(x 1 ,x 2 ) < n(t,xi,x 2 ) - h. 

Following Boussinesq (see [15, page 37]) the strategy of proof is to perform a sym- 
metrization process (following the process which is illustrated on Figure 3 in the 
introduction). The two main points in this section are to reduce the proof of The- 
orem 7.3 to the proof of the following result and to show how in turn Theorem 7.5 
can be proved by a rather straightforward adaptation of the proof of Theorem 1.2. 

Theorem 7.5. Let s > 2 (notice that the space dimension of our free surface is 
d = 2 and consequently s > 1 + d/2). Let 

H S (T x R) = H S+ ^(T x R) x H S+ ^(T x R) x H S (T x R) x H S (T x R), 

and assume that the initial data satisfy 

(5?o,4^Bo)^ s (TxR) 

are periodic with respect to the x\ variable, the Taylor sign condition is satisfied 
at time t = and rjQ > b + h for some positive constant h. Then there exists a 
time T > and a unique solution (rj,v) of the system (1.3) (1.4) such that 

i) ( Vo ,4>o,Vo,B ) e C([0,T]-n s (T x R)), 

ii) the Taylor sign condition is satisfied at time t and n > b + h/2. 

Before proving this result, let us show how it implies Theorem 7.3 by a simple 
reflection procedure. 
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7.2. Reduction to a periodic setting. Without additional assumptions, the 
reflection procedure should yield in general a Lipschitz singularity. However, here 
the possible singularities are weaker according to the above physical observation 
about the right angles at the interface which implies that 

d Xl rj(t, 0, x 2 ) = and d Xl rj(t, 1, x 2 ) = 0. 

For a function v defined on (0,+oo), define v ev and v od to be the even and odd 
extensions of v to (— oo, +oo) defined by 

v ev (y) = 



v(-y), if y < 
v(y) if y > 0. 



-v(-y), if y < 
v(y) if y > 0. 



• Assume that s G (— T/ien i/ie map v i-> t; oc( is continuous from 



(7.5) 

^ od (y) 

We have the following result 
Proposition 7.6. We have 

• Assume that s G (— |, ^) U (^,|). TTten i/te map v i-» ■u e? ' is continuous 
from # s (0,+oo) to F S (R). 

• Assume that s G (§,§) U (|,|) and /ei if s (0,+oo) &e i/ie closed subspace 
H s (0,+oo) of functions whose derivatives vanish at 0. T/ien i/ie map w i— > 

v ev is continuous from H s (0, +oo) to iiP(R). 
l i> 

iF(0,+oo) to iF(R). ' 

• Assume that s G (^, |) U (|, |) and Zei # s (0, +oo) &e i/ie closed subspace 
H s (0, +oo) of functions vanishing at 0. 27ien i/ie map u i-> t; e " is continuous 
from tf s (0,+oo)) to iF(R). 

• Assume that s G (|,|) U (|,§) and let H s (0,+oo) be the closed subspace 
H s (0, +oo) o/ functions vanishing at suc/i that d x u \ x =q= 0. T/ien the 
map v i—)- w 6 '' is continuous from H S (Q, +oo) to iP(R). 

Proof. For s G (— ^|), it is well known that the space Cq°(I) is dense in H S (I) 
for any intervalle I C R, hence 

w G H s (0, +oo) !->■ ifl ?/ >o 

is bounded from iT s (0, +oo) to iP(R) and the result follows for s G (— 3, §)■ For 
s G (5, |), the jump formula gives 

d y v cv = d y v ev + \v ev }(0) 5 y=0 = d y v ev 

where 

[ v cv ](0) = v cv (0+) -v cv (0-) 

denotes the jump of the function v cv at (which exists since s > \ and is equal to 
since the function v ev is even. As a consequence, since 

K V II^(R) < c(\\v cv \\ L 2 + \\d y v cv \\ HS - im ), 

and the result for s G (3, §) follows directly from the result for s G —5, 5), and the 
fact that Cq°([0, +00)) is dense in H S (R). For s G (|, |), another application of the 



jump formula and the fact that by assumption 

[dyV CV ](0)= 2dyV(0) =0, 
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(7.6) 



gives 

d 2 y v cv = oy v + [d y v ev ] (o) ® ^=0 = d y v ev 

and the result for s G (|, |) follows from the result for s £ (3, §)• A last application 
of the jump formula (and the fact that the function v ev is even) gives the result for 
s G (§, X). The proof for the map v 1— > v od is similar. □ 

Let us now for a function u on (0, 1) x R d define the even (resp odd) periodic 
extensions on T x H d , u ev (resp. u od ), by 

v{— xi,x'), if — 1 < x\ < 

v ev (xi,x') = < v(xi,x'), if 1 > xi > 0, 

v{xi-2k,x'). if xi -2A; G (-1,1) 

—v(—xi,x'), if — 1 < x\ < 
f od (xi,x') = < v(x\,x'), if 1 > x\ > 0, 

k u(xi + 2,x'). if a?i - 2fc € (-1, 1) 

Corollary 7.7. • Assume that s G (— |, |) U (|, |). T/ien toe map d h-> u e,; 

is continuous from H s ((0, 1) x R d ) to # S (T x K d ). 

• Assume that s G (|,|)U(|,|) and Ze£ H s ((0, 1) x R d ) 6e £/ze closed subspace 
-fP((0, 1) x R d ) of functions whose derivatives vanish at and 1. T/ien toe 
map v 1 — ^ u e " is continuous from H s ((0, 1) x R d ) to .£P(T x R d ). 

• Assume that s G ( — 5,5)- T/ien toe map v 1— > i> od is continuous from 
H s ((0, 1) x R d ) to iF(T x R d ). 

• Assume that s G (|,|)U(|,|) and Ze£ H s ((0, 1) x R d ) 6e £/ie closed subspace 
H S ((0,1) x R rf ) of functions vanishing at 0. TTien toe map t> 1— > is 
continuous from H s ((0, 1) x R d ) to H S (T x R rf ). 

• Assume that s G (|,I)U(X,|) and Zei H s ((0, 1) x R d ) 6e toe closed subspace 
H s ((0, 1) x R rf ) of functions vanishing at and 1 suc/i that d x u | x= o i i= 0. 
Then the map v 1— > v ev is continuous from H s ((0, 1) x R d ) to H S (T x R d ). 

Indeed, since (£>!» ev = D*,(v ev ), (D%v) od = D%(v od ) the result is clearly a one 
dimensional result and it is enough to prove it for d = 0, in which case it is a direct 
consequence of Proposition 7.6 and a localization argument. 

We can now show how Theorem 7.5 implies Theorem 7.3. Consider initial data 
v o) satisfying the assumptions in Theorem 7.3. The function «o is defined in the 
domain (0, 1) x R, while the function vq is defined on the domain 

f2 = {(%i,%2,y) G R x (0, 1) x R; 6(xi,x 2 ) < y < rjofei, ^2)}- 

Let (rjo,b) = (77^ , &q v ) be the even periodic extensions of {r]o,b) to T x R defined 
as in (7.6), and V = (v 0}X1 ,v 0jX2 ,v 0} y) = (vQ d Xl ,v^ X2 ,v^ y ). Now v Q is defined on the 
domain 

Q = {(xi,x 2 ,y) G R x T x R;6 (xi,x 2 ) < V < Vo(xi,x 2 )}. 
Notice that the boundary of f^o is 

(7.7) ^o = S ur^, r^ = {(x 1 ,x 2 ,y)GRxTxR;6 (x 1 ,x 2 )=y} 

while the boundary of 0,q is given by (7.1). Denote by 

S = {(xi,x 2 ,y) G R x T x R;y = rj }. 
Let us finally define Vo, -Bo the traces of the horizontal and vertical velocities on So- 
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Clearly, we have 

V , Xl = V ° X1 , V 0jX > = V^ x ,, B = B% v . 

Let us now consider the velocity potential (fro, ifro its trace on Eo, <fro its even periodic 
extension to £lo, and ipo the even periodic extension on tfro on T x R. Notice that 
the function <fro is harmonic in Qq. Indeed, according to the jump formula 

A^ = [d Xl $o] <8> 6fr, 

where [d xl (fro] = 2<9</>o denotes the jump of the function d Xl cfro through T2. But since 
d n <fr \dn = 0) this jump is vanishing. 

We notice now that the following properties are satisfied for e = 0, 1 (as soon as the 

functions are smooth enough so that the traces make sense) 

(7.8) 

dxA^o) \ Xl =e = d Xl 4>{e,x',7](e,x')) + d y <fr(e, x', rj(e, x'))d Xl i]o(e, x') = 
9 X1 (B Q ) \ X1=£ = d Xl d y <fr(e,x' ,r](e,x')) + dy<fr(e, x', rj(e, x'))d Xl r] (e, x) = 
d xl (V 0)X >) \ xi=e = d Xl d x ><fr(e,x' ,rj(e,x')) 

+ d y d x '(j)(e,x',r](£,x'))d x tr]o(e,x') = 
V 0>X1 \ Xl=e = d Xl (fr(e,x' ,rj(e,x')) = 0, 
d Xl V , Xl \ Xl =s = d^(f)(e,x',ri(e,x')) + dyd^^e, x', r)(e, x'))d Xl rj(e, x') 

+dyd Xl (fr(e, x', r)(e, x'))(d Xl r](e, x')) 2 + d y d Xl cfr(e, x', r?(e, x'))d y d Xl r](e, x') = 0, 

where in the last equality, we used that since <fr is harmonic d Xl <fr = —(d x , + dy)d xi (fr 
Now we deduce from Proposition 7.1, assumptions i) ii) in Theorem 7.3 and Corol- 
lary 7.7 that the extensions we just defined satisfy 

(^o,^o,^o,5 ) G H S+ ^(T x R) x H S+ ^(T x R) x H S (T x R) x H S (T x R). 

Let (77, v) be the solution of the free surface water waves system given by The- 
orem 7.5. Since the initial data (770, V'o, ^o', ^0) are even while Vb )X i is odd, our 
uniqueness result guaranties that the solution satisfies the same symmetry property 
(because if we consider our solution, the function obtained by symmetrization is also 
a solution with same initial data). As a consequence if we define v, rj, P as the trace 
of v, rj, P on (0, 1) x R, we get that they satisfy trivially the system free boundary 
Euler equation 

d t v + v ■ V x , y v + V XtV P = -gey, div x>y v = in 17, 

(7.9) dtv = V 1 + l Vr ?l 2 v ■ v 011 s , 

P = on S, 

and to conclude on the existence point in Theorem 7.3, it only remains to check that 
the "solid wall condition" 

(7.10) vn = o, onr = riur 2 

is satisfied. On T\ it is a straightforward consequence of the condition v-T = 0, while 
on T2 it is simply consequence of the fact that the component of the velocity field 
along xi, v xi is odd and 2-periodic. To prove the uniqueness part in Theorem 7.3, 
starting from a a solution of (7.9), (7.10), on the time interval [—T,T], if we define 
the function v, rj at each time t following the same procedure, we end up with a 
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solution of (1.3), (1.4) in the domain {(t,x,y);t G (—T,T),(x,y) G O(i)}, at the 
same level of regularity. Indeed, the jump formula gives 

d t v + v ■ V XiV v + V x , y P = ~ge y + [v Xl ■ d Xl v] ® 5r 2 = -ge y , 

where in the last equality we used that the component of the velocity field along x\ 
vanishes on IV The uniqueness part in Theorem 7.3 consequently follows from the 
uniqueness part in Theorem 7.5 



7.3. Sketch of proof of Theorem 7.5. The proof of our well posedness result 
in the periodic setting is the same word to word as in the non periodic case, once one 
is able to develop a suitable para-differential calculus and show that all properties 
in Section 2 remain true. It is quite standard that the theory of pseudo-differential 
operators (or more generally para-differential operators) can be developed in the 
framework of manifolds. Here, since we are working in the context of T x R, the 
presentation is particularly simple. The torus is endowed with the following atlas: 
let Cjii = 1 , 2 be defined by 

Kl : (-l,l) mo d2 9 X I-)- X G (-1, 1) 
{>Al) k 2 : (0,2) mod2 Bx^xe (0,2). 

and consider the operators on functions 

KjU(x) = u(kJ 1 (x), Kj t *v(x) = v(Kj(x)). 

Let 1 = Ci 2 (^0 + C2 2 ( x ~) be a partition of unity compatible with this atlas (i.e. 
supp(Ci) C (— 1, l)mod2, and supp(C2_) C (0, 2) mod2 , and £ 2 two functions supported 
respectively in (—1, l) mo d2 ; and (0, 2) mod2 such that C Cj = Cj- We use the classical 
definition of pseudo-differential operators on manifolds (i.e. these operators are the 
operators which in the maps of the atlas are usual pseudo-differential operators). In 
our context, the definition is very simple 

Definition 7.8. Given p € [0, 1] and m G R, T^(T x R) denotes the space of 
locally bounded functions a( T a/,£) on T x R x (R d \ 0), which are C°° with respect 
to £ for £ 7^ and such that, for all a G N rf and all £ 7^ 0, the function (x,x') 1— > 
d^a(x,x' ,£) belongs to W P '°°(T x ~R)and there exists a constant C a such that, 

(7-12) V|f| > \, ||^a(.,fl|| w , )ao(TxR) < C (l + |flr-N. 

We shall use the notation 

(l + |£|)W- m 9?a(t;.,0 



Mq (a) := sup sup sup 

te[o,T] \ a \<f +1+p \£\>l/2 



L°°(TxR) 



Definition 7.9. Given a symbol a G T™(T X R), we define the paradifferential 
operator T a associated to a and acting on functions on T x R by 

(7.13) f a u= K J >(«y(? i )r /e *£. a /s5(^w)^ 

For two functions, we also define the para-product 

T u v= K j,*(^*j(l j ) T K*Cu'^i(i v ) 



=1,2 



uv = T u v + T v u + R(u, v) 
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with 

R(u,v)= K j,*( K j(Cp^(^Ci u ' K i^)))' 
(7.14) uv = T u v + T v u + R(u,v). 

In the sequel, to take benefit from the very simple structure of T, we shall identify 
functions on T x R with 2-periodic functions (with respect with the first variable) 
on R 2 . As a consequence, the definition (7.15) becomes 

(7-15) f a u = E Tfc ( E r <>0 n )' 

ke2Z ,7=1,2 

where Tfc(/)(x) = f(x + k), (j(x) = C -( K 7 ( x )) are smooth functions supported on 
(— 1, 1) and (0,2) respectively such that 

E t *(Ci + C2 2 ) = i, 

fee2Z 

and the functions u, a are simply 2-periodic function (with respect to the first vari- 
able) on R 2 . Notice however that in this context, the H s norms have to be taken 
only on a period (say (—1,1) x R or (0, 2) x R). 

Remark 7.10. Though we will not need it, it is worth noticing that the definition 
above does not depend, modulo operators of orders — oo on the choice of the cut- 
off functions Q and Q. Also, for any smooth function bounded as well as all its 
derivatives, £, since — £ is of order — oo, we get that — £ is also of order — oo. 

Taking benefit of the particularly simple structure of our atlas, it is very easy to 
show that all the properties of the calculus of para-diufferential operators proved on 
R rf remain true in this context (and are actually consequences of these properties), 
under the same assumptions (replacing of course R 2 by T x R). On the more 
general framework of an abstract compact manifold, this is still true but require the 
invariance of para-differential operators by diffeomorphisms, while here we use the 
invariance by translations. 

Proposition 7.11. All the estimates and symbolic calculus of Sections 2.2, 2.3 
and 2.4 remain mutatis mutandi true in the context o/T x R. 

Proof. Estimates (2.4), (2.12) (2.13), (2.14) and (2.23) from the definition of 
the paradifferential operator (7.15) and the fact that these estimates hold for the 
usual para-differential operators on R 2 . Similarly, relations (2.9), (2.10) and (2.11) 
follow from (7.9). Let us prove the composition rules (2.5). We have 

(7.16) f a f bU = E r *( E o'VO-CE r «( E c^o))) 

ke2Z 3=1,2 qe2Z i=l,2 

Let us denote by Q\ the contribution of j = 1 in the sum above. Due to the support 
properties of the functions £1 and r q Q, in the expression defining Qi, the sum is 
restricted to q = for (j,i) = (1,1) and q = 0, — 2 for = (1,2). We also have 
according to (2.5) (on R 2 ), 

\\lxi,T Qb ]\\ H ^ m _ m _ m , +p < KMf(( 3 b) < KMf{b). 

Using that — £ is of order —00, we get, modulo an operator which is bounded 
from H^(T x R) to H^~ m - m ' + P(T x R) by 

^M m (Cia)M™'(06) < KM™{a)M™\b), 
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(7.17) Q 1= Y T k (^ 1 T( ia (T C 2 b ( 1 u + T C 2 b ( 1 u + r_ 2 (r C 2 6 Ci u)) 

fcG2Z 

= E^(ci r c ia (r ( 2 2 2)b ci«) 

fcG2Z y 

Applying again (2.5) (on R 2 ), we get that modulo an operator which is bounded 
from H»(T x R) to #m-™-™'+p(t x R) by 

K(M™{a)M?\b) + M^a)M?{b)) 

Qi = E^(Ci r ClX ( cf+C | +r _ 2C |) a6 Ci«) 

fc£2Z V 7 

Noticing now that on the support of (j, we have (£ 2 + (f r -2Cf) = 1) an d estimating 
similarly the contribution of j = 2 in (7.16), this concludes the proof of (2.5) for 
T x R. The proof of (2.6) follows similar lines. 

The most delicate part is to prove that Lemma 2.16 holds in this context. We shall 
actually prove that it is a consequence of the result on R 2 . Let us bound the L 2 
norm of 



(7.18) M E Q T Qp(j( d t u + E ^((Cl^vCl + C 2 T C2 yC 2 )Vt 

fce2Z j=l,2 <je2Z 

- [ d t E Tfc (E or 0p ciw) 

fc£2Z j=l,2 

+ E ^((Ci^vCi + C2T ?2 vC2)v y ^(E o r c^ 

ge2Z fcG2Z 3=1,2 

The contributions of j = 1,2 in the expression above are estimated separately and 
using the same method. Let us estimate 



(7.19) Yl T k((iT ClP (i[dtu + Y T *{[(i T CxV<i + C 2 T ?2V <2)Vu 

fce2Z <j£2Z 

- [d t Y Tk(CiT ClP Ciu) + Y ^((Ci^vCi + C2T C2 yC 2 )V Y Tk(CiT Cll 
fce2Z qe2Z fee2Z 

We compute the L 2 -norm on ( — 1, 1). Since for k ^ 0, Tfc(d) is disjoint from ( — 1, 1) 
and for k ^ 0,-2, Tk(C,2) is disjoint from (—1,1), while for k ^ 0,2 Tk((i) has its 
support disjoint from £2; we need only to estimate the L 2 -norm of 



(7.20) Ci^pCi [dtu + (CiT^vdVu + C 2 T C2 yC 2 Vn + r_i ((2^(2) Vu y 

- 3 t CiT ?lP Ci^ + (CiTc 1 vCiV(CiTc 1 pCi'") 

+ E r g(C2^C2vC2V( ^ r fc ^Cir Clp Ci^, 

9=0,-2 fc=0,2 

We shall say that a term in the expression above is admissible if its L 2 norm on 
(—1, 1) x R is bounded by 

KM^(p)\\V(t)\\ cl+e \\u(t)\\ Hm{TxR) . 

According to Proposition 2.11, we can replace in (7.20) the multiplication by any 
smooth function £ by in the expression above, modulo admissible remainder 
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terms R. According to (2.5), we can also replace in (7.20) T^T^y by T^y and we 
can commute and T^y modulo remainder terms. As a consequence, we end up 
estimating (using that CjCj = Cj an d that t^V = V, t^u = u) 



(7.21) CiT ClP (a t <>- ( r ( f?+f | +T _ 2C| )v V ^ 



C: 



d t T ClP Ciu - [T ( 2 V V o T Clp ( lU + Yj T r 9 Cl^ Vo Yl T T q +kC lP T q +k{(;i)u 

g=0,-2 fc=0,2 



Notice here that in the sum above we can keep only the terms such that q + k = 0, 
because modulo admissible remainder terms we have 

(7.22) Tq^T^vCiV o TkfcT^pCxu)) = T q+k CiT q (c 2 T( 2 vC2V ° T fc (T Cl pCin 
which is vanishing on (—1,1) ifg + Zc^O. Finally, we end up with 



(7.23) CiT (lP [dtCiu+ (T (Ci2+c|+r _ aCf)v VCiu 



- Ci ($T flP Ci« + ( r ( cf +C | +T _ 2Cf ) v V r Ci P Ci« 
= C 



T( 1P , d t + ( T, 2+c|+r _ 2C| ) v V o T ClP 



Citt 



We can now apply Lemma 2.16 in the context of R 2 to get that the L 2 norm of this 
expression is bounded by 



(7.24) K(MZ(( lP ) || (Ci 2 + Cl + T„ 2 (l)V(t)\\ ci+£ 
+ Mtffdtihp) + (Ci 2 + C 2 2 + r_ 2 C 2 2 )^ • V(Cip)) || (Ci 2 + C 2 2 + r_ 2 c|)y(t)|| Lo 

IICl«(*)llfl«(Ra) • 

We now use that on the support of Ci, (Ci +C 2 + r - 2 C 2 )^ = V an d obtain that (7.24) 
is bounded by 

(7.25) k{m^( p ) \\V{t)\\ cl+e + M^fdtp + V-Vp) \\V(t)\\ Loo ) \\u(t)\\ Hm{Txn) . 
This concludes the proof of Proposition 7.11. □ 

Remark 7.12. Let us remark that the method developed in this section (including 
the geometric result in Proposition 7.1) could be easy adapted to deal with the 
totally periodic framework T x T (the case of wave pools) 
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